AFIT/GGC/EE/7  9-5 


USING  VERTICAL  GUST  ALLEVIATION  TO 
IMPROVE  THE  TARGET  TRACKING  CAPABILITY 
OF  THE  CONTROL  CONFIGURED  YF-16 


THESIS 

AFIT/GGC/EE/79-5  Eric  R.  Molner 

1  Lt  US  AF 


-D'D  C 

mpaiP.noM 

'W  RB  11  WO  ■  j 


irTTT-'TL'  \J 

A 


Approved  for  public  release;  distribution  unlimited 


AFIT/GGC/EE/79-5 


USING  VERTICAL  GUST  ALLEVIATION  TO 


IMPROVE  THE  TARGET  TRACKING  CAPABILITY 
OF  THE  CONTROL  CONFIGURED  YF-16  * 


Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 


in  Partial  Fulfillment  of  the 


NIIS  GRUI 
EE>C  I  A3 
Unannouncrd 
Ju~t if i c.  ■ : 


Eric  R./Molner/  flpfe 
1st  Lieutenant  USAF 


Approved  for  public  release;  distribution  unlimited 


Pre  face 


This  thesis  was  undertaken  to  investigate  the  possibility 
of  designing  a  flight  control  system  for  the  YF-16  that  is 
capable  of  reducing  aircraft  pitch  response  to  unwanted 
disturbances,  thereby  improving  target  tracking  performance. 
It  was  realized  at  the  start  of  this  effort  that  a  system  so 
designed  generally  tends  to  provide  underdamped  responses  to 

pilot  commands  as  well,  resulting  in  degraded  overall 

/ 

handling  qualities.  It  was  anticipated  that  this  problem 
could  be  eliminated  by  augmenting  the  control  afforded  by  the 
horizontal  tail  with  a  second  longitudinal  aerodynamic 
control  surface.  This  was,  indeed,  a  workable  approach.  As  is 
demonstrated,  coordinated  deflection  of  both  the  horizontal 
tail  and  the  trailing-edge  flaperons  can  provide  both  gust 
alleviation  and  improved  aircraft  response  to  pilot  commands. 

I  would  like  to  thank  my  advisor,  Dr.  C.  H.  Houpis,  for 
his  continued  patience  and  guidance  from  the  inception  to  the 
completion  of  this  thesis.  I  am  also  greatly  indebted  to  my 
sponsor,  Kr .  Finley  Barfield,  who  suggested  the  project,  and 
who  was  never  too  busy  to  take  the  time  to  help  me  find  some 
of  the  answers.  I  must  also  thank  my  wife,  Debra,  for  her 
unceasing  and  much  appreciated  love,  understanding,  assis¬ 
tance,  and  tolerance  during  the  many  months  that  my  attention 
was  focused  on  this  project. 
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Abstract 


V 

This  thesis  investigates  the  use  of  multiple  aerodynamic 
control  surfaces  as  a  part  of  a  flight  control  system  con¬ 
figured  to  reduce  aircraft  pitch  response  to  atmospheric 
disturbances.  Coordinated  deflection  of  two  surfaces  provides 
the  desired  gust  alleviation  with  no  decrease  in  system, 
damping.  At  the  same  time,  the  ccmmanded  aircraft  performance 

is  improved  relative  to  the  performance  of  the  baseline 

/ 

flight  control  system  using  a  single  control  surface. 

The  aircraft  chosen  for  this  analysis  is  the  General 
Dynamics  YF-16.  This  particular  aircraft  has  a  fly-by-wire 
electronic  flight  control  system  -  an  arrangement  that 
permits  the  control  surfaces  tc  be  electronically  activated. 

A  means  of  applying  discrete  vertical  gusting  to  a 
modeled  system  is  presented  to  permit  evaluation  of  system, 
performance.  Using  a  simulation  computer  program,  the  YF-16 
with  baseline  flight  control  system  and  Maneuver  Enhancement 
augmentation  are  examined  to  determine  gust  response  and 
pilot  response  characteristics. -With  this  analysis  completed, 
a  basis  for  comparison  has  been  established.  To  begin  the 
design  process,  a  system  configuration  is  selected,  and  its 
equation  derived.  This  matrix  equation  is  used  to  solve  for 
transfer  functions  in  terms  of  unspecified  system  variables. 
With  the  aid  of  a  second  computer  program,  two  design 
solutions  are  derived  in  an  iterative  manner,  each  of  which 
satisfies  the  design  requirements.  Simulations  of  these 
systems  are  run  to  determine  actual  performance  relative  to 
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the  defined  target  tracking  criteria  and  the  current  YF-16 
flight  control  systems. 


USING  VERTICAL  GUST  ALLEVIATION  TO 
IMPROVE  THE  TARGET  TRACKING  CAPABILITY 
OF  THE  CONTROL  CONFIGURED  YF-16 


I .  I  n troduct  ion 


Background 

In  January  of  1975,  a  full-scale  development  contract  was 
awarded  to  General  Dynamics  to  produce  the  first  lightweight 
fighter  aircraft  in  the  Air  Force  inventory  -  the  YF-16  (Ref 
3).  With  a  length  of  just  over  46  feet,  a  wingspan  of  31 
feet,  and  a  combat  weight  of  less  than  16,600  pounds,  the 
YF-16  is  considerably  lighter  and  substantially  smaller  than 
most  present-day  fighters.  As  an  example,  it  is  only  half  the 
weight  of  the  F-4E  which  it  is  destined  to  replace  (Ref  17). 
In  addition,  with  its  single  Pratt  &  Whitney  F100-PW-100 
afterburning  turbofan  engine  producing  25,000  pounds  of 
thrust,  the  aircraft  has  a  power-to-weight  ratio  greater  than 
unity.  This  results  in  incredible  acceleration  and  takeoff 
rolls  as  short  as  1000  feet  with  full  afterburner.  The  YF-16 
has  a  maximum  speed  of  Mach  2+  at  40,000  feet  (Ref  21). 

One  of  the  reasons  for  the  YF-16 's  superior  performance 
is  its  light  weight.  Its  supersonic  performance,  however,  is 
due  to  what  has  been  called  "relaxed  longitudinal  static 
stability".  This  simply  means  that  the  aircraft  has  been 
designed  to  be  unstable  in  sub-sonic  flight.  This  situation 
is  created  by  locating  the  center  of  gravity  (CC)  aft  of  the 
aerodynamic  center  (AC).  Traditionally*,  aircraft  have  been 


designed  to  be  aerodynamically  stable  in  the  longitudinal 
mode.  This  provides  not  only  a  trim  condition,  but  tends  to 
bring  the  aircraft  back  to  straight  and  level  attitude  when 
disturbed  from  a  trimmed  flight  condition.  Were  it  not  for 
this  feature,  the  pilot  would  continuously  be  adjusting  and 
compensating  to  keep  the  aircraft  flying;  in  many  instances, 
the  aircraft  could  not  be  flown  at  all.  As  an  aircraft  breaks 

*  0 

through  the  transonic  region  (Mach  0.9  to  Mach  1.1)  and  goes 
supersonic,  however,  it  becomes  "super  stable",  as  the  AC 
effectively  shifts  rearward.  Stability  is  increased  to  the 
extent  that  maneuvering  becomes  difficult.  Generally,  only 
slow  and  gradual  maneuvers  are  possible  -  due,  in  part,  to 
the  supersonic  flow  of  air  over  the  control  surfaces. 

In  order  to  make  an  aircraft  maneuverable  at  supersonic 
speeds,  the  longitudinal  stability  margin  must  be  reduced.  No 
net  improvement  can  be  obtained,  however,  without  making  the 
aircraft  unstable  in  sub-sonic  flight  -  a  situation  which, 
until  recently,  was  not  tolerable.  Advances  in  automatic 
control  system  technology  have  provided  a  solution  to  this 
problem.  The  answer  is  the  electronic  fly-by-wire  (FBK) 
flight  control  system.  Instead  of  a  complex  network  of 
mechanical  linkages,  this  system  uses  wires  to  relay  pilot 
‘  commands.  Using  various  transducers  and  gyros,  all  pilot  and 

aircraft  control  information  is  made  available  in  the  form  of 
electrical  signals.  With  such  a  system,  it  is  a  relatively 
simple  task  to  design  a  control  system  using  electronic 
circuits  to  provide  compensation  and  closed-loop  control. 
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The  YF-16  uses  such  a  control  system  for  primary  and 
secondary  control.  With  its  basic  airframe  designed  to  be 
unstable  for  sub-sonic  flight,  the  aircraft  relies  upon  the 
control  system  for  stabilization.  The  advantage  is  that  in 
sub-sonic  flight,  the  aircraft  can  be  made  to  fly  in  a  normal 
(stable)  manner,  and  in  super-sonic  flight,  the  stability  can 
be  electronically  relaxed  to  allow  greater  maneuverability. 
In  addition,  overall  drag  is  reduced  as  a  result  of  the 
unstable  airframe  which  reduces  fuel  requirements. 

The  YF-16  is  clearly  a  technologically  advanced  aircraft. 
Because  of  this,  it  was  selected  over  the  F-4  to  be  used  in 
conducting  further  research  in  the  area  of  Control  Configured 
Vehicles  (CCVs) .  A  CCV  is  a  vehicle  in  which  advanced  control 
technology  as  well  as  aerodynamics,  structures,  and  propul¬ 
sion  are  all  employed  in  the  initial  configuration  definition 
process  (Ref  10).  To  the  YF-16,  this  has  meant  the  addition 
of  vertical  canards  mounted  to  the  underside  of  the  inlet 
structure,  isolated  fuel  tanks  for  manual  control  of  the  CG, 
and  control  system  modifications  and  additions.  Under  a  $6 
million  contract  from  the  Air  Force,  General  Dynamics  has 
modified  the  No.  1  prototype  YF-16  for  this  purpose. 

There  are  several  potential  advantages  to  the  YF-16  CCV 
concept.  As  it  incorporates  multiple  control  surfaces  about 
the  longitudinal  and  lateral-directional  axes,  it  is  capable 
of  maneuvers  that  most  other  fighter  aircraft  cannot  achieve. 
In  the  longitudinal  mode,  the  control  system  uses  flaperons 
(combination  of  trailing-edge  flaps  and  ‘  ailerons)  and  a 


YF-16  Control  Configured  Vehicle. 


horizontal  tail  to  enable  direct  lift  control,  fuselage  pitch 
pointing,  and  vertical  translation.  In  the  lateral-direction¬ 
al  mode,  the  rudder  and  vertical  canards  are  used  to  enable 
direct  side  force,  fuselage  yaw  pointing,  and  lateral  trans¬ 
lation.  It  has  been  shown  that  in  a  combat  situation,  these 
features  serve  to  improve  tracking  ability  and  increase 
firing  opportunities.  Figure  1  shows  the  CCV  configuration. 

An  additional  mode  offered  by  the  CCV  control  system  is 
known  as  Maneuver  Enhancement  (ME).  In  this  longitudinal 
mode,  the  trailing-e dge  flaperons  are  controlled  by  the 
difference  between  pilot  G-command  and  aircraft  normal 
acceleration  response.  This  results  in  faster  pitch  rate  and 
normal  acceleration  transient  response,  and  provides  a  reduc¬ 
tion  in  aircraft  response  to  gust  disturbances.  It  has  been 
determined  through  flight  testing  that  the  ME  augmentation 
makes  it  noticeably  easier  to  perform  precision  tasks  such  as 
a  constant-G  turn  at  constant  altitude.  The  end  result  is 
increased  maneuverability  and  improved  tracking  ability. 

Problem  Statement 

Because  the  YF-16  is  a  small,  light  weight  aircraft,  it 
is  more  susceptible  to  the  perturbing  effects  of  wind  gusting 
than  most  other  fighters.  For  target  tracking  and  weapon 
delivery,  it  is  essential  to  minimize  the  effect  of  any 
external  disturbances.  General  Dynamics  has  shown  that  by 
feeding  back  normal  acceleration  to  the  flaperons,  the  amount 
of  unwanted  G-response  at  the  pilot  station  as  a  result  of 
vertical  gusting  is  substantially  reduced.  This  type  of 


feedback  does  little,  however,  to- reduce  the  perturbations  in 
pitch  that  also  occur.  Reduction  of  normal  acceleration 
response  provides  mostly  ride  control.  For  the  tracking  task, 
perturbations  in  pitch  angle  generally  have  a  more  disturbing 
effect  on  the  aircraft. 

The  problem  addressed  in  this  study  is  the  development  of 
a  flight  control  system  that  improves  the  tracking  perfor¬ 
mance  and  disturbance  rejection  properties  of  the  aircraft, 
without  degrading  overall  /flying  qualities.  Disturbance 
induced  perturbations  in  pitch  angle  are  of  primary  interest. 

Order  of  Presentation 

When  setting  out  to  design  an  aircraft  flight  control 
system,  there  are  many  prerequisites  that  must  be  satisfied 
prior  to  starting  the  actual  design  phase.  To  begin  with, 
regardless  of  the  control  system  configuration,  the  objective 
is  to  somehow  modify  the  open-locp  dynamics  of  the  aircraft 
to  provide  the  required  performance  improvement.  To  accom¬ 
plish  this  task,  some  means  of  mathmatically  modeling  the 
physical  aircraft  must  be  determined.  The  required  non-dimen¬ 
sional  longitudinal  equations  of  motion  are  presented  in 
Chapter  II.  Stability  derivatives  for  these  equations, 
derived  from  wind  tunnel  data  for  various  flight  conditions, 
are  given  in  Appendix  B. 

As  the  flight  control  system  to  be  designed  is  intended 
to  reduce  the  effect  of  a  given  disturbance  on  the  aircraft, 
a  means  of  applying  that  disturbance  to  the  aircraft  model 

•4 

must  be  determined,  as  well  as  the  shape  or  nature  of  the 
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disturbance  itself.  Also,  some  means  of  measuring  the  effect 
of  the  unwanted  disturbance  on  the  desired  performance  of  the 
aircraft  must  be  established  to  permit  evaluation  of  the 
flight  control  system.  These  requirements  are  also  satisfied 
in  Chapter  II. 

An  analysis  of  the  YF-16  baseline  and  Maneuver  Enhance¬ 
ment  flight  control  systems  is  presented  in  Chapter  III.  This 
is  done  to  provide  an  assessment  of  the  performance  of  these 
systems  and  to  establish  a  b,asis  for  comparison.  Both  time 
response  and  frequency  response  data  is  given,  all  of  which 
is  obtained  through  the  use  of  a  digital  computer  program, 
developed  as  a  part  of  this  thesis.  A  description  of  the 
features  and  capabilities  of  this  program,  as  well  as  an 
example  of  its  use,  is  given  in  Appendix  K. 

The  design  and  analysis  of  the  desired  flight  control 
system  is  presented  in  Chapter  IV.  The  chapter  begins  by 
restating  the  design  objective.  This  objective  is  expanded  to 
include  the  elimination  of  several  undesirable  characteris¬ 
tics  of  the  current  YF-16  flight  control  systems,  as  deter¬ 
mined  in  the  previous  chapter.  Also  addressed  are  the  sub¬ 
jects  of  design  approach  alternatives  and  the  limitations  of 
the  rigid-body  assumption. 

The  actual  design  process  begins  with  the  selection  of  a 
system  configuration.  Alternate  configurations,  considered 
but  determined  to  be  less  appropriate  in  light  of  the  design 
criteria,  are  presented  in  Appendix  F.  A  matrix  equation  is 
derived  for  the  selected  system  configuration,  «nd  is  used  to 
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determine  system  transfer  functions.  These  transfer  functions 
are  used  in  conjunction  with  a  second  computer  program, 
described  in  Appendix  L,  to  determine  two  design  solutions. 
The  chapter  concludes  with  an  evaluation  of  the  performance 
of  those  design  solutions.  This  is  accomplished  through 
simulation  and  comparison  to  the  performance  of  the  current 
YF-16  flight  control  systems. 

Lastly,  conclusions  drawn  from  this  effort  and  recommend¬ 
ations  for  continued  research  ,are  presented  in  Chapter  V. 


( 


II.  Modeling  and  Simulation  Techniques 

Chapter  Overview 

In  this  chapter,  a  set  of  linear  equations  of  irotion  is 
presented  along  with  the  assumptions  used  in  their  deriva¬ 
tion.  The  use  of  angle-of-attack  perturbation  for  simulation 
of  vertical  gusting  is  discussed,  as  is  the  technique  for 

using  it  in  conjunction  with  the  particular  equations  of 

/ 

motion  that  are  chosen.  The  modeling  of  air  turbulence  as  a 
discrete  "1-cos"  gust  is  then  addressed.  Finally,  longitu¬ 
dinal  target  tracking  error  is  defined  as  the  means  of 
evaluating  disturbance  rejection  performance  in  this  thesis. 

Longitudinal  Equations  of  Motion 

In  order  to  investigate  how  various  types  of  control 
affect  the  disturbance  rejection  properties  of  a  given 
aircraft  and  flight  control  system  configuration,  a  modeling 
scheme  must  be  chosen  to  replicate  the  aircraft  dynamics. 
This  model  must  be  simple  enough  to  allow  its  use  in  design 
and  analysis  work,  but  complex  enough  to  provide  a  suffi¬ 
ciently  realistic  simulation  of  the  aircraft  dynamics  for  the 
intended  purpose.  By  far,  the  most  popular  models  for  analyt¬ 
ic  design  are  linear.  A  linear  model  allows  the  designer  to 
use  the  myriad  of  linear  systems  design  techniques  developed 
over  the  years,  greatly  simplifying  the  design  process.  In 
addition  to  assuming  linearity,  many  modeling  schemes  neglect 
flexible  body  modes,  and  perform  a  so-cal  lec.  r  igid-body 
analysis.  For  a  longitudinal  study  of  high  performance  jet 
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fighter  aircraft,  this  simplification  is  generally  quite 
valid  . 

As  the  dynamic  characteristics  of  an  aircraft  are 
actually  highly  non-linear,  a  linear  aircraft  model  is  only 
acccurate  and  valid  over  a  small  operating  range  about  an 
initial  condition.  This  method  of  modeling  and  analysis, 
known  as  perturbation  analysis,  usually  provides  sufficient 
accuracy  for  most  design  purposes  if  the  perturbations  from 
the  initial  (trimned)  flight  /condition  are  kept  "reasonably" 
small.  In  most  cases,  depending  upon  the  particular  aircraft 
and  flight  condition  in  question,  reasonably  small  means 
that  angular  perturbations  are  restricted  to  a  few  degrees  or 
less. 


For  longitudinal  analysis,  as  in  this  report,  a  three- 
degree-of-f r eedom  linear  model  is  generally  used.  Depending 
upon  the  intent  of  the  analysis,  the  aircraft  can  be  modeled 
relative  to  any  one  of  several  reference  frames,  such  as 
stability-axis  or  wind-axis  systems  (Ref  7).  Stability-axis 
equations  are  chosen  for  this  study.  This  set  of  equations 
models  the  longitudinal  motion  of  the  aircraft  relative  to 
the  earth,  where  the  earth  is  considered  to  be  an  inertial 
reference.  The  two  standard,  and  equally  valid,  forms  of  this 
model  are  dimensional  and  non-dimensional.  The  choice  cf 
which  to  use  is  usually  made  according  to  what  type  of 
aircraft  dat*  is  available,  as  one  form  may  be  more  directly 
applicable.  Both  forms  produce  the  same  model  dynamics,  as 


only  the v  gain  constants  are 


affected  to  account  for  the 


dimensions  of  the  states.  Thus,  the  modeling  technique  chosen 
for  this  report  is  a  thr ee-degree-of-f reedom,  non-dimensional 
stability-axis  set  of  equations  (Ref  1).  As  the  derivation  of 
these  equations  is  somewhat  lengthy  and  involved,  no  attempt 
is  made  to  derive  their,  here.  There  are  excellent  books  on  the 
subject  (Refs  1,  7) ,  and  the  reader  unfamiliar  with  this 

modeling  practice  is  encouraged  to  consult  the  bibliography 
of  this  thesis. 

As  these  linear  equation^  attempt  to  model  a  non-linear 
system,  several  assumptions  are  made  in  their  derivation.  It 
is  important  to  understand  these  assumptions  to  ensure  that 
the  model  is  valid  for  a  given  analysis,  and  that  the  limita¬ 
tions  on  the  accuracy  of  that  analysis  are  known.  These 
assumptions  are  listed  below  (Ref  1): 

1.  The  X  and  Z  (stability)  axes  lie  in  the  plane  of 
symmetry,  and  the  origin  of  the  axis  system  is  at  the 
center  of  gravity  (CG)  of  the  aircraft. 

2.  The  mass  of  the  aircraft  is  constant  (fuel  consump¬ 
tion  is  not  considered) . 

3.  The  aircraft  is  a  rigid  body  (hence  the  term 
"rigid-body  analysis"). 

4.  The  earth  is  an  inertial  reference. 

5.  The  perturbations  from  equilibrium  are  small. 

6.  The  flow  of  air  is  quasi-steady. 

In  addition,  the  equations  require  that  the  X-axis  be  aligned 
with  the  aircraft  velocity  vector  while  the  vehicle  is  in 
equilibrium  flight. 


These 


equations,  in  terms  of  stability  derivatives  and 
constants  (the  values  of  which  depend  upon  flight  condition 
and  aircraft  structure)  are  given  below  in  Egns  (1),  (2),  and 

(3)  (Ref  1) .  The  non-dimensional  stability  derivatives  are 
determined  and  listed  for  various  flight  conditions  in 
Appendix  B. 


'  u  -  Cxu  '  uj  -  Cx£ '  Ot  Cw  (  cos  0O  )  0  ^Fxa  ^  ^ 

-czu‘u  +  |j|^  -  -^jCZ(iJ  ’a  -  cZa’aj 

+  -  Cw(sineo)eJ=  cFza  (2) 

[■|jCind,6t  “  Cina,a]  +  [jgqc®  ■  =  Cma  (3) 


If  the  initial  conditions  are  assumed  to  be  zero,  as  is  the 
case  for  the  perturbation  variables  'u,  'a/  a™3  0»  these 

equations  can  be  written  directly  in  Laplace  notation  as 
follows: 

s  -  CxuJ  'u(s)  -  Cxjj'a(s)  -  Cw(  cos  0O  )  9(s)  -  CFxa  (4) 

-0^  'u(s)  ♦  |j5i  -  |,c2cij  S  -  c^j  a(s) 

+  [[-^  -  5  -  °»(  Si"  >]8(8)  =  °F^a  <5) 

[‘Sja">as  -  Sjl  '«(»)  +[^B'  '  ^C"5S]e(5)  =  c"a  ,6) 
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Figure  2.  Equilibrium  and  disturbed  aircraft  longitudinal  stability  axes  (Ref  l). 


The  definition  of  the  axis  system  and  the  perturbation 
velocities  and  angles  is  shown  in  Figure  2.  In  the  figure, 
XE,  YE,  and  ZE  define  the  earth  reference  axes,  XG,  Y0 ,  and 
ZD  define  the  equilibrium  aircraft  axes  (orientation  of  the 
aircraft  in  trimmed  flight),  and  X,  Y,  and  Z  define  the 
disturbed  aircraft  axes.  The  Y  axis  is  defined  out  of  the 
right  wing  of  the  aircraft,  and  for  a  longitudinal  analysis, 
Y  =  Y0  =  Ye.  The  X  axis  of  the  aircraft  is  initially  aligned 
with  the  Xq  axis  (into  the  Relative  wind),  after  which  it 
remains  fixed  to  the  aircraft  during  the  study  of  perturba¬ 
tions  from  the  initial  flight  condition  defined  by  the  Xc 
axis. 

The  aircraft  velocities  of  interest  in  a  longitudinal 
analysis  are  the  linear  forward  velocity,  U,  the  linear 
downward  velocity,  W,  and  the  angular  pitching  velocity,  C* 
In  a  perturbation  analysis,  these  velocities  are  observed  as 
they  change  from  an  initial  steady  state  condition.  As  it  is 
only  the  changes  that  are  of  interest,  these  velocities  can 
be  expressed  as  the  sum  of  the  initial  condition  and  the 
change  from  that  condition,  such  that 

U  =  UG  +  u 

W  =  WD  +  w 

Q  =  Qo  +  q  (7) 

where  UD ,  WQ ,  and  Qc  are  the  equilibrium  velocities,  and  u, 
w,  and  q  are  the  changes  in  these  velocities  resulting  from 
some  d isturbance .j Once  the  initial  flight  condition  has  been 
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established,  U0 ,  W0 ,  and  Co  are  taken  as  constants.  Because  a 
tr  inured  aircraft  is  in  unaccelerated  flight,  the  equilibr  ium 
pitching  velocity,  Q0 ,  and  the  equilibrium  downward  velocity, 
W0 ,  are  equal  to  zero;  thus  0  =  q,  and  W  =  w.  There  is, 
however,  a  non-zero  equilibrium  forward  velocity,  UD .  If  the 
perturbation  forward  velocity,  u,  is  divided  by  this  equili¬ 
brium  forward  velocity,  the  non-dimensional  perturbation 
forward  velocity,  'u,  results. 

The  angles  of  interest  are  the  pitch  angle,  0,  the  attack 
angle,  a,  and  in  some  instances  the  flight  path  angle,  y.  The 
initial  flight  path  angle,  y0  •  may  he  defined  as  the  angle 
between  the  horizontal  and  the  equilibrium  velocity  vector  of 
the  aircraft,  VTC •  Note  that  in  using  stability  axes,  the 
initial  pitch  angle,  0O ,  is  equal  to  the  initial  flight  path 
angle.  As  with  the  velocities  previously  described,  the  total 
pitch  angle  is  equal  to  the  sum  of  the  initial  pitch  angle 
and  the  perturbation  pitch  angle,  0. 

The  angle-of-at tack ,  a,  is  defined  as  the  angle  between 
the  velocity  vector  of  the  aircraft  and  the  chord  of  the 
wing.  This  angle  is  the  sum  of  the  trimmed  attack  angle,  (the 
attack  angle  in  initial  equilibrium  flight),  and  the  pertur¬ 
bation  attack  angle,  'a.  From  Figure  2,  this  perturbation 
attack  angle  is  defined  as 

'a  =  tan"*  (  w/  U  )  (8) 

The  ratio  of  the  two  velocities  is  a  dimensionless  quantity, 
or  radian  measure.  Since  an  underlying  assumption  of  pcrtur- 


bation  analysis  is  that  the  variation  in  velocities  and 
angles  is  kept  small,  this  definition  can  be  approximated 
with  sufficient  accuracy  as 

•a  =  w/U  (9) 

It  follows,  then,  that  the  time  rate  of  change  of  'a  is 
defined  approximately  by 

'a  =  w/  U  (10) 

/ 

Control  Surface  Inputs 

In  Eqns  (4),  (5),  and  (6),  the  inputs  are  defined  as 

CFxa  ,  CpZa  ,  and  ,  respectively.  Cpx  represents  the  effect 
of  the  sum  of  the  forces  in  the  X  direction  resulting  from 
the  deflections  of  all  control  surface  inputs.  In  this 
analysis,  these  inputs  are  the  horizontal  tail  and  the 
tr ailing-edge  flaperons.  The  force  created  by  these  surfaces 
in  the  X  direction  is  a  change  in  total  drag.  Similarly,  Cp^ 
represents  the  effect  of  the  sum  of  the  forces  in  the  Z 
direction  resulting  from  the  deflection  of  the  control 
surfaces.  This  total  force  is  a  change  in  lift.  C,^  ,  the 
input  to  the  remaining  equation,  represents  the  effect  of  the 
movement  of  the  control  surfaces  on  the  pitching  moment. 

For  an  aircraft  with  both  horizontal  tail  and  trailing- 
edge  flaperon  control  surfaces,  the  above  input  terms  are 
defined  as 


CFxa  =  Cx6ht6ht  +  Cx6tef  6tef 
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CFza  =  Cz6ht6^t  +  Cz&tef6tef  (12> 

cma  =  Cm6ht 6ht  +  Cin6tef &tef  (13* 

where  6ht  and  &tef  are  the  actual  angular  deflections  of  the 
control  surfaces.  There  are  certain  sign  conventions  associ¬ 
ated  with  these  angular  deflections,  and  they  are  depicted  in 
Figure  3.  Note  that  a  positive  deflection  produces  a  positive 
change  in  lift  on  the  deflecting  surface.  A  positive  deflec¬ 
tion  of  the  horizontal  tail  creates  lift  on  the  tail,  causing 
it  to  move  upward  relative  to  the  fuselage.  The  result  is  a 
negative  change  in  pitch  angle,  and  a  nose  down  attitude. 
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causing  the  aircraft  to  lose  altitude.  A  positive  deflection 
of  the  trailing-edge  flaperons  creates  lift  on  the  wing, 
causing  the  aircraft  to  gain  altitude.  There  may  also  be  a 
change  in  pitching  moment,  depending  on  the  size  and  location 
of  the  flaperons. 

Simulation  of  Vertical  Turbulence 

Atmospheric  disturbances  occur  in  a  random  fashion  as 

quantities  of  air  move  about  through  space.  In  a  three 
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dimensional  world,  the  velocity  of  this  moving  air  can  be 
thought  of  as  the  sum  of  three  velocity  components  relative 
to  some  axis  system.  In  the  analysis  of  aircraft  control, 
these  velocity  components  are  defined  as  vertical  gusting, 
longitudinal  gusting,  and  lateral  gusting,  relative  to  the 
aircraft  axis  system. 

For  the  longitudinal  analysis  performed  in  this  report, 
it  is  assumed  that  the  lateral  gust  can  be  ignored,  as  it  is, 
by  the  nature  of  the  aircraft  axis  system,  perpendicular  to 
the  vertical  plane  (wings-level  flight) .  In  addition,  as  the 
small  change  in  the  relative  forward  velocity  of  the  aircraft 
caused  by  a  longitudinal  gust  results  in  an  assumed  negli¬ 
gible  change  in  lift,  this  gust  velocity  component  is  not 
considered.  The  vertical  gust,  however,  has  a  direct  effect 
on  total  lift,  and  is  the  gust  velocity  component  used  for 
the  gust  response  analyses  in  this  thesis. 

A  vertical  gust  causes  a  change  in  the  net  velocity  of 
the  air  surrounding  the  aircraft.  An  upward  gust,  producing 
positive  lift,  causes  a  positive  change  in  angle-of-attack  as 
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it  changes  the  direction  of  the  relative  wind.  Because  of 
this,  vertical  gusting  can  be  directly  modeled  as  a  change  in 
attack  angle.  If  a  step  upward  gust  could  occur,  an  instanta¬ 
neous  positive  change  in  attack  angle  would  result  before 
there  could  be  any  change  in  aircraft  motion.  Note,  however, 
that  this  is  apparently  in  conflict  with  the  previous  defini¬ 
tion  of  the  perturbation  attack  angle  given  in  Eqn  9.  In  this 
equation,  "w"  is  the  downward  velocity  of  the  aircraft  - 
which  implies  that  if  the  aircraft  doesn't  move,  there  can  be 
no  change  in  angle-of- at  tack. 

It  is  important  to  realize  that  the  stability-axis 
equations  model  aircraft  motion  relative  to  an  inertial 
reference  (approximated  as  the  earth).  The  motion  of  the 
surrounding  air  is  not  considered.  If  the  analysis  or  design 
to  be  accomplished  does  not  include  the  use  of  attack  angle 
feedback,  the  manner  by  which  angle-of-at  tack  is  modeled  is 
not  important  providing  the  resulting  response  is  correctly 
interpretted.  The  YF-16,  however,  uses  such  feedback,  and 
therefore  the  correct  quantity  must  be  employed.  This 
quantity,  which,  in  this  report,  is  called  true  perturbation 
attack  angle,  or  ,  is  defined  as  the  perturbation  from  trim 
in  the  angle  between  the  relative  wind  and  the  chord  of  the 
wing.  Note  that  this  is  not  the  same  as  'a,  which  is  defined 
as  the  perturbation  in  the  angle  between  the  velocity  vector 
of  the  aircraft  (in  inertial  space)  and  the  chord  of  the 
wing.  The  subtle  difference  in  these  two  definitions  is  that 
the  velocity  vector  of  the  aircraft  in  inertial  space  is  only 
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parallel  to  the  relative  wind  when  there  is  no  air  movement 
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The  velocity  of  moving  air  is  perhaps  most  easily  viewed 
from  the  earth  reference  frame.  A  vertical  gust  can  then  be 
represented  as  a  single  velocity  component  parallel  to  the  Ze 
axis  as  defined  in  Figure  2.  If  both  the  moving  air  and  the 
motion  of  the  aircraft  are  defined  in  this  same  axis  system, 
it  becomes  apparent  that  the  vertical  velocity  of  the  air¬ 
craft  relative  to  the  surrounding  air  is  simply  the  differ¬ 
ence  of  the  two  velocities  relative  to  the  earth.  If  w&/A  is 
defined  as  the  vertical  velocity  of  the  aircraft  relative  to 
the  surrounding  air,  is  defined  as  the  vertical  velocity 
of  the  aircraft  relative  to  the  earth,  and  wA/E  is  defined  as 
the  vertical  velocity  of  the  air  relative  to  the  earth,  then 

"a/*  =  "a/E  ‘  WA/E  (14) 


Dividing  by  the  equilibrium  forward  velocity  yields 


Wa/A  _  Wa/E  WA/E 
U0  =  ~o 


Since  in  a  perturbation  analysis,  U0  is  assumed  to  be  con¬ 
stant,  Eqn  15  implies  that 


%,/A  ~  %/E  ~  aA/E  (16) 

where  the  subscripts  are  as  defined  in  Eqn  14.  Since  is 
the  attack  angle  of  the  aircraft  relative  to  the  surrounding 
air,  and  'z  is  the  attack  angle  of  the  aircraft  relative  to 
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the  earth,  this  equation  may  be  expressed  as 


dm  =  ’a  +  cl. 


(17) 


where  ctg  is  a  shorthand  notation  for  -aA/E.  The  change  of  sign 
is  made  to  allow  an  upward  gust  to  be  defined  as  positive. 
The  error  associated  with  using  'a  instead  of  as  a  feed¬ 
back  source  is  demonstrated  in  Chapter  III,  where  an  analysis 
of  the  YF-16  CCV  longitudinal  flight  control  system  is 
performed.  , 


Modeling  a  Vertical  Gust 

As  atmospheric  disturbance  is  random  in  nature,  it  would 
seem  appropriate  that  it  should  be  modeled  as  such.  Several 
models  have  been  proposed,  most  of  them  statistically  based. 
The  disturbance  velocity  field  is  generally  assumed  to  be  a 
zero-mean  Gaussian  random  process  for  the  analysis  and 
simulation  of  linear  systems.  Although  there  is  sufficient 
evidence  to  prove  this  assumption  false,  most  analyses 
conducted  on  aircraft  systems  are  concerned  only  with  rela¬ 
tive  performance  (Ref  13).  For  this  reason,  MIL-F-8785B 
allows  the  use  of  either  of  two  statistically  based  random 
models  -  the  Dryden  or  von  Karmam  spectral  forms.  For  simula- 


tion 

purposes, 

many 

times  a 

filtered 

Gaussian  distributed 

noise 

source  is 

used. 

and  is 

considered  to  be 

an  acceptable 

model 

.  It  does 

not , 

however, 

produce 

enough 

extreme  gusts 

when  compared  to  actual  atmospheric  turbulence  data.  For  this 
reason,  MIL-F-8785B  requires  that  the  random  turbulence 
generated  in  this  way  be  supplemented  with  discrete  gusting. 
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For  the  purpose  of  comparing  aircraft  flight  control 
systems,  a  discrete  gust  model  may  be  more  appropriate  than  a 
random  gust  model.  In  addition  to  permitting  measurement  of 
peak  perturbations  in  aircraft  motion,  the  discrete  gust  is 
very  useful  in  determining  control  surface  actuator  rate  and 
displacement  saturation  effects.  For  these  reasons,  the 
discrete  gust  is  used  in  the  analyses  in  this  thesis. 

The  Discrete  Gust.  The  form  of  discrete  gust  selected  is 
as  it  appears  in  MIL-F-878 5B.'  This  form,  depicted  in  Figure 
4,  is  referred  to  as  a  "1-cos"  type  gust.  Its  name  comes  from 
the  fact  that  it  is  formed  by  displacing  an  inverted  cosine 
wave  a  unit  distance  above  the  distance  axis.  This  model 
provides  a  continuous  waveform  beginning  and  ending  at  zero 
velocity.  For  an  aircraft  flying  at  a  relatively  constant 
forward  velocity,  the  length  of  the  gust  may  be  thought  of  in 
terms  of  time  instead  of  distance.  For  the  simulations  in 
this  thesis,  the  duration  of  all  gusts  is  measured  in 
seconds . 

There  are  two  important  shaping  factors:  the  peak  value, 
or  maximum  gusting  velocity,  and  the  duration  of  the  gust.  It 
is  required  that  several  gust  durations  be  used  such  that  the 
gust  is  tuned  to  each  of  the  natural  frequencies  of  the 
aircraft  and  flight  control  system.  This  provides  for  a 
worst-case  analysis  in  addition  to  allowing  control  power 
effectiveness  to  be  evaluated  at  critical  frequencies.  The 
maximum  gusting  velocity  is  an  equally  critical  parameter  in 
such  an  analysis.  Although  a  linear  system  insures  a  one-to- 


cos 


Discrete  (l-cos)  vertical  gust  (Ref  13). 


one  relationship  between  input  and  output,  larger  inputs 


eventually  cause  rate  and/or  limit  saturation  of  the 
actuators. 

Throughout  this  report,  a  50  feet  per  second  ( fps)  peak 
velocity  is  used  for  the  discrete  gust.  At  the  selected 
flight  condition,  Mach  0.8  at  Sea  Level,  this  velocity 
corresponds  to  an  approximate  3.2  degree  peak  change  in 
relative  wind  direction.  The  standard  duration  used  is  one 
second  (one-half  second  to  peak  velocity),  except  for  simula¬ 
tions  at  critical  frequencies. 

Definition  of  Longitudinal  Tracking  Error 

In  performing  this  investigation,  it  is  necessary  to 
examine  the  use  of  various  feedback  configurations  in  an 
effort  to  produce  a  flight  control  system  that  makes  the 
target  tracking  task  easier  for  the  pilot.  To  determine  the 
performance  of  a  given  FCS,  a  measurement  tool  needs  to  be 
devised.  The  gust  alleviation  performance  of  the  Maneuver 
Enhancement  (ME)  system  of  the  YF-16  CCV  was  evaluated  by 
General  Dynamics  by  measuring  the  normal  acceleration 
response  at  the  pilot  station  (Fef  10).  There  was  no  atten¬ 
tion  given,  however,  to  the  effect  on  pitch  response.  It  is 
clear  that  aircraft  angular  rotation,  or  pitch  angle  res¬ 
ponse,  is  important  if  a  gun  (non-seeking  weapon)  is  fired  at 
a  distant  target.  As  an  example,  for  a  target  at  6000  feet,  a 
change  in  pitch  attitude  of  a  single  degree  results  in  a  miss 
distance  of  approximately  105  feet  if  inertia  and  accelera- 
tion  factors  are  neglected. 
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An  accurate  tracking  error  measurement  would  be  composed 
of  all  factors  defining  the  motion  of  both  target  and  track¬ 
ing  vehicle,  as  well  as  the  characteristics  of  the  projec¬ 
tile.  The  target  tracking  measurement  used  in  this  report  is 
greatly  simplified  compared  to  such  a  comprehensive  defini¬ 
tion.  It  does,  however,  provide  a  good  indication  of  tracking 
performance  for  an  aircraft  undergoing  minor  acceleration 
firing  at  a  fixed  target.  This  first-order  longitudinal 
approximation,  shown  in  Figure  5,  is  made  up  of  two  compo¬ 
nents.  The  tangent  of  the  pitch  angle  is  multiplied  by  the 
target  distance  to  calculate  the  tracking  error  due  to 
aircraft  rotation.  Then,  the  change  in  aircraft  vertical 
position  is  added  to  this  quantity  to  provide  for  miss 
distance  due  to  aircraft  vertical  translation.  The  sum  of 
these  two  components  yields 

TE  =  h  +  d  tan  6  (18) 
and  defines  target  tracking  error  as  used  in  this  report. 
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Definition  of  longitudinal  target  tracking  error 


There  are  essentially  two  significant  reasons  for  the 
analysis  performed  in  this  chapter.  First,  in  order  to 
evaluate  the  design  solution  presented  in  Chapter  IV,  it  is 
appropriate  to  compare  its  performance  to  the  performance  of 
the  current  YF-16  CCV  flight  control  systems.  Second,  by 
running  simulations  on  the  existing  flight  control  systems, 
the  simulation  results  can  be  checked  against  those  document¬ 
ed  by  General  Dynamics  (Refs  9,  10)  to  validate  the  simula¬ 
tion  approach.  A  rather  lengthy  digital  computer  program  is 
used  to  perform  the  simulations  using  a  multiple  state 
transition  matrix  technique  (see  Appendix  K  for  a  complete 
description),  and  a  verification  of  its  capability  provides 
assurance  that  it  performs  a  sufficiently  accurate  simula¬ 
tion. 

In  the  pages  that  follow,  the  baseline  YF-16  flight 
control  system  is  described,  as  is  the  integration  of  the 
additional  circuitry  needed  to  provide  the  CCV  Maneuver 
Enhancement  function.  These  systems  are  then  simplified  and 
linearized  to  facilitate  simulation.  A  system  time  response 
is  then  performed,  employing  both  pilot  and  gust  disturbance 
inputs.  As  the  flight  control  system  uses  angle-of-attack 
feedback,  simulations  are  run  using  both  true  perturbation 
attack  angle  and  stability-axis  perturbation  attack  angle,  so 
that  the  results  may  be  compared  (see  Chapter  II).  An  empiri- 


cal  frequency  response  analysis  is  also  performed  for  the 
system  in  response  to  discrete  vertical  gusts.  The  chapter 
concludes  with  a  brief  discussion  of  the  results  obtained. 

Control  System  Functional  Description 

The  longitudinal  flight  control  system  of  the  YF-16  is 
illustrated  in  the  functional  block  diagram  of  Figure  6.  The 
solid  line  connecting  paths  represent  the  baseline  YF-16  FCS, 

and  the  dashed  line  connecting  paths  represent  the  implemen- 

/ 

tation  of  the  Maneuver  Enhancement  augmentation.  In  the 
following  paragraphs,  the  function  of  each  block  is  briefly 
described  in  an  effort  to  provide  the  reader  with  a  basic 
understanding  of  the  design  philosophy. 

Baseline  YF-16  FCS.  In  the  design  of  the  longitudinal 
axis  portion  of  YF-16  flight  control  system.  General  Dynamics 
chose  to  use  three  measurable  quantities  as  feedback  signals: 
angle-of-attack  to  allow  for  variance  of  static  stability  by 
scheduled  gain,  pitch  rate  for  damping,  and  normal  accelera¬ 
tion  to  provide  g-command  capability.  Electrical  signals  are 
transmitted  from  the  redundant  sensors,  gyros,  and  accelero¬ 
meters  that  measure  these  quantities  to  the  analog  Basic 
Flight  Control  Computer.  The  pilot  inputs  longitudinal 
•  commands  to  this  computer  via  a  force-sensing  sidestick 

controller.  As  shown  in  Figure  6,  this  input  is  reshaped  by 
nonlinear  stick  force  gradients.  These  nonlinearities  tailor 
the  force-stick  output  by  allowing  for  mechanical  preloading 
and  electrical  breakout,  as  well  as  a  pitch  command  gradient. 


28 


Figure  6.  YF-16  CCV  flight  control  system  (baseline  FCS  with  Maneuver  Enhancement). 


In  addition  to  the  force-stick,  the  pilot  is  also  permitted  a 
tr  iir  adjustment. 


The  command  prefilter  is  simply  a  unity-gain  first-order 
lag  that  is  used  tc  filter  pilot  command  transients  to  a 
desired  rate.  The  variable  command  limit,  a  nonlinear  gain 
computed  as  a  function  of  ratios  of  static  and  impact  pres¬ 
sures,  safeguards  the  aircraft  from  control  inputs  that  would 
result  in  exceeding  the  permissible  wing-root  bending  moment 
at  critical  flight  conditions.  Depending  on  the  flight 
condition,  this  command  limit  varies  frcm  a  minimum  cf 


approximately  5.5  Gs  to  a  maximum  of  8  Gs.  The  output  of  this 
limiter  may  be  considered  as  the  effective  system,  input. 

Variable  forward-loop  gain  is  provided  as  a  function  cf 
impact  pressure  which  schedules  the  gain  according  to 
airspeed.  A  proportional  fcrward-locp  path  provides  for  quick 
response  tc  commanded  inputs,  while  an  integrated  path  forces 
a  zero  steady-state  error  condition  to  exist  for  those  seme 
inputs.  As  the  feedback  summed  in  prior  tc  that  integration 
consists  of  normal  acceleration  and  washed-out  pitch  rate. 


the  steady-state  response  is  equal  to  the  commanded  Gs.  This 
type  of  flight  control  system  provides  the  pilot  with  what  is 
essentially  a  C*  command  aircraft  (see  Appendix  J  for  a 
discussion  of  the  C*  concept).  During  command  transients,  the 
aircraft  responds  with  pitch  rate  changes  not  unlike  the 
response  from  a  pitch  rate  commanded  FCS.  In  steady-state 
however,  the  pilot  is  commanding  Gs,  as  the  output  of  the 
washout  in  the  pitch  rate  feedback  loop  goes  to  zero. 
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addition,  this  design  provides  a  series  trim  to  maintain  a 
1-G  trimmed  condition  in  the  absence  of  a  pilot  input.  The 
lag  and  lead  compensation  in  the  pitch  rate  and  normal 
acceleration  feedback  paths  provide  the  phase  relationships 
which  result  in  desired  system  performance  characteristics. 

Attack  angle  feedback,  used  to  control  static  stability, 
is  summed  in  after  the  forward-loop  integration.  The  gain  in 
this  feedback  path  is  scheduled  according  to  Mach  number.  It 
is  this  gain  that  allows  the  Longitudinal  static  stability  to 
be  increased  for  sub-sonic  flight,  and  decreased  for  super¬ 
sonic  flight  as  discussed  in  Chapter  I.  Attack  angle  feedback 
is  also  used  to  provide  angle-of-attack  limiting  for  automat¬ 
ic  spin  prevention.  As  the  angle-of-attack  approaches  a 
critical  point  in  flight,  a  nonlinear  multiple  of  attack 
angle  feedback  is  used  to  reduce  the  pilot's  commanded  input 
to  prevent  further  increases  in  angle-of-attack.  A  variable 
amount  of  pitch  rate  input  to  this  limiter  is  used  to  provide 
the  additional  system  damping  that  is  desirable  at  these  high 
angle-of-attack  conditions. 

The  output  of  the  flight  control  computer  is  routed  to 
the  electromechanical  horizontal  tail  command  servos.  The 
mechanical  outputs  of  these  servos  drives  the  hydraulic 
horizontal  tail  power  actuators.  The  horizontal  tail  is 
designed  to  permit  independent  left-hand  and  right-hand 
deflections,  which  allow  for  smoother  roll  control.  For 
longitudinal  control,  however,  the  two  command  servos  are  fed 
symmetr ically. 
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Maneuver  Enhancement  Augmentation,  The  Maneuver  Enhance¬ 
ment  system  is  based  on  commanding  a  symmetrical  deflection 
of  the  trailing-edge  flaperons  as  a  consequence  of  a  non-zero 
differential  between  commanded  Gs  and  normal  acceleration.  As 
a  portion  of  the  CCV  flight  control  package,  this  system  is 
implemented  in  an  analog  Auxiliary  Flight  Control  Computer. 

Commanded  input  is  obtained  from  the  Basic  Flight  Control 

Computer  using  the  signal  at  the  output  of  the  command 

/ 

prefilter.  This  quantity  is  compared  to  the  measured  normal 
acceleration  to  provide  an  error  signal.  The  algebraic  signs 
of  the  comparator  inputs  are  opposite  to  those  in  the  basic 
FCS,  which  produces  an  error  signal  of  opposite  polarity  to 
that  driving  the  horizontal  tail.  Recall  that  positive 
rotation  of  the  aircraft  is  obtained  by  an  upward  (negative) 
horizontal  tail  deflection.  Positive  translation,  on  the 
other  hand,  is  acquired  as  a  result  of  downward  (positive) 
flap  deflection. 

The  flap  G-error  signal,  multiplied  by  a  scheduled  gain, 
is  used  in  three  ways.  It  is  used  directly  to  command  trail¬ 
ing-edge  flaperon  deflections  by  symmetrically  driving  the 
flaperon  command  servos.  These  servos  operate  the  hydraulic 
power  actuators  which  in  turn  deflect  the  flaperon  surfaces. 
It  is  also  used  to  provide  a  horizontal  tail/trailing-edge 
flaperon  interconnect  by  symmetrically  driving  the  horizontal 
tail  command  servos  through  a  scheduled  gain.  This  intercon¬ 
nect  is  intended  to  allow  sufficient  tail  deflection  to 
offset  the  pitching  moment  created  by  flaperon  deflection. 
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This  allows  a  direct  lift  effect  with  no  change  in  attack 
angle  during  the  early  transient  period  following  a  commanded 
stick  input.  Finally,  the  error  signal  is  used  to  provide  a 
pitch  rate  bias.  This  bias  signal,  summed  in  with  the  pitch 
rate  feedback  to  the  basic  FCS,  is  intended  to  negate  this 
feedback  and  prevent  an  opposing  horizontal  tail  deflection 
during  the  early  transient  period  (Ref  9). 

Before  leaving  this  subject,  it  is  important  to  note  that 
since  the  horizontal  tail/tra }ling-edge  flaperon  interconnect 
and  pitch  rate  bias  signals  are  generated  in  an  open-loop 

manner,  the  intended  function  of  these  two  signals  is  clearly 
dependent  upon  the  accuracy  of  the  scheduled  gains.  Also,  as 
previously  mentioned,  operation  at  critical  attack  angles 
causes  a  non-zero  output  from  the  angle-of-attack  limiter  in 
the  basic  FCS,  which  results  in  reduced  pilot  command 

authority.  As  this  occurs,  the  flaperon  G-error  signal 
increases,  commanding  large  flaperon  deflections.  These  large 
deflections  result  in  increased  drag  which  degrades  aircraft 
performance. 

Linearization  for  Analysis 

Baseline  YF-16  FCS.  To  permit  a  linear  simulation,  the 
YF-16  flight  control  system  nonlinearities  must  first  be 
simplified.  Most  of  these  nonlinearities  have  a  piecewise- 
linear  characteristic,  that  allows  reduction  to  a  simple  gain 
if  signal  level  restrictions  are  adhered  to.  The  simplest 

case  is  the  hard  or  soft  limiter  with  a  dead  zone.  If  the 

input  level  is  assumed  to  remain  within  the  dead  zone,  the 
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output  is  zero,  and  the  function  may  be  removed  completely. 
Scheduled  gains,  on  the  other  hand,  may  be  easily  determined 
by  chosing  a  flight  condition.  These  gains  are  defined  in 
Appendix  D.  Finally,  as  the  trimmed  angle-of-attack  at  the 
selected  flight  condition  is  relatively  small  (see  Appendix 
B),  it  is  assumed  that  the  angle-of-attack  limiter  output  is 
negligible.  General  Dynamics  apparently  makes  a  similar 
assumption  in  their  linear  analysis  performed  in  Reference 
10.  , 

The  flight  condition  selected  for  the  analyses  in  this 
chapter,  as  well  as  in  Chapter  IV,  is  0.8  Mach  at  Sea  Level. 
This  FC  is  chosen  for  two  reasons.  First,  it  is  near  the 
transonic  region,  and  is  a  condition  quite  sensitive  to 
terrain- induced  and  atmospheric  turbulence.  Second,  it  is  the 
flight  condition  chosen  for  the  disturbance  analysis  accom¬ 
plished  by  General  Dynamics  (Ref  10).  Use  of  this  identical 
flight  condition  enables  a  comparison  of  results,  and  permits 
verification  of  the  simulation  technique. 

The  linearized  model  of  the  basic  YF-16  flight  control 
system  is  shown  in  Figure  7.  The  forward  loop  consists  of  the 
command  prefilter,  the  scheduled  gain  F3  and  P+1  network,  and 
the  command  servo  and  power  actuator.  There  are,  as  shown  in 
Figure  6,  actually  two  command  servos  and  two  power  actua¬ 
tors.  However,  as  they  are  driven  and  operate  symmetrically 
and  both  sides  of  the  tail  must  move  to  constitute  a  horizon¬ 
tal  tail  deflection,  modeling  them  as  a  single  servo  and 
actuator  is  valid.  Note  also  in  comparing  this  figure  to 
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Linearized  YF-16  baseline  flight  control  system. 


Figure  6  that  algebraic  signs  at  comparator  inputs  are 
changed,  and  the  command  servo  gain  is  shown  as  negative. 
These  changes  are  made  to  present  the  system  in  a  more 
familiar  form. 

Although  not  shown  in  the  figure,  it  is  important  to  note 
that  the  power  actuator  rate  saturates  at  60  degrees  per 
second  with  a  maximum  displacement  of  25  degrees  in  either 
direction.  Also  of  interest  is  the  fact  that  an  electrical 
resistance  proportional  to  /displacement  is  available  from 
both  the  command  servo  and  the  power  actuator.  Appendix  1 
details  the  design  of  a  first-order  observer  to  permit  access 
to  the  internal  state  of  the  command  servo. 

The  feedback  paths  consist  of  the  appropriate  sensor  and 
simple  first-order  networks,  with  only  the  attack  angle 
feedback  path  containing  a  scheduled  gain.  The  model  for  the 
angle-of- attack  sensor  is  the  first-order  lag 


TF,  0  x  = 
(  a  Sens  ) 


66.7 
s  +  66.7 


(19) 


The  models  for  the  pitch  rate  gyro  and  normal  accelerometer 
are  identical,  and  are  given  as  the  second-order  lag 


TF(  q  Gyro)  TF  (  AN  Acc  ) 


98596 _ 

s2  +  376.8  s  +  98596 


(20) 


The  system  of  Figure  7  is  further  simplified  by  combining 
blocks,  returning  all  feedback  paths  to  the  same  point,  and 
moving  the  command  prefilter  inside  the  feedback  closure.  The 


result  is  shown  in  Figure  8.  The  only  minor  complication  is 
the  occurance  of  an  unrealizable  transfer  function  in  the 
pitch  rate  feedback  path.  For  simulation  purposes,  the  free 
differentiator  was  removed  and  pitch  acceleration  feedback 
used  instead  of  pitch  rate.  This  option  in  the  simulation 
program  allows  a  more  general  case  FCS  to  be  simulated  (see 
Appendix  K  for  further  discussion). 

Maneuver  Enhancement  Augmentation.  The  linearized  model 
of  the  YF-16  FCS  with  the  Maneuver  Enhancement  system  engaged 
is  somewhat  more  complicated  as  indicated  in  Figure  9.  The 
upper  portion  of  this  figure  is  identical  to  Figure  7,  with 
the  lower  portion  showing  the  Maneuver  Enhancement  augmenta¬ 
tion.  Although  the  models  of  the  flaperon  channel  command 
servo  and  power  actuator  are  dynamically  identical  to  those 
in  the  horizontal  tail  channel,  there  are  a  few  minor 
differences  that  should  be  addressed.  As  indicated  in  the 
figure,  the  sign  of  the  servo  gain  in  the  flaperon  channel  is 
positive,  and  is  modeled  as  such  to  reflect  the  described 
sign  convention  (Chapter  II).  Also,  the  power  actuator  rate 
saturates  at  56  degrees  per  second,  and  is  displacement 
limited  to  20  degrees  in  either  direction. 

The  system  of  Figure  9  may  be  greatly  simplified  by 
combining  blocks  and  moving  feedback  and  feedforward  paths  to 
common  nodes.  Leaving  out  the  details,  this  reduction  is 
given  in  Figure  10.  The  system  may  be  further  simplified  by 
moving  the  resultant  fourth-order  "prefilter"  inside  the 
horizontal  tail  feeotack  locp  closure.  The  factored  form  of 
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Linearized  YF-16  CCV  flight  control  system. 


10.  Simplified,  linearized  YF-16  GCV  flight  control  system. 


this  prefliter  is 


56.59  (s  +  .4504)  (s  ♦  2.062)  (s  +  11.85) 

(s  +  l)(s  +  5) (s  +  8.3) (s  +  15) 

and  the  factored  form  of  the  remaining  horizontal  tail 
forward  loop  block  is 


TF 


-21929  (s  +  5) _ 

(s  +  20)  (s  +  36.4  ±  J37.14)s 


(22) 


In  combining  these  two  blocks,  a  single  cancellation  is  made, 
yie lding 


Tp  _  -1240962  (s  +  .4504)  (a  +  2.062)  (B  +  11.85) _ 

(s  +  l)(s  +  8.3)  (s  +  15)  (s  +  20)  (s  +  36.4  ±j  37.l4)s  (23) 

Due  to  the  tact  that  the  simulation  computer  program  cannot 
currently  handle  transfer  functions  in  excess  of  fifth  order, 
some  simplifying  approximation  must  be  made.  One  straight¬ 
forward  approach  is  to  assume  that  the  second  order  dynamics 
of  the  servo  are  non-dominant.  As  there  is  a  pole  at  -1  on 
the  real  axis  with  no  zero  nullifying  its  dominance,  this 
assumption  is  valid.  There  is  sufficient  cancellation  of 
poles  and  zeros  for  the  feedback  transfer  functions,  there¬ 
fore  no  similar  approximations  are  necessary.  The  final 
system  model  is  given  in  Figure  11. 

System  Time  Responses 

With  the  aircraft  and  flight  control  system  suitably 
modeled,  a  simulation  of  the  composite  system  can  now  be 
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Simulation  configuration  of  the  linearized  YF-16  CCV  PCS 


effected  for  the  selected  flight  condition  (0.8  Mach  §  SL) . 
The  simulation  includes  system  state  responses  to  pilot  and 
disturbance  inputs  for  the  YF-16  FCS  with  the  Maneuver 
Enhancement  augmentation  both  out  of  circuit  (Figure  8)  and 
in  circuit  (Figure  11). 

Computer  generated  plots  of  selected  quantities  are 
presented  as  simulation  output.  Angle-of-a ttack  is  shown  in 
two  ways,  since  two  forms  of  this  quantity  have  been 
identified.  Those  plots  labeled  "ANGLE  OF  ATTACK"  show  the 
stability-axis  perturbation  attack  angle  response.  Plots 
labeled  "OUTPUT  OF  ATTACK  ANGLE  SENSOR"  reflect  the  form  of 
attack  angle  that  is  being  used  as  feedback,  modified  by  the 
specified  dynamics  of  the  attack  angle  sensor. 

Horizontal  tail  deflection  and  trailing-edge  flaperon 
deflection  are  also  among  the  quantities  plotted.  In  addition 
to  the  plotted  response,  maximum  achieved  rate  of  deflection 
and  rate  limiter  disposition  are  printed  on  these  plots.  If  a 
rate  saturation  limit  for  the  actuator  is  specified,  the  flag 
"LIMIT, ON"  is  printed.  If  no  limit  is  defined,  the  flag 
•LIMIT, OFF"  is  printed.  Recall  that  rate  saturation  of  the 
horizontal  tail  actuators  occurs  at  60  degrees  per  second, 
and  for  the  trailing-edge  flaperon  actuators,  at  56  degrees 
per  second.  For  those  simulations  wnere  a  rate  limit  is 
specified,  these  values  are  utilized. 

Response  to  Pilot  Step  Input.  To  simulate  system  response 
to  a  pilot  command,  a  1-G  step  input  is  applied.  System  re¬ 
sponse  with  Maneuver  Enhancement  switched  off  is  shown  in 


Figure  12,  while  response  with  Maneuver  Enhancement  switched 
on  is  shown  in  Figure  13.  The  input  dimensionality  is  given 
in  volts.  This  is  done  for  convenience  only,  since,  depending 
on  the  system  being  simulated,  the  commanded  response  can  be 
any  state  or  a  linear  combination  of  states.  As  the  actual 
input  to  the  FCS  is  a  voltage,  the  label  "volts"  is  chosen  as 
a  representative  dimensionality. 

Examination  of  the  system  response  characteristics 
portrayed  in  these  figures  reveals  the  intended  function  of 
Maneuver  Enhancement  -  to  improve  transient  performance.  Note 
that  the  onset  of  both  pitch  rate  and  normal  acceleration  are 
substantially  quickened  relative  to  the  basic  YF-16  FCS 
responses.  This  initial  transient  improvement  is  due  to  the 
deflection  of  the  tra iling-edge  flaperons.  With  horizontal 
tail  control  only,  changes  in  normal  acceleration  are  accom¬ 
plished  by  rotation  of  the  aircraft.  The  rate  of  response  is, 
therefore,  dependent  on  the  rate  of  aircraft  rotation.  The 
tr ailing-edge  flaperons,  however,  directly  increase  lift  on 
the  wing,  and  permit  a  more  rapid  increase  in  normal 
acceleration. 

Approximately  one  second  after  application  of  the  step 
input,  the  normal  acceleration  response  of  both  systems  is 
essentially  the  same  -  a  slow  increase  over  the  succeeding 
few  seconds  to  the  commanded  level.  For  the  Maneuver  Enhance¬ 
ment  system,  the  transition  from  flaperon  induced  translation 
to  horizontal  tail  induced  rotation  is  clearly  defined  at  the 
0*3  second  point.  This  interim  pause  in  G-onset  is  a  somewhat 


undesirable  handling  characteristic  caused  by  the  extremely 
sudden  change  in  flaperon  deflection  rate. 

Although  the  improvement  in  initial  transient  response  is 
clearly  visible,  it  is  important  to  note  that  the  responses 
shown  are  for  a  1-G  step  command.  Even  at  this  low  input 
level,  trailing-edge  flaperon  deflection  rate  reaches  approx¬ 
imately  one-half  the  actuator  saturation  rate.  As  the  air¬ 
craft  is  capable  of  pulling  in  excess  of  7  Gs  over  straight 
and  level  flight,  an  actuatpr  rate  of  close  to  200  degrees 
per  second  would  be  required  to  obtain  the  same  transient 
response  characteristics  shown  in  Figure  13.  As  the  flaperons 
are  clearly  responsible  for  the  improvement  in  transient 
response  obtained  with  Maneuver  Enhancement,  a  flaperon 
deflection  rate  limit  of  56  degrees  per  second  decidedly 
impacts  the  amount  of  improvement  available  during  high- 
performance  maneuvers. 

Response  to  the  Discrete  Vertical  Gust.  To  simulate  sys¬ 
tem  response  to  atmospheric  disturbances,  discrete  (1-cos) 
vertical  gust  inputs  are  applied.  Peak  gust  velocity  is  50 
feet  per  second,  with  a  duration  of  one  second.  This  gust  was 
chosen  to  be  consistant  with  that  used  in  the  disturbance 
analysis  performed  by  General  Dynamics. 

To  enable  a  direct  comparison  of  the  system  responses 
that  are  obtained  as  a  result  of  feeding  back  each  of  the  two 
previously  defined  forms  of  angle-of-attack,  simulations  are 
run  for  each  case.  In  addition,  responses  are  shown  for  the 
FCS  with  and  without  Maneuver  Enhancement. 
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System  response  with  ME  off  and  'a  feedback  is  shown  in 
Figure  14.  In  this  case  ,  the  measured  angle-of-attack  is 
assumed  to.be  'a  modified  by  the  Sensor  dynamics.  Therefore, 
the  two  traces  are  almost  identical.  System  response  with  ME 
off  and  %  feedback  is  shown  in  Figure  15.  The  difference 
between  the  measured  attack  angle  and  'a  is  apparent. 

Comparing  the  two  figures,  the  peak  perturbation  in 
aircraft  pitch  angle  is  actually  almost  80%  greater  when 
using  feedback  than  that  pbtained  by  incorrectly  using  'a 
as  the  feedback  quantity.  This  results  in  a  corresponding 
increase  in  target  tracking  error.  There  is  also  a  change  in 
the  number  of  negative  Gs  pulled  as  a  result  of  the  gust,  as 
well  as  an  increase  in  horizontal  tail  deflection. 

Figures  16  and  17  shew  the  system  response  with  ME  on, 
using  'a  and  dm  feedback,  respectively.  As  in  the  previous 
case  for  ME  off,  use  of  Qg  as  a  feedback  quantity  provides  a 
more  accurate  simulation  of  actual  system  performance. 

The  gust  alleviation  properties  of  Maneuver  Enhancement 
can  be  recognized  from  a  comparison  of  Figures  15  and  17.  For 
a  discrete  gust  of  one  second  duration,  pitch  angle  pertur¬ 
bation  and  peak  target  tracking  error  are  reduced  by  about 
10%,  and  normal  acceleration  perturbation  is  improved  by  more 
than  40%.  The  large  improvement  in  normal  acceleration 
response  is  the  natural  result  of  feeding  back  only  normal 
acceleration  to  the  trailing-edge  flaperons.  Any  deviation 
from  the  commanded  value  forces  a  flaperon  deflection  that 
attempts  to  nuil  only  that  error. 


Figure  13b.  YF-16  CCV  FCS  (with  Maneuver  Enhancement)  response  to  a  pilot  step  input 


Figure  lUa.  YF-16  baseline  FCS  with  'a  feedback  response  to  a  discrete 


14b.  YF-16  baseline  FCS  with  'a  feedback  response  to  a  discrete  gust 


Figure  15^«  YF-16  baseline  FCS  with  otm  feedback  response  to  a  discrete  gust. 


GUST  DISTURBANCE  INPUT  .  OUTPUT  OF  ATTACK  ANGLE  SENSOR 


FCS  with  'a  feedback  response  to  a  discrete  gust 


l6to.  YF-16  CCV  FCS  with  'a  feedback  response  to  a  discrete 


OUST  DISTURBANCE  INPUT  .  OUTPUT  OF  ATTACK  ANOLE  SENSOR 


Figure  17a.  YF-16  CCV  FCS  with  feedback  response  to  a  discrete  gust. 


Frequency  Response  to  Discrete  Vertical  Gusts 

To  permit  a  more  encompassing  assessment  of  system 
disturbance  response  characteristics,  gusts  of  various 
frequencies  should  be  applied  to  the  system.  A  single  gust  of 
one  second  duration  dees  not  give  a  true  indication  of  the 
performance  of  the  system  in  the  presence  of  gusts  that  are 
of  shorter  or  longer  duration.  Digesting  page  after  page  of 

simulation  output  to  get  that  overall  picture  is,  however,  a 
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formidable  task.  A  type  of  frequency  response  analysis  is  the 
obvious  answer. 

To  determine  the  frequency  response  analytically,  the 
closed-loop  system  transfer  functions  must  be  known.  Although 
the  transfer  function  for  the  system  without  Maneuver 
Enhancement  can  be  calculated  with  little  difficulty,  the 
same  cannot  be  said  about  the  system  with  Maneuver  Enhance¬ 
ment.  The  added  complication  of  running  feedback  loops  to  an 
additional  input  -  when  the  measurable  outputs  are  dependent 
on  both  inputs  -  increases  the  complexity  by  several  orders 
of  magnitude.  Even  if  the  transfer  functions  could  be  deter¬ 
mined,  it  would  perhaps  be  more  informative  to  examine  system 
response  to  discrete  gusts  of  varying  durations,  as  opposed 
to  applying  a  sine  wave  input  of  varying  frequency. 

One  way  of  obtaining  such  data  is  to  first  run  several 
simulations  applying  discrete  gusts  of  varying  durations  to 
the  system.  Then,  for  each  run,  the  absolute  value  of  the 
maximum  perturbation  of  a  given  quantity  can  be  plotted 
against  the  frequency  of  the  gust  applied.  The  gust  frequency 
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is  simply  the  multiplicative  inverse  of  the  gust  duration. 
There  is,  of  course,  no  periodic  frequency  associated  with  a 
discrete  gust,  but  this  approach  permits  the  data  to  be 
plotted  in  a  more  familiar  form. 

For  this  analysis,  ten  simulations  are  run  for  each  of 
the  four  system  configurations  -  'a  or  otg,  feedback  with  and 
without  Maneuver  Enhancement.  Gust  durations  covering  a  range 
of  .03  to  30  seconds  are  applied,  and  in  each  case,  maximum 
gust  velocity  is  50  feet  per  second.  The  results  are  shown 
for  seven  quantities  of  interest  in  Figures  18  through  24. 
The  effect  of  using  'a  feedback  is  clearly  visible,  as  is  the 
influence  of  Maneuver  Enhancement.  A  small  marker  placed  on 
the  frequency  axis  in  each  figure  indicates  the  gust  duration 
that  was  used  in  the  previous  time  response  analysis. 


Discussion 

It  has  been  demonstrated  that  aircraft  performance  can 
indeed  be  improved  through  the  use  of  the  trailing-edge 
flaperons  as  an  active  control  surface.  In  response  to  pilot 
commands,  initial  transient  response  is  significantly  enhanc¬ 
ed  by  the  flaperon  deflections  which  supplement  the  normal 
aircraft  rotation.  In  turbulence,  the  flaperon  movements 
result  in  substantially  reducing  undesirable  perturbations  in 
normal  acceleration. 

Also  demonstrated  was  the  importance  of  using  the  correct 
measure  of  angle-of-at tack  as  a  feedback  quantity  to  model 
the  actual  system.  In  the  case  presented,  'a  is  not  the 
attack  angle  physically  measured  by  the  sensing  device.  As 
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demonstrated,  'a  is  only  correctly  used  when  the  analysis 
being  performed  excludes  the  effect  of  external  disturbances 
on  the  aircraft  (i.e.  ,  when  determining  aircraft  response  to 
pilot  inputs).  It  is  shown  in  Equation  17  that  when  there  is 
no  movement  of  surrounding  air  (  Og  =  0),  then,  and  only  then, 
is  'a  equal  to  0^. 

Although  the  Maneuver  Enhancement  augmentation  designed 
by  General  Dynamics  does  function  as  a  gust  alleviation 
system  which  provides  improvement  in  normal  acceleration 
response,  it  provides  little  or  no  improvement  in  pitch 
response.  As  a  result,  it  provides  no  appreciable  reduction 
in  target  tracking  error  as  defined  herein.  The  major  reason 
for  this  is  that  the  system  was  designed  primarily  to  enhance 
aircraft  response  to  pilot  commands  -  thus.  Maneuver  Enhance¬ 
ment.  Its  gust  alleviation  performance  was  essentially  a 
fallout  of  the  chosen  configuration  employing  normal  acceler¬ 
ation  feedback  to  the  flaperons. 

It  is  likely  that  aircraft  response  to  turbulence  can  be 
further  reduced  if  that  is  the  design  objective.  One 
possibility  of  interest  might  be  the  use  of  attack  angle 
and/or  pitch  rate  feedback  to  the  flaperons.  This  is  the 
subject  of  the  next  chapter. 
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FREQUENCY  (Radians/  Second 


In  this  chapter,  a  flight  control  system  design  is 
presented.  The  design  objective  is  first  addressed,  followed 
by  a  discussion  of  the  design  approach  alternatives  that  were 
attempted.  The  range  of  applicability  of  the  rigid-body 

assumption  for  the  longitudinal  axis  of  the  YF-16  is  then 
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briefly  discussed.  This  is  done  to  provide  some  insight  as  to 
the  limitations  which  must  be  imposed  on  selection  of  the 
dominant  closed-loop  eigenvalues. 

The  chosen  design  methodology  utilizes  a  matrix  represen¬ 
tation  of  the  aircraft  equations  of  motion,  and  therefore  a 
set  of  state-space  equations  -  including  all  the  required 
measurable  outputs  -  is  determined.  Of  the  several  possible 
feedback  configurations  that  were  examined,  a  pitch  rate 
command  system  is  selected  and  modeled  in  matrix  format. 
Detailed  modeling  of  the  rejected  systems  is  contained  in 
Appendix  F. 

The  transfer  functions  of  interest  are  determined  by 
employing  Cramer's  Pule,  and  then  briefly  analyzed.  An 
iterative  approach  to  finding  acceptible  solutions  to  these 
transfer  functions  is  used  to  derive  the  required  gains  in 
two  design  examples.  The  primary  design  is  simulated  with 
both  pilot  and  disturbance  inputs,  anc.  compared  to  the 
b tmu 1  at  ions  in  Chapter  III.  Both  analytical  and  empirical 
equency  response  data  is  also  presented.  The  secondary 
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design  is  simulated  for  discrete  gust  input  only,  to  compare 
its  performance  with  the  primary  design. 


Design  Objective 

The  first  step  in  the  design  of  any  control  system  must 
be  the  definition  of  system  performance  requirements  -  a 
design  objective.  In  this  thesis,  that  design  objective  is 
the  development  of  a  flight  control  system  that  improves  the 

tracking  performance  and  disturbance  rejection  properties  of 

/ 

the  aircraft,  without  degrading  overall  flying  or  handling 
qualities. 

There  are  a  few  undesirable  characteristics  of  the 
previously  examined  Maneuver  Enhancement  system  that  should 
be  avoided  in  the  design,  if  possible.  As  mentioned  in 
Chapter  III,  operation  of  the  aircraft  at  critical  attack 
angles  causes  a  non-zero  output  from  the  angle-of-attack 
limiter,  resulting  in  increasing  reductions  in  pilot  command 
authority.  As  the  ME  system  is  not  affected  by  this  limiter, 
increased  pilot  command  levels  result  in  large  flaperon 
deflections.  These  large  deflections  result  in  increased 
overall  drag  at  a  flight  condition  where  maximum  aircraft 
performance  is  essential.  It  is  therefore  desirable  to 
eliminate  this  problem  in  the  design. 

Also  mentioned  in  Chapter  III  is  the  fact  that  the 
tr ailing-edge  flaperon  power  actuators  are  too  quickly  rate 
saturated.  A  pilot  command  in  excess  of  two  Gs  results  in 
actuator  rate  saturation  and  subsequent  degradation  of 
aircraft  performance  with  increased  pilot  command  levels.  It 


would  be  most  desirable  to  maintain  flaperon  effectiveness 
over  the  entire  range  of  pilot  command  authority. 

A  third  difficulty  with  the  ME  system  is  the  interim 
pause  in  G-onset  caused  by  a  sudden  change  in  flaperon 
deflection  rate.  As  previously  discussed,  this  is  the  point 
of  transition  from  flaperon  induced  translation  to  horizontal 
tail  induced  rotation  of  the  aircraft.  A  softening  of  this 
transition  point  is  a  desirable  design  feature. 


Design  Approach  Alternatives 


The  intent  of  this  section  is  to  provide  some  background 
and  rationale  for  the  design  approach  presented  in  this 
thesis.  Many  hours  were  invested  in  attempts  to  use  other 
methods  -  each  attempt  resulting  in  failure  for  one  reason  or 
another.  The  selected  approach,  straightforward  though 
iterative  in  nature,  was  the  only  approach  that  produced  an 
acceptable  design.  In  the  paragraphs  that  follow,  each  of  the 
methods  attempted  is  briefly  discussed. 


Root  Locus .  A  classical  design  approach,  this  method  was 
tried  as  a  first  attempt  toward  determining  an  acceptible 
design.  It  provides  visibility  into  the  locations  of  poles 
and  zeros  and  permits  easy  incorporation  of  dynamic  feedback 
compensation.  However,  although  it  can  be  applied  to  multi¬ 
ple-output  systems  by  appropriate  modeling  and  successive 
feedback  loop  closure,  it  does  not  lend  itself  well  to 
multiple-input  systems.  For  a  two-input  system,  each  half  of 
the  system  must  be  designed  separately,  and  the  designer  has 
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little  or  no  feel  for  the  performance  of  the  composite 
system. 

Modal  Control.  This  approach  can  be  easily  applied  to 
systems  with  any  number  of  inputs  or  outputs.  Also,  the 
desired  closed-loop  eigenvalues  can  be  selected  according  to 
a  model  system.  One  problem  encountered  with  this  method, 
however,  is  the  need  to  make  all  states  available  for  feed¬ 
back.  This  requires  the  design  of  a  first-order  observer  to 
"observe"  the  internal  state  of  the  second-order  servos  (see 
Appendix  I).  In  addition,  it  is  very  difficult  to  implement 
dynamic  compensation  -  if  desired  -  using  this  method.  The 
fatal  difficulty,  hcwever,  is  the  fact  that  there  is  no 
positive  control  over  zero  locations.  As  a  result,  although 
it  was  not  difficult  to  design  a  system  that  responded  well 
to  pilot  commands,  no  way  was  found  to  design  a  system  to 
reject  disturbances. 

Optimal  Control.  This,  method,  also  directly  applicable  to 
multiple-input  multiple-output  systems,  was  attempted  using  a 
Quadratic  Performance  Index.  As  with  modal  control,  access  to 
all  system  states  is  required  -  in  this  case,  for  a  solution 
to  the  algebraic  Ficatti  Equation.  This  method,  also,  is 
based  upon  feedback  of  fixed-gain  multiples  of  system  states, 
and  therefore  dynamic  compensation  to  increase  closed-loop 
system  flexibility  is  extremely  difficult  to  implement. 

A  frequent  obstacle  encountered  when  attacking  the  opti¬ 
mal  servo  problem  is  the  selection  of  the  weighting  matrices. 
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A  related  drawback  is  the  lack  of  positive  control  over 
placement  of  zeros,  as  well  as  the  limited  range  of  locations 
available.  These  problems,  when  combined  perhaps  with  an 
insufficient  knowledge  of  the  design  or  application  theory, 
prevented  the  design  of  a  system  that  met  the  specified 
disturbance  rejection  requirements. 

State  Variable  Feedback  (Single  Input).  This  method  not 

only  permits  the  assignment  of  eigenvalues  by  designing  to  a 
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selected  model,  but  also  allows  the  designer  to  observe  and 
predict  the  effect  of  system  gain  changes  on  zero  locations. 
Analysis  of  the  transfer  functions  of  interest,  when  derived 
in  general  terms,  indicates  the  relative  influence  of  each 
gain  on  each  factor  of  the  characteristic  and  numerator 
equations.  The  obvious  drawback  of  this  approach,  as  with  the 
root  locus  method,  is  the  problem  of  designing  each  "half"  of 
the  two-input  system  separately  -  resulting  in  a  composite 
system  with  unknown  characteristics. 

Limited  State  Feedback  (Dual  Input) .  This  technique  is 
simply  a  reduced  form  of  full-state  feedback  with  the  added 
complexity  introduced  by  modeling  the  complete  two- input 
system  with  feedback  to  each  input.  This  method  may  or  may 
not  allow  assignment  of  eigenvalues,  depending  upon  the 
feedback  enployed.  It  does,  hcwever,  generally  allow  suffi¬ 
cient  flexibility  to  provide  acceptable  pole  locations.  In 
addition,  transfer  functions  can  be  derived  in  general  terms 
to  permit  the  assessment  of  the  impact  of  system  gain  values 
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on  zero  locations. 

The  primary  drawback  of  this  approach  is  the  sheer  volume 
of  algebraic  manipulations  that  are  required  to  solve  for  the 
transfer  functions  -  the  extent  of  effort  increasing  exponen¬ 
tially  with  the  system  order.  The  method  does,  however,  allow 
the  addition  of  cascade  and  feedback  compensation,  which 
provides  additional  flexibility  -  at  the  cost  of  increasing 
the  order  of  the  system. 

The  limited  state  feedback  approach  is  the  design 
technique  that  produced  the  best  results.  No  attempt  is  made 
to  exactly  specify  eigenvalue  locations;  instead,  a  solution 
is  sought  in  an  iterative  manner  with  the  intent  of 
determining  the  system  gains  that  result  in  poles  and  zeros 
being  placed  in  an  acceptable  envelope.  The  remaining  portion 
of  this  chapter  details  how  this  approach  is  used  in 
designing  a  flight  control  system  that  meets  the  performance 
requirements. 

Limitation  of  the  Rigid-Body  Assumption 

When  deciding  upon  a  desired  or  model  closed-loop  charac¬ 
teristic  equation,  the  placement  of  the  dominant  eigenvalues 
is  perhaps  the  most  important  consideration.  The  position  of 
these  eigenvalues  largely  determines  the  passband  of  the 
system  -  discounting  the  affect  of  zeros,  of  course.  If 
dominant  pole  placement  is  such  that  this  passband  extends 
sufficiently  high  in  frequency  to  allow  sensitivity  to 
unmodeled  body  bending  modes,  rigid-body  simulation  results 
are  no  longer  valid. 
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In  designing  a  flight  control  system  for  a  high  perfor¬ 
mance  aircraft,  there  is  a  tendency  to  increase  the  open-locp 
natural  frequency  to  permit  a  faster  response  time.  The 
farther  the  dominant  poles  are  moved  from  the  origin,  the 
faster  the  system  responds  to  an  input.  There  are  several 
physical  limits  that  must  be  observed,  however,  such  as 
aircraft  structure,  control  surface  hinge  moments,  maximum 
actuator  rates,  handling  quality,  and  sensitivity  to  body 
bending  modes.  . 

There  are  many  generally  accepted  guidelines  to  aid  in 
the  selection  of  dominant  eigenvalues  for  the  purpose  of 
attaining  reasonable  handling  qualities.  One  method,  express¬ 
ed  as  a  mapping  of  short  period  natural  frequencies  against 
damping  ratio,  is  known  as  the  Cornell  Aeronautical  Lab  (CAL) 
"Thumbprint"  (Ref  22) .  This  envelope,  shown  in  Figure  25,  was 
derived  from  extensive  handling  quality  data  obtained  with 
variable  stability  aircraft.  It  is  intended  for  use  in  the 
design  of  flight  control  systems  for  fighter  aircraft. 
Caution  must  be  exercised  in  its  use,  however,  for  if  a 
flight  control  system  is  designed  such  that  it  is  inadver¬ 
tently  sensitive  to  a  bending  mode,  the  aircraft  could  become 
unstable  under  certain  conditions.  Therefore,  it  is  important 
to  ensure  that  the  frequencies  of  the  structural  bending 
modes  are  well  outside  the  proposed  closed-loop  system 
passband . 

A  partial  list  of  structural  bending  modes  for  the  YF-16 
is  given  in  Table  I  (Ref  8).  Note  that  the  frequency  of  the 
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SHORT  PERIOD  DAMPING  RATIO 


Figure  25.  Cornell  Aeronautical  Laboratory  (CAL)  Thumbprint  (Ref  22). 


first  wing  bending  mode  is  approximately  31  radians  per 
second,  and  the  frequency  of  the  first  fuselage  bending  mode 
is  approximately  77  radians  per  second.  For  longitudinal 
analyses,  the  fuselage  bending  modes  are  generally  of  great¬ 
est  interest,  and,  in  many  instances,  the  second  structural 
modes  may  be  dominant.  For  this  aircraft,  the  second  and 
third  modes  are  dominant,  indicating  a  relatively  sound 
fuselage  (Ref  8).  From  this  data,  it  is  assumed  that 
placement  of  dominant  eigenvalues  within  the  CAL  Thumbprint 
will  not  result  in  unacceptable  sensitivity  to  any  of  the 
bending  modes. 


State-Space  Equations  for  the  Aircraft 


For  the  purposes  of  design,  a  two-degree-of-f reedom,  or 


Short  Period  model  for  the  aircraft  is  used.  This  simplifi¬ 


cation  is  valid  for  modern,  high-speed  fighter  aircraft,  as 
small  pitch  perturbations  have  a  negligible  effect  on  forward 
velocity  (Ref  1).  For  the  flight  condition  of  interest,  this 
model,  derived  in  Appendix  C,  can  be  written  as 


(24) 


The  longitudinal  quantities  that  are  measurable  and  available 
for  feedback  are  angle-of-attack,  pitch  rate,  and  normal 
acceleration.  These  outputs  must  now  be  expressed  in  terms  of 
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the  two  states  and  input  quantities  of  Eon  24.  The  simplest 
of  the  three  is  pitch  rate,  which  is  one  of  the  two  aircraft 
states.  It  may  be  written  directly  as 
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Normal  acceleration  can  be  expressed  as  a  linear  combina¬ 
tion  of  the  aircraft  states.  Referring  to  Figure  2  in  Chapter 
II,  the  sum  of  the  forces  in  the  Z-direction  may  be  expressed 
as  (Ref  1), 


m(w  -  UDq) 


(26) 


where  'm'  is  the  mass  of  the  aircraft.  Pecall  that  the  force 
acting  on  an  object  is  equal  to  the  product  of  the  mass  of 
that  object  and  its  acceleration.  If  both  sides  of  this 
equation  are  divided  by  the  mass,  an  expression  for  accelera¬ 
tion  in  the  Z-direction  results.  Setting  'hz  '  equal  to  the 
sum  of  these  accelerations,  the  expression  becomes 


A 


z 


(27) 


Referring  to  Eqn  10,  if  Uc  is  factored  from  the  right  side  of 
the  expression,  then 


Az  =  U0('a  -  q)  (28) 

defines  normal  acceleration,  positive  downward,  at  the  CG  of 
the  aircraft.  The  normal  accelerometers  in  the  YF-16, 
however,  are  not  located  at  the  CG,  but  rather  12.95  feet 
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further  forward  (see  Appendix  E  for  this  calculation).  As 
positive  aircraft  rotation  results  in  an  upward  movement  of 
these  accelerometers,  they  sense  pitching  accelerations  in 
addition  to  normal  accelerations.  The  linear  magnitude  of 
these  angular  pitch  accelerations  at  the  accelerometer 
location  is  equal  to  the  product  of  the  acceleration,  q,  and 
the  length  of  the  moment  arm,  12.95  feet.  Peversing  the  terms 
on  the  right  side  of  Eqn  28,  this  composite  measurement  - 
which  is  still  referred  to  as  normal  acceleration  -  is  now 
defined  as 

aNs  =  uo(<l  -  ’a)  +  l2-954  (29) 

and  the  acceleration  is  positive  upward. 

Equation  29  may  now  be  expanded  in  terms  of  the  states 
and  input  quantities  given  in  Eqn  24.  With  U0  =  893.5  ft/sec 
for  this  flight  condition,  this  expansion  results  in 

An  =  893.59  -  893.5  (-2.5194 'a  +  -98255q  -  .28ll5&h 

+  .3637 6f  -  2.5194  ctg)  +  12.95(12.339  ’a  -  2.48539 
-  42.317  6h  +  6.1322  6f  +  12.339  Og)  (30) 

Performing  the  multiplications  and  summing  like  terms  yields 

aNs  =  2410.9  'a  -  16.595  9  -  296.80  6h  -  245.56  6f  +  2410.9  (Xg  (31) 

From  Eqn  29,  the  dimension  of  Ans  at  this  point  is 
f t-deg/sec-sec,  as  all  angles  are  in  degrees,  and  the  forward 
velocity  is  in  feet  per  second.  As  the  normal  accelerometers 
sense  Ans  in  Gs,  it  is  appropriate  to  multiply  Eqn  31  by  some 
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conversion  factor.  Since  one  radian  is  equal  to  57.296 


degrees,  and  at  sea  level,  g  =  32.177  ft/sec-sec,  the 

dimensional  conversion  factor  is 


CF 


r _  « _ .] 

ft-deg  /  sec-sec 


=  (32.177) (57.296)  =  1843.6 


(32) 


Dividing  by  this  factor,  normal  acceleration  (in  Gs)  may  be 
written  as 


0S]  =  [1.3077  -.0090014] 

+  [-.16099  -.13320] 
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[l.3077] 
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(33) 


The  third  output  to  be  derived  is  attack  angle.  As  the 
aircraft  is  modeled  to  fly  through  other  than  still  air, 
measurable  angle-of- at  tack  is  not  the  equivalent  of  stabil¬ 
ity-axis  perturbation  attack  angle  as  given  in  Eqn  24.  The 
relationship  between  these  angles,  as  previously  derived,  is 


Om  =  'a  +  Qg 


(17) 


There  is  an  additional  complication  due  to  the  affect  of 
the  location  of  the  sensor  itself.  As  with  the  normal  accel¬ 
eration  defined  in  Eqn  28,  attack  angle  has  been  defined  as  a 
measurement  at  the  CG  of  the  aircraft.  The  angle-of-attack 
sensors  on  the  YF-16  are  located  on  the  nose,  however, 
approximately  22.28  feet  forward  of  the  CG,  as  determined  in 
Appendix  E.  Thus,  these  sensors,  which  measure  the  angle  of 
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incidence  of  the  moving  air,  are  affected  by  the  pitch  rate 
(angular  velocity)  of  the  aircraft  due  to  rotation. 

One  way  to  determine  an  expression  for  the  effect  of 
pitch  rate  on  the  sensor  is  to  view  that  motion  in  terms  of 
velocities.  The  linear  vertical  velocity  contribution  of  the 
angular  velocity  at  the  sensor  location  resulting  from 
aircraft  rotation  is  simply  the  product  of  the  magnitude  of 
that  velocity  and  the  distance  from  the  CG .  The  vertical 
velocity  of  the  sensor,  positive  downward,  relative  to  the 
earth,  may  now  be  given  as 

ws  =  w  -  22.28q  (34) 

where  pitch  rate  is  measured  in  radians  per  second.  Dividing 
each  term  in  this  expression  by  the  equilibrium  forward 
velocity,  UG  ,  results  in  an  equation  def  ning  the  angular 
measurement  of  the  sensor,  relative  to  the  earth.  This  may  be 
written  as 


Og  =  'a  -  22.28  (q  / UQ)  (35) 

where  attack  angle  is  dimensioned  in  radians.  In  a  manner 
analogous  to  the  derivation  of  Eqn  17,  the  measurment  of 
attack  angle  at  the  sensor  location  relative  to  the  surround¬ 
ing  air  can  be  expressed  as 

Og  =  'a  -  . 02494 q  +  Og  (36) 

where  the  velocity  of  the  aircraft  is  substituted  for  U0»  and 
angle-of- attack  is  again  dimensioned  in  radians.  As  pitch 
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rate  in  the  system  is  actually  measured  in  degrees  per 
second,  and  angle-cf-attack  is  measured  in  degrees,  the 
equation  is  appropriately  dimensioned  for  both  states  with  no 
need  for  dimensional  conversion.  Combining  Equations  25,  33, 
and  36,  the  system  measurable  output  equation  may  now  be 
expressed  as 


In  the  following  sections,  a  system  configuration  is 
selected  and  modeled,  and  the  transfer  functions  of  interest 
for  that  system  are  derived.  So  that  the  model  and  associated 
transfer  functions  may  be  determined  in  a  more  general  form, 
Eqns  24  and  37  are  used  with  symbolic  entries  -  requiring 
substitution  of  the  actual  values  for  the  selected  flight 
condition  (0.8  Mach  @  SL)  at  the  appropriate  time.  Use  of 
unspecified  entries  allows  the  derived  transfer  functions  to 
be  applied  to  any  flight  condition  -  or  even  any  aircraft 
that  is  modeled  in  short  period  form.  It  permits,  in  general, 
the  substitution  of  any  two-state  system,  which  includes  a 
discrete  state  model.  By  converting  a  continuous  system  model 
to  a  discrete  system  model  for  some  specified  sampling  rate 
(using  the  FGRTFAN  pregram  YFCCV, 


as  described  in  Appendix 


K),  the  technique  presented  in  this  Chapter  may  be  used  to 
accomplish  a  direct  digital  design.  The  generalized  aircraft 
state  and  output  matrix  equations  are  given  in  Eons  38  and 
39. 
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(39) 


Selected  System  Configuration 

Three  basic  system  configurations  were  investigated  in  an 
attempt  to  determine  which  would  be  most  effective  in  meeting 
the  design  criteria.  The  system  chosen  on  that  basis  is  shown 
in  Figure  26.  For  reference  purposes,  details  of  the  rejected 
system  variations  are  contained  in  Appendix  F,  along  with  a 
discussion  explaining  why  each  did  not  meet  the  design 
requirements. 

The  selected  configuration,  the  q-command  system  in 
Figure  26,  contains  cascade  compensation,  and  makes  use  of 
simple  gain  feedback  to  produce  the  desired  closed-loop 
performance.  All  gains  and  compensation  are  given  in  general 
form,  and  are  determined  in  the  design  process.  The  power 
actuators  have  been  eliminated  from  the  system  model  for  the 
purpose  of  simplification.  They  are,  however,  not  forgotten 
or  assumed  away,  as  they  are  included  in  the  system  for  all 
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simulations  presented  in  the  latter  portion  of  this  chapter. 

Each  of  the  two  servos  is  modeled  as  a  first-order  lag 
with  an  unspecified  eigenvalue.  Recall  that  the  servos  and 
actuators  each  have  a  resistive  output  available  which  can  be 
used  to  produce  a  voltage  that  is  proportional  to  position. 
As  currently  on  the  aircraft,  each  servo  is  actually  an 
integrating  device,  with  the  required  feedback  provided  by  a 
mechanical  linkage.  Removal  of  this  linkage  would  return  the 
servo  to  its  integrating  posture,  and  an  electrical  feedback 
could  close  the  loop,  modifying  the  servo  dynamics  so  as  to 
produce  a  desired  first-order  lag.  Although  this  flexibility 
exists,  the  original  dynamics  of  the  servos  are  maintained  in 
the  design  solution. 

One  important  feature  of  the  system  configuration  is  the 
washout  circuit  that  is  used  to  feed  the  trailing-edge 
flaperon  servo.  The  purpose  of  this  compensation  is  to 
prevent  a  steady-state  flaperon  deflection,  thereby 
eliminating  the  problem  of  increased  drag  in  high 
angle-of-attack  conditions.  This  technique  permits  flaperon 
deflection  during  transients  only,  such  as  in  turbulence  or 
in  response  to  pilot  commands,  but  ensures  that  the  horizon¬ 
tal  tail  is  the  primary  control  surface. 

Also  of  interest  is  the  use  of  both  a  proportional  and  an 
integral  feed-forward  path  to  the  horizontal  tail.  The 
proportional  path  allows  quick  response  to  pilot  commands, 
and  the  integral  path  ensures  zero  steady-state  error  in  the 
commanded  pitch  rate.  The  trailing-edge  flaperons  are  com- 
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manded  proportionally,  •  through  an  unspecified  gain.  Pitch 
rate  is  also  fed  back  inside  of  the  P+I  network  to  allow  for 
increased  design  flexibility.  Lastly,  an  unspecified  first- 
order  lag  command  prefilter  is  included  to  smooth  transients 
in  commands  and  improve  handling  characteristics. 

System  Equation 

Given  the  selected  configuration,  the  system  must  now  be 

modeled  as  a  single  matrix  expression.  This  can  be  accom- 

/ 

plished  in  a  fairly  straightforward  fashion,  step-by-step, 
beginning  with  the  modeled  short  period  aircraft  dynamics  and 
working  back  towards  the  system  input.  This  effort  is 
described  in  Appendix  G,  with  only  the  final  system  model 
presented  here.  The  system  state  equation  is  given  in  Egn  40, 
and  the  system  measurable  output  equation  is  given  in  Egn  41. 
Notice  that  the  first-order  command  prefilter  is  not  included 
in  the  state  equation.  This  is  because  the  prefilter  is  a 
simple  cascade  signal  conditioner  that  is  intended  to  provide 
rate-limiting  of  the  commanded  input,  and  therefore  has 
little  or  no  impact  on  the  design  of  the  remaining  system 
dynamics.  It  should  also  be  noted  that  the  attack  angle 
sensor,  pitch  rate  gyro,  and  normal  accelerometer  are  not 
modeled  in  the  state  equation.  They  have  not  been  included 
for  design  purposes  as  the  dynamics  of  these  devices  are 
essentially  non-dominant,  and  therefore  they  have  a  negligi¬ 
ble  affect  on  the  design  solution.  They  are  shown  in  Figure 
26  for  reference  only.  They  are,  however,  included  in  the 
system  model  for  all  simulations. 
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System  Transfer  Functions 

With  the  entire  system  expressed  in  a  matrix  format, 
Cramer.  *s  Rule  may  now  be  directly  applied  to  solve  for  the 
required  transfer  functions.  As  the  selected  system  is  a 
pitch  rate  control  system,  the  transfer  function  relating 
pitch  rate  to  pilot  command  is  of  great  interest.  Also, 
because  the  design  criteria  require  the  minimization  of  pitch 

response  to  disturbance  inputs,  the  transfer  function  relat- 

/ 

ing  pitch  rate  to  a  vertical  gust  input  is  needed. 

Details  of  the  process  of  deriving  the  transfer  functions 
are  contained  in  Appendix  H,  where  the  final  product,  given 
in  general  form,  is  the  system  characteristic  equation  and 
the  two  numerator  equations.  As  the  particular  equations  for 
the  chosen  flight  condition  are  now  required,  the  values  of 
the  entries  in  Eqns  24  and  37  are  substituted  into  them.  The 
resulting  expression  for  the  closed-loop  system 
characteristic  equation  is  given  in  Eqn  42.  The  numerator 
expression  for  pitch  rate  output  with  a  pilot  command  input 
is  given  in  Equation  43,  and,  the  numerator  expression  for 
pitch  rate  output  with  a  disturbance  input  is  given  in 
Equation  44. 


Normally,  when  designing  a  state  variable  feedback 
system,  a  desired  system  characteristic  polynomial  is  deter¬ 
mined  and  set  equal  to  the  derived  system  characteristic 
polynomial  on  a  term  by  term  basis.  In  this  way,  each  feed¬ 
back  gain  is  resolved  to  a  particular  value  such  that  the 
eigenvalues  of  the  system  being  designed  are  set  equal  to 


[  5.0047  +  (p  +  Y  +  R)  +  2.664KfKg  ]  ss  +  [  -5.8622  +  5. 0047  (  0  +  y  +  R)  + 


[-5.8622R^y  +  2201. ?(Rp  +  z0  +  Rz)KhKcK5  +  846 ^RpzK^C^  +  220i.7RpKhK3  +  645.55RpKhKi  - 
2201.7DRpKhK1  +  . 33393zKhKfKCK6K5  "  5466.6zKhKfKcK5K2 ]  6  +  [  2201.7zRpKhK(J<5] 


[846.34KhKc  -  i22.64KfKp] 


those  of  the  desired  or  model  system.  For  a  single  input 
system,  this  is  relatively  straightforward,  as  the  number  of 
feedback  loops  is  generally  equal  to  the  order  of  the  system. 
In  the  case  of  a  two  input  system  with  non-independent 
outputs,  however,  a  multiple  of  each  system  state  may  be  fed 
back  two  or  more  times  -  necessitating  the  solution  of  non¬ 
linear  equations  to  determine  the  feedback  gain  values. 

For  the  system  proposed  in  this  chapter,  five  feedback 
loops  are  closed  around  a  fifbh-order  system.  As  three  eigen¬ 
values  are  modeled  as  completely  variable,  the  feedback  loops 
are  used  to  modify  only  the  second-order  dynamics  of  the 
airframe.  With  the  addition  of  a  sixth  feedback  loop  to 
define  the  controlled  state,  the  resulting  characteristic 
equation  for  this  configuration  contains  highly  nonlinear 
polynomial  coefficients.  A  glance  at  Eqn  42  confirms  this 
statement.  Knowing  that  the  equation  is  sixth  order,  and 
that  there  are  thirteen  variables,  only  six  of  these  vari¬ 
ables  can  be  independent  -  leaving  the  remaining  seven  as 
dependent  variables.  There  is  the  possibility  that  more  than 
six  degrees  of  freedom  exist,  implying  a  potential  capability 
to  solve  for  a  desired  characteristic  polynomial  and  numera¬ 
tor  polynomial  simultaneously.  This  approach  was  attempted 
for  the  numerator  polynomial  relating  pitch  rate  output  to  a 
vertical  gust  input,  but  no  way  was  found  to  solve  the 
resulting  set  of  nine  nonlinear  equations. 

When  designing  a  system  to  reject  a  certain  input  (a 
vertical  gust  disturbance,  in  this  case),  the  numerator 
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equation  is  very  important.  Even  in  this  system,  where  the 
numerator  gain  is  fixed  (see  Eqn  44),  the  placement  of  zeros 
relative  to  poles  can  cause  large  variations  in  the  system 
response  to  that  input.  The  mere  placement  of  eigenvalues  is 
not  sufficient.  The  design  is  further  complicated  by  the  need 
to  ensure  that  the  system  responds  in  a  desired  manner  for  a 
second  (pilot)  input.  An  attempt  to  assign  locations  for  the 
zeros  of  both  numerator  polynomials  and  the  poles  of  the 
characteristic  polynomial,  however,  would  certainly  require 
more  degrees  of  freedom  than  are  available.  Therefore,  even 
if  a  technique  for  handling  the  resulting  nonlinear  equations 
could  be  found,  there  would  be  no  solution. 

Assuming  that  it  is  not  possible  to  specify  a  desired 
location  for  both  the  poles  and  zeros  of  both  transfer  func¬ 
tions,  it  is  still  very  possible  that  a  combination  of  system 
gains  can  be  found  that  results  in  acceptable  system 
performance.  Given  the  proposed  closed-loop  system,  the 
flexibility  afforded  by  the  number  of  system  variables  makes 
the  existence  of  such  a  solution  quite  probable. 

A  tr ial-and-error  approach  is  utilized  to  search  for  an 
acceptable  design  solution.  This  process  is  simplified 
considerably  by  the  use  of  an  interactive  computer  program. 
Designated  YFSOLVE ,  the  program  allows  the  user  to  input  a 
set  of  system  variables,  and  then  calculates  the  closed-loop 
transfer  functions  of  interest  in  both  polynomial  and  factor¬ 
ed  form.  One  or  more  of  the  system  variables  can  then  be 
altered;  permitting  the  user  to  note  the  corresponding  change 
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in  each  of  the  three  equations.  This  method,  however, 
requires  each  resulting  set  of  equations  to  be  evaluated  to 
determine  if  the  current  iteration  produces  a  better  system 
than  the  previous  iteration.  With  the  aid  of  program  YFCCV, 
several  system  solutions  were  simulated  in  the  process  of 
developing  the  ability  to  recognize  a  good  design  solution 
from  the  transfer  functions. 


System  Design  Solutions 

/ 

In  the  remaining  portion  of  this  chapter,  two  design 
solutions  are  analyzed  to  determine  their  respective 
performance  characteristics.  Both  solutions  use  the  system 
configuration  shown  in  Figure  26.  The  primary  design,  which 
is  called  the  Pitch  Performance  Enhancement  (PPE)  flight 
control  system,  is  designed  to  minimize  pitch  response  to 
vertical  gusting  -  using  both  the  horizontal  tail  and  the 
tr ailing-edge  flaperons  to  nullify  pitch  perturbations.  As  a 
result,  the  flaperons  initially  deflect  in  a  positive  direc¬ 
tion  (down)  in  response  to  an  upward  vertical  gust,  and 
assist  the  horizontal  tail  in  controlling  pitching  moment. 

The  second  design,  designated  YFF57 ,  is  included  for 
comparison  purposes.  This  alternate  system  makes  use  of  a 
different  set  of  system  gains  to  allow  a  negative  flaperon 
deflection  in  response  to  an  'upward  vertical  gust.  This  is 
done  in  an  attempt  to  reduce  normal  acceleration  while 
simultaneously  minimizing  pitch  response  with  horizontal  tail 
deflections.  It  is  shown  that  for  the  pitch  rate  command 
system  proposed,  only  a  small  improvement  in  normal  accel- 


eration  can  be  realised,  and  only  at  the  expense  of  largely 
degrading  pitch  response. 

The  values  determined  for  the  system  variables  for  the 
PPE  and  YFF57  flight  control  systems  are  given  in  Table  2. 
Note  that  the  feedback  and  feed-forward  gains  are  largely  the 
same.  The  differences  lie  in  the  flaperon  forward-loop  gain, 
the  relative  amount  of  attack  angle  feedback  to  the  flaper- 
ons,  and  the  amount  of  pilot  input  driving  the  flaperon 
channel.  In  each  case,  the  eigenvalues  of  both  servos  are 
unchanged  from  those  of  the  existing  servos  on  the  YF-16.  As 
only  slight  improvement  is  obtained  by  modifying  them,  the 
original  eigenvalues  are  retained  so  as  not  to  require  any 
mechanical  modifications  to  the  aircraft.  Also,  for  both 
designs,  an  attempt  is  made  to  keep  the  value  of  each 
variable  within  reasonable  bounds;  i.e.,  no  extremely  high 
feedback  or  feed-forward  gains. 

The  transfer  functions  for  each  of  these  systems,  as 
output  by  program  YFSOLVE,  are  given  in  Tables  3  and  4. 
Examining  the  transfer  functions  modeling  pitch  rate  for  a 
pilot  input  (Q/P),  recall  that  the  objective  in  each  case  is, 
of  course,  to  make  the  aircraft  flyable  as  well  as  insensi¬ 
tive  to  disturbances.  The  dominant  eigenvalues  of  the  PPE 
system  are  located  to  be  overdamped,  making  use  of  two 
dominant  zeros  and  a  second  pair  of  underdamped,  but  less 
dominant  poles  to  enhance  the  initial  transient  response. 
Acceptable  response  characteristics  are  obtained  with  system 
YFfi7,  on  the  other  hand,  by  placing  dominant  poles  and  zeros 
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in  close  proximity,  with  overall  response  also  affected  by  a 
pair  of  less  dominant  poles. 

Examination  of  the  pitch  rate  for  disturbance  input  (Q/W) 
transfer  functions  reveals  that  the  system  configuration 
results  in  two  free  differentiations.  This  permits  zero 
steady-state  pitch  angle  perturbation  for  a  step  disturbance 
input.  In  each  case,  an  attempt  is  made  to  locate  the 
remaining  three  zeros  in  such  a  way  as  to  minimize  the  tran¬ 
sient  output  from  the  transfer  functions.  It  should  not  be 
forgotten  that  freedom  to  locate  poles  and  zeros  in  specific 
locations  is  limited.  As  there  are  a  total  of  thirteen 
alterable  poles  and  zeros,  at  least  an  equal  number  of 
degrees  of  freedom  must  exist  to  freely  place  them.  The 
selection  of  transfer  functions  is  therefore  not  only  con¬ 
strained  in  this  manner,  but  also  by  the  tr ial-and-er ror 
technique  that  must  be  employed. 

System  Time  Responses 

Using  the  values  determined  for  the  system  variables, 
time  response  simulations  are  run  for  each  of  the  two 
designs.  As  in  the  previous  chapter,  simulation  data  is  in 
the  form  of  computer  generated  plots.  As  the  primary  design 
solution,  the  PPE  flight  control  system,  is  of  greater 
interest,  simulations  are  run  for  both  a  discrete  vertical 
gust  input,  and  a  pilot  step  input.  The  alternate  design 
solution,  YFF57 ,  is  included  only  to  show  the  effect  of 
requiring  the  trailing-edge  flaperons  to  deflect  so  as  to 
oppose  vertical  translation  of  the  aircraft  in  response  to 
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vertical  gusting.  Therefore,  this  system  is  simulated  for  a 
discrete  vertical  gust  input  only. 

Response  to  Pilot  Step  Input.  To  simulate  system  response 
to  a  pilot  command,  a  one  degree  per  second  step  input  is 
applied.  The  resulting  aircraft  response  (PPE  system)  is 
shown  in  Figure  27.  Pilot  stick  input,  as  indicated,  is 
dimensioned  in  volts  to  account  for  a  general  flight  control 
system,  as  previously  discussed.  The  applied  unit  step  input 
is  actually  dimensioned  in  degrees  per  second,  as  the  system 
is  a  pitch  rate  FCS. 

The  simulation  results  shewn  in  Figure  27  may  now  be 
compared  to  the  simulation  results  for  the  YF-16  flight 
control  systems  presented  in  the  last  chapter  in  Figures  12 
and  13.  In  making  this  comparison,  it  becomes  apparent  that 
the  systems  are  not  responding  to  the  same  magnitude  of 
input.  As  the  YF-16  flight  control  systems  are  G-command 
systems,  a  unit  step  input  produces  a  1-G  normal  acceleration 


output. 

and 

roughly  a  2.05  degree  per 

second 

pitch  rate 

output. 

For 

a  pitch  rate  command 

system 

,  a  unit 

step  input 

produces 

a 

one  degree  per  second 

p  i  ten 

rate  output,  and  in 

this  case,  roughly  a  0.485  G  normal  acceleration  output. 
Therefore,  the  commanded  system  responses  essentially  differ 
by  a  factor  of  two. 

One  way  of  judging  system  performance  from  time  response 
data  is  to  determine  some  of  the  classical  response  charac¬ 
teristics,  such  as  rise  time,  peak  time,  and  settling  time. 
These  measurements  generally  provide  a  good  overall  indica- 


tion  of  system  transient  performance.  The  calculated  time 
response  characteristics  for  the  measurable  aircraft  outputs 
are  given  in  Table  5.  Rise  time  is  computed  in  two  ways:  the 
period  of  time  during  which  the  output  increases  from  5%  to 
80%,  and,  similarly,  from  5%  to  95%  of  its  maximum  value.  The 
former,  less  standard  measure  is  made  necessary  by  the 
overdamped  nature  of  the  current  YF-16  flight  control 
systems. 

Attack  angle  response,  in  this  instance,  is  a  measure  of 
how  quickly  the  aircraft  attains  the  required  attitude  with 
respect  to  the  relative  wind.  At  a  given  flight  condition,  an 
increase  in  pitch  rate  or  normal  acceleration  requires  a 
specific  increase  in  attack  angle.  This  is  only  true,  of 
course,  for  commanded  maneuvers  that  are  achieved  as  a  result 
of  aircraft  rotation,  as  is  the  case  for  each  of  these 
systems.  It  is  true  that  the  Maneuver  Enhancement  system  and 
the  PPE  FCS  use  the  direct  lift  created  by  the  flaperons  to 
enhance  normal  acceleration  response.  The  effect  of  such 
flaperon  deflection,  however,  is  limited  to  the  early  tran¬ 
sient  portion  of  the  commanded  response,  as  discussed  in 
Chapter  III.  The  improvement  in  attack  angle  response  realiz¬ 
ed  by  the  PPE  design  is  quite  apparent  from  the  table.  Rise 
time  is  much  faster  than  the  Maneuver  Enhancement  system, 
which  is,  surprisingly,  even  slower  than  the  baseline  FCS. 

Pitch  response  is  also  much  improved  -  especially  beyond 
80%  of  final  value.  While  the  baseline  and  Maneuver  Enhance¬ 
ment  system  responses  slowly  approach  a  maximum,  the  PPE  FCS 


damped  to  prevent  measurement  of  rise  time. 


output  settles  within  2%  in  three-quarters  of  a  second.  The 
amount  of  overshoot  that  occurs  is  well  within  an  acceptable 
range  for  pitch  rate  command  systems,  and  is  therefore  not  a 
problem. 

Comparing  normal  acceleration  responses,  it  is  interest¬ 
ing  to  note  that  the  baseline  system  generates  a  faster 
response  than  the  Maneuver  Enhancement  system  -  at  least 
above  about  80%  of  the  commanded  value.  The  overall  PPE 
system  response  is,  however,  more  than  three  times  as  rapid 
as  the  baseline  system.  In  addition,  the  transition  from 
translation  to  rotation,  which  occurs  at  about  the  0.25 
second  point,  is  much  more  subtle.  Overall  performance  is, 
clearly,  much  improved. 

Response  to  the  Discrete  Vertical  Gust.  To  simulate  sys¬ 
tem  response  to  atmospheric  disturbances,  discrete  (1-cos) 
vertical  gusts  are  applied.  To  enable  a  direct  comparison  to 
the  simulation  results  in  Chapter  III,  the  discrete  gust  peak 
velocity  is  again  50  feet  per  second,  with  a  duration  of  one 
second . 

Simulation  results  for  the  PPE  system  are  shown  in  Figure 
28,  and  the  results  for  system  YFF57  are  shown  in  Figure  29. 
As  was  intended,  the  trailing-edge  flaperons  deflect  in 
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opposite  directions  for  each  system.  As  discussed,  the  PPE 
system  is  designed  such  that  the  flaperons  augment  the 
horizontal  tail  in  opposing  pitch  moment,  and  therefore,  as 
both  control  surfaces  are  located  aft  of  the  CG,  both 

4 

initially  deflect  in  the  same  direction. 
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Comparing  the  normal  acceleration  plots  for  both  systems, 
it  is  clear  that  system  YFF57  has  resulted  in  only  about  a  7% 
improvement  -  a  reduction  from  2.48  to  2.32  Gs  perturbation. 
At  the  same  time,  however,  pitch  angle  response  is  degraded 
by  about  470%  -  from  0.017  to  0.08  degrees  perturbation. 
Although  this  pitch  response  is  still  very  acceptable,  the 
small  improvement  afforded  in  normal  acceleration  response 
with  this  system  leads  to  the  selection  of  the  PPE  FCS  as 
more  appropriate  relative  to  the  design  objective. 

Comparing  the  simulation  results  of  the  PPE  FCS  (Figure 
28)  to  the  results  obtained  with  the  YF-16  flight  control 
systems  (Figures  15  and  17)  in  the  previous  chapter,  the 
differences  are  readily  determined.  The  most  obvious  is  the 
desirable  improvement  in  pitch  response  and  tracking  error. 
Peak  pitch  angle  perturbation  is  reduced  by  a  factor  of  about 
140  over  the  baseline  YF-16  FCS,  and  by  a  factor  of  about  130 
over  the  YF-16  FCS  with  Maneuver  Enhancement.  Maximum  longi¬ 
tudinal  tracking  error  is  reduced  by  factors  of  about  9.5  and 
9.0,  respectively. 

A  greater  relative  improvement  in  tracking  error  can  be 
realized  with  a  target  at  increased  distance.  This  is  because 
the  tracking  error  for  the  PPE  system  results  mostly  from  the 
vertical  movement  of  the  aircraft  -  not  the  pitch  perturba¬ 
tion.  Exactly  the  opposite  is  true  for  the  YF-16  flight 
control  systems.  Equally  important  is  the  shape  of  the  actual 
tracking  error  response.  It  is  very  probable  that  the  track¬ 
ing  error  with  the  PPE  system  could  be  more  easily  controlled 


by  the  pilot 


Normal  acceleration  response,  on  the  other  hand,  is 
improved  less  significantly  compared  to  the  baseline  YF-16 
PCS.  In  Figure  15,  it  is  shown  that  the  YF-16  with  baseline 
FCS  responds  by  pulling  roughly  2.0  positive  Gs,  followed  by 
about  2.7  negative  Gs.  This  results  in  a  peak-to-peak  change 
of  4.7  Gs  over  a  one-half  second  time  period.  Also,  negative 
Gs  are  more  difficult  to  withstand  than  positive  Gs.  The  PPE 
FCS  responds  with  a  peak-to-peak  change  of  about  4.2  Gs  over 
the  same  time  period  -  an  improvement  of  approximately  10%. 
In  addition,  the  peak  negative  perturbation  is  improved  by 
more  than  35%. 

The  normal  acceleration  response  of  the  YF-16  FCS  with 
Maneuver  Enhancement,  however,  is  substantially  better  than 
the  PPE  FCS.  The  .peak-to-peak  variation  is  only  2.4  Gs  -  a 
reduction  of  almost  50%  over  the  baseline  YF-16  FCS  and  more 
than  40%  over  the  PPE  system.  Also,  the  peak  perturbation  in 
negative  Gs  is  reduced  by  roughly  55%  over  the  baseline 
system,  and  more  than  30%  over  the  PPE  FCS.  This  is  expected, 
however,  since  the  Maneuver  Enhancement  system  is  designed  to 
control  normal  acceleration. 

Comparing  control  surface  deflections,  peak  horizontal 
tail  movement  is  reduced  by  more  than  70%  over  the  baseline 
YF-16  FCS,  and  about  75%  over  the  YF-16  FCS  with  ME.  In 
addition,  maximum  horizontal  tail  deflection  rate  is  reduced 
by  roughly  75%  and  80%,  respectively.  Trailing-edge  flaperon 
deflection  and  rate  of  deflection  are  reduced  by  90%  and  92% 
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over  the  YF-16  FCS  with  Maneuver  Enhancement. 

To  ensure  a  valid  test  of  the  ability  of  this  system  to 
reject  discrete  disturbances,  it  is  of  some  interest  to 
excite  the  system  at  its  pitch  resonant  frequency  and  observe 
the  results.  A  gust  of  approximately  0.3  seconds  duration  was 
determined  to  be  this  worst  case  input.  The  resultant  system 
response  is  shown  in  Figure  30.  Peak  pitch  angle  perturbation 
is  still  less  than  0.06  degrees,  and  maximum  target  tracking 
error  is  reduced  to  less  thab  eight  feet.  In  addition,  the 
probability  of  a  gust  of  such  short  duration  occurring  with  a 
magnitude  equal  to  that  of  the  longer  duration  discrete  gust 
used  for  the  previous  analysis  is  quite  low  (Ref  4) .  In  light 
of  this  fact,  system  performance  at  the  pitch  resonant 
frequency  is  very  good. 

System  Frequency  Response 

To  complete  the  assessment  of  system  performance  ,  both 
analytical  and  empirical  frequency  response  analyses  are 
presented.  The  analytical  method  is  the  standard  process  of 
applying  sine  wave  inputs  of  various  frequencies  to  the 
system,  plotting  output  phase  and  magnitude  relative  to  the 
input.  The  empirical  method,  as  described  in  Chapter  III, 
indicates  the  maximum  perturbation  of  system  states  for 
constant  magnitude  discrete  gust  inputs  of  various  durations. 
This  analysis  is  included  to  permit  a  further  comparison  of 
the  PPE  and  YF-16  flight  control  systems. 


Analytical  Frequency  Response.  In  the  design  process. 


the  system  was  modeled  in  standard  matrix  format  (Fqns  40  and 
41).  These  matrices  were  used  directly  by  program  TOTAL  (Pef 
23)  to  calculate  the  analytical  frequency  response  plots 
shown  in  Figures  31  through  34.  Notice  that  the  resonant 
frequency  of  the  system  is  approximately  23  radians  per 
second.  The  system  is  therefore  maximally  excited  by  an  input 
of  that  frequency,  roughly  equivalent  to  a  0.3  second  dura¬ 
tion  discrete  gust. 

The  system  pitch  rate  Response  to  a  pilot  command  is 
essentially  flat  (within  3dB)  up  to  50  radians  per  second. 
Normal  acceleration  response  is  fairly  good  (within  6dB)  up 
to  30  radians  per  second.  The  required  deflection  of  control 
surfaces  decreases  monotonically  below  the  system  resonant 
frequency  to  maintain  system  performance  without  actuator 
displacement  saturation. 

The  system  pitch  angle  response  to  a  disturbance  input  is 
no  less  than  34  dE  down  over  the  entire  frequency  spectrum. 
This  accounts  for  the  excellent  disturbance  rejection  proper¬ 
ties  of  this  system.  Normal  acceleration  response  ,  on  the 
other  hand,  falls  off  at  20dB  per  decade  only  below  about  3 
radians  per  second.  Above  that  frequency,  there  is  little  or 
no  rejection.  Again,  the  required  deflection  of  control 
surfaces  decreases  below  the  system  resonant  frequency,  which 
reduces  the  possibility  of  actuator  displacement  saturation. 

Empirical  Frequency  Response.  Using  the  technique  de¬ 
scribed  in  the  previous  chapter,  empirical  frequency  response 
data  is  compiled  for  the  system.  Peak  gust  magnitude  is  held 
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constant  at  50  feet  per  second.  The  resulting  peak  perturba¬ 
tion  data  is  plotted  in  Figures  35  through  41.  To  permit  an 
instant  comparison  of  this  information  with  that  previously 
determined  for  the  YF-16  flight  control  systems,  the  data  for 
those  systems  is  also  plotted  in  each  figure. 

The  ability  of  the  system  to  minimize  pitch  perturbations 
resulting  from  atmospheric  disturbances  is  depicted  quite 
clearly  in  Figure  35.  Peak  perturbation  is  no  more  than  0.06 
degree  for  a  gust  of  any  frequency.  Maximum  tracking  error  is 
shown  in  Figure  36.  A  marked  improvement  is  obtained  for  all 
gust  frequencies  above  0.6  radians  per  second.  The  large 
tracking  error  associated  with  the  lower  frequencies  is  a 
result  of  insufficient  control  of  normal  position.  Because 
normal  acceleration  is  not  fed  back  to  the  horizontal  tail 
channel,  this  control  surface  does  not  directly  control 
vertical  translation.  Also,  the  washout  circuit  in  the 
tr ailing-e dge  flaperon  channel  zeros  flaperon  deflection 
prior  to  subsidence  of  gusting  at  the  lower  frequencies.  It 
should  be  well  within  the  pilot's  capabilities,  however,  to 
control  vertical  translation  under  these  conditions. 

Figure  37  verifies  that  the  PPE  system  provides  little  or 
no  improvement  over  the  baseline  YF-16  FCS  in  controlling 
normal  acceleration  peak  perturbation.  It  does  not  show, 
however,  the  improvement  in  peak-to-peak  acceleration  change, 
or  the  reduction  in  peak  negative  Gs  as  discussed  previously. 
Overall  performance  is,  therefore,  somewhat  better  than  is 
indicated  by  the  figure. 
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Reduction  in  control  surface  deflection  is  confirmed  in 
Figures  38  and  40.  Control  surface  deflection  rate  is  also 
decreased,  except  at  very  high  frequencies,  as  shown  in 
Figures  39  and  41.  Actuator  rate  for  the  horizontal  tail 
never  exceeds  18  degrees  per  second.  The  trailing-edge 
flaperon  actuator  rate-satura tes  at  a  gust  frequency  of  35 
radians  per  second,  as  opposed  to  8  radians  per  second  for 
the  YF-16  Maneuver  Enhancement  system. 


PILOT  STICK  INPUT  .  OUTPUT  OF  ATTACK  ANOLE  6EN80R 


Figure  27a.  FPE  FCS  response  to  a  pilot  step  input. 
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OUST  DISTURBANCE  INPUT  _  OUTPUT  OF  RTTRCK  ANGLE  SENSOR 


Figure  28a.  FPE  FCS  response  to  a 


Figure  28b.  FPE  PCS  response  to  a  discrete  vertical  gust. 


0U8T  DISTURBANCE  INPUT  .  OUTPUT  OF  ATTACK  ANOLE  SENSOR 


29a.  YFF57  FCS  response  to  a  discrete  vertical  gust. 


Figure  29b.  YFF57  FCS  response  to  a  discrete  vertical  gust. 


30b.  FPE  FCS  response  to  a  discrete  gust  at  the  pitch  resonant  frequency 


WKWITUOE  (DECIBELS)  PHASE  SHIFT  (OEOREE3) 


Figure  3ib.  Frequency  response  of  the  FPE  FCS. 


PHASE  SHIFT 


NORMAL  ACCELERATION  OUTPUT,  PILOT  INPUT 


NORMAL  ACCELERATION  OUTPUT.  PILOT  INPUT 


nflONITUOE 


124 


Figure  33a.  Frequency  response  of  the  PPE  FCS 


nnONtTUOE  ( OECIBELS ) 


Figure  36.  Comparison  of  peak  target  tracking  error  response  of  YF-16  CCV  and  FPE  flight 
control  systems. 


Figure  37.  Comparison  of  peak  normal  acceleration  response  of  YF-16  CCV  and  PPE  flight  control 


Comparison  of  peak  horizontal  tail  deflection  of  YF-16  CCV  and  PPE  flight  control 


Figure  39.  Comparison  of  peak  horizontal  tail  deflection  rate  of  YF-16  CCV  and  PPE  flight 
control  systems. 


V.  Conclusions  and  Fecomme ndat ions 


Conclusions 

The  primary  objective  of  this  investigation  was  to  deter¬ 
mine  the  possibility  of  reducing  the  magnitude  of  disturbance 
induced  perturbations  in  the  pitch  response  of  the  YF-16 
through  simultaneous  and  cccrdinated  activation  of  the  flap- 
erons  and  horizontal  tail.  It  has  been  clearly  demonstrated 
that  by  careful  selection  of  /the  system  configuration,  it  is 
possible  to  both  decrease  aircraft  response  to  disturbances 
and  improve  response  to  pilot  commands. 

An  iterative  approach  using  limited  state  feedback  was 


found 

to  produce  the 

best 

results,  as  this  method 

permitted 

the 

required 

con  trol 

of 

zero  locations 

by  clearly 

modeling 

the 

effect  of 

each  of 

the 

system  gains. 

Primary  feedback  to 

the 

horizon  tal 

tail 

consisted  of  pitch 

rate  and 

angle-of- 

attack  measurements.  A  secondary  pitch  rate  feedback  loop  was 
closed  around  a  proportional  plus  integral  (P+I)  compensator 
to  form  the  desired  pitch  rate  control  system.  The  use  of  the 
cascaded  P+I  compensation  resulted  in  an  excellent  system 
transient  response  to  pilot  commands,  in  addition  to  func¬ 
tioning  as  a  series  trim.  Feedback  to  the  trailing-edge 
flaperons  consisted  of  pitch  rate,  normal  accleration,  and 
angle-of- at  tack.  The  problem  of  unnecessary  drag  induced  by 
residual  flaperon  deflection  was  overcome  by  placing  a 
washout  circuit  in  the  forward  loop. 

In  the  process  of  chosing  a  system  configuration,  it  was 


observed  that  in  order  to  reduce  the  effect  of  disturbances 


on  pitch  response,  a  pitch  rate  command  system  is  required. 
The  use  of  normal  acceleration  feedback  to  the  horizontal 
tail  to  permit  a  G-command  system  clearly  degraded  the  pitch 
response  to  turbulence.  As  a  possible  compromise  solution,  a 
C*  command  configuration  was  attempted.  The  results,  however, 
were  no  better  than  with  the  pure  G-command  system. 

One  other  outcome  of  this  thesis  that  should  be  mentioned 
is  the  recognition  that  stability-axis  perturbation  attack 
angle,  'a  ,  only  has  meaning  as  an  attack  angle  when  the 
surrounding  air  is  assumed  motionless.  If  atmospheric 
disturbances  are  present,  'a  is  not  the  attack  angle  that  is 
measured  by  the  angle-of-a  t  tack  sensors.  Since  attack  angle 
is  used  as  a  feedback  quantity,  system  response  is  impacted 
by  the  manner  in  which  it  is  defined.  This  was  demonstrated 
clearly  in  Chaptf'  III. 

Re commend at  ions 

There  is  much  follow-on  research  related  to  this  investi¬ 
gation  that  can  be  of  great  interest  to  pursue.  As  with  any 
flight  control  system  design  that  has  been  examined  at  a 
single  flight  condition,  it  is  important  to  determine  if  the 
desired  performance  can  be  achieved  throughout  the  entire 
flight  envelope.  Non-dimensional  stability  derivatives  for 
several  flight  conditions  have  been  determined  for  the 
purpose  of  continued  research,  and  are  included  in  Appendix 
B.  Such  an  analysis  would  additionally  provide  the  necessary 
data  to  establish  any  required  gain  scheduling  as  a  fuction 
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of  flight  condition.  It  rray  also  reveal  that  a  particular  set 
of  fixed  system  gains  could  be  sufficient,  with  no  need  for 
gain  scheduling. 

A  related  area  that  should  be  investigated  is  the  affect 
of  changes  in  feedback  and  feed-forward  gains  on  system 
stability  and  overall  performance.  A  sensitivity  analysis  and 
failure  modes  and  effects  analysis  could  provide  a  good 
indication  of  the  magnitude  of  performance  degradation  that 
can  be  anticipated  as  a  result  of  unexpected  gain  changes. 
The  results  of  such  analyses  might  specify  a  requirement  for 
certain  components  to  be  redundant  with  either  simultaneous 
balanced  feed  or  a  voting  capability. 

For  the  purpose  of  reducing  the  amount  of  effort  required 
to  solve  for  the  system  transfer  functions,  the  configuration 
of  the  PPE  flight  control  system  did  not  include  the  second- 
order  power  actuators.  Although  they  were  included  for 
simulation  purposes,  thus  validating  the  assumption  that  they 
can  be  removed  for  the  design  process,  it  may  be  interesting 
to  determine  if  the  dynamics  of  the  actuators  can  be  used  to 
some  advantage.  To  support  such  an  effort,  a  first-order 
observer  has  been  designed  to  access  the  actuator  internal 
state,  and  is  included  in  Appendix  I. 

One  final  area  for  continued  research  might  be  to  inves¬ 
tigate  the  use  of  the  presented  approach  in  the  design  of  a 
discrete-time  digital  controller  for  the  aircraft.  In  Chapter 
IV,  the  system  equations  were  intentionally  developed  in  a 
general  form  to  support  such  an  effort.  It  is  assumed  that  a 


digital  design  can  be  implemented  in  a  manner  analogous  to 
that  presented,  by  specifying  a  sampling  rate  and  simply 
determining  transfer  functions  in  the  z-plane  instead  of  the 
s-plane.  A  successful  design  would  have  the  advantage  of 
permitting  the  use  of  a  digital  flight  control  computer.  Such 
a  system  is  not  only  generally  smaller  and  lighter  than  its 
analogue  counterpart,  but  more  flexible  and  less  failure 
prone  as  well.  Control  law  implementation  can  be  software 
derived,  as  opposed  to  harow/ired,  providing  the  capability 
for  reprogramming. 
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Appendix  A 


List  of  Abbreviations  and  Symbols 


Abbreviations 


AC 

A/C 

AFCC 

CAL 

CCV 

CDC 

CF 


CG 

CSA 

FEfc 

FC 

FCC 

FCS 

f  pS 

FS 

GK 

HT 

Hz 

IC 

L/ll 

MAC 

ME 

P+1 


Aerodynamic  Center 
Aircraft 

Aux  Flight  Control  Computer 

Cornell  Aeronautical  Laboratory 

/ 

Control  Configured  Vehicle 

Control  Data  Corporation 

Conversion  Factor 

Center  of  Gravity 

Ccmpensator/Se  r vo/Ac tua  tor 

Fly-by-wire 

Flight  Condition 

Flight  Control  Computer 

Flight  Control  System 

Feet  Per  Second 

Fuselage  Station 

Gross  Ke  ight 

Horizontal  Tail 

Hertz 

Initial  Condition 
Left  Hand 

Mean  Aerodynamic  Chord 
Maneuver  Enhancement 
Proportional  Plus  Integral 
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PPE 

- 

Pitch  Performance  Enhancement 

R/H 

- 

Right  Hand 

C/S 

- 

Radians  per  Second 

SL 

- 

See  Level 

SS 

- 

Span  Station 

TEF 

- 

Trailing-edge  Flapercr. 

TF 

Transfer  Function  (Laplace) 

Undef 

- 

Undefined 

Unk 

- 

Unknown  / 

WL 

Waterline 

u 


o 
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Symbols 


A 

AN 


B 

C 


c 


CL 

% 

Cma 

% 

C% 


Continuous- t ime  system  dynamics  matrix 

Normal  Acceleration,  defined  positive  upward, 
measured  at  the  normal  accelerometer  location 

Normal  Accelerometer  output  (positive  upward) 

Normal  Acceleration,  defined  positive  downward, 
at  the  center  of  gravity  of  the  aircraft 

Continuous-time  system  control  matrix 

System  output  matrix 

/ 

Mean  aerodynamic  chord  length 
Coefficient  of  drag 

The  change  in  the  forces  in  the  X  direction  as 
a  result  of  the  sum  of  the  deflections  of  all 
control  surfaces 

The  change  in  the  forces  in  the  Z  direction  as 
a  result  of  the  sum  of  the  deflections  of  all 
control  surfaces 

Coefficient  of  lift 

The  change  in  the  coefficient  of  lift  as  a 

result  of  a  pitching  velocity 

The  change  in  the  coefficient  of  lift  as  a 

result  of  a  change  in  angle-of-attack  caused  by 
a  change  in  w 

The  change  in  pitching  moment  as  a  result  of 
the  sum  of  the  deflections  of  all  control 
sur  f aces 

The  change  in  pitching  moment  as  a  result  of  a 
pitching  velocity 

The  change  in  pitching  moment  as  a  result  of  a 
change  in  angle-of-attack  caused  by  a  change  in 
w 

The  change  in  pitching  moment  as  a  result  of 
the  time  rate  of  change  of  angle-of-attack 
caused  by  w 
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The  change  in  the  pitching  moment  due  to  the 
deflection  of  a  control  surface 


^tef 


The  change  in  the  pitching  moment  due  to  the 
deflection  of  the  horizontal  tail 

The  change  in  the  pitching  moment  due  to  the 
deflection  of  the  trailing-edge  flaperons 

Coefficient  due  to  the  acceleration  of  gravity 

The  change  in  the  force  in  the  X  direction  as  a 
result  of  a  change  in  forward  velocity 

The  change  in  the  force  in  the  X  direction  as  a 
result  of  a  change  in  angle-of-at tack  caused  by 
a  change  in  w  , 

The  change  in  the  force  in  the  X  direction  due 
to  a  control  surface  deflection 

The  change  in  the  force  in  the  X  direction  due 
to  a  horizontal  tail  deflection 

The  change  in  the  force  in  the  X  direction  due 
to  a  trailing-edge  flaperon  deflection 

The  change  in  the  force  in  the  Z  direction  as  a 
result  of  a  change  in  forward  velocity 

The  change  in  the  force  in  the  Z  direction  as  a 
result  of  a  pitching  velocity 

The  change  in  the  force  in  the  Z  direction  as  a 
result  of  a  change  in  angle-of-attack  caused  by 
a  change  in  w 

The  change  in  the  force  in  the  Z  direction  as  a 
result  of  the  time  rate  of  change  of 
angle-of-attack  caused  by  w 

The  change  in  the  force  in  the  Z  direction  due 
to  a  control  surface  deflection 


'z6tef 


The  change  in  the  force  in  the  Z  direction  due 
to  a  horizontal  tail  deflection 

The  change  in  the  force  in  the  Z  direction  due 
to  a  trailing-edge  flaperon  deflection 

C-star  longitudinal  response  criteria 


System  feedforward  matrix 


J 


Initial  distance  from  aircraft  to  target 

Coefficient  of  pitching  velocity  specifying  the 
amount  that  is  sensed  by  the  angle-of-at  tack 
sensor  as  a  result  of  its  position  relative  to 
the  CG  of  the  aircraft 

Decibel 

Distance  to  peak  vertical  gust  velocity  (feet) 

Discrete*- time  system  dynamics  matrix 

Force  in  the  Z  direction 

Forward-loop  scheduled  gain 

Attack  angle  feedback  scheduled  gain 

Gravitational  acceleration 

Vertical  acceleration 

D iscre te-t ime  system  control  matrix 

g  error  signal 

Vertical  displacement  of  the  aircraft  caused  by 
an  external,  disturbance 

Altitude,  in  feet 

Mass  moment  of  inertia  about  the  Y  axis  of  the 
aircraft 

System  feedback  matrix 

Horizontal  tail  command  forward-loop  gain 

Tr ailing-edge  flaperon  forward-loop  gain 

Horizontal  tail  forward-loop  gain 

C*  feedback  gain  of  normal  acceleration 
measured  at  the  CG  of  the  aircraft 

Trailing-edge  flaperon  command  forward-locp 
gain 

C*  feedback  gain  of  pitch  rate  based  upon 
crossover  velocity 

4  C *  feedback  gain  of  pitch  acceleration  based 
upon  the  distance  from  the  aircraft  CG  to  the 


pilot  station 

Attack  angle  feedback  gain  to  the  HT  channel 

Attack  angle  feedback  gain  to  the  TEF  channel 

Pitch  rate  feedback  gain  to  the  HT  channel 

Pitch  rate  feedback  gain  to  the  TEF  channel 

Pitch  rate  command  feedback  gain  to  the  HT 
channel 

Normal  acceleration  feedback  gain  to  the  TEF 
channel 

Lift  force  t 

The  distance,  in  feet,  from  the  wing  MAC  to  the 
horizontal  tail  MAC 

Mass  of  the  aircraft,  in  slugs 

System  disturbance  feedforward  matrix 

System  disturbance  input  matrix 

Static  pressure  at  Sea  Level 

Static  pressure  at  altitude 

The  change  in  pitching  velocity  resulting  from 
some  disturbance 

Dynamic  pressure 

Total  angular  pitching  velocity 

Equilibrium  angular  pitching  velocity 

Impact  pressure 

Pitch  rate  gyro  output 

TEF  washout  eigenvalue 

Laplace  variable 

Planform  area  of  the  wing,  in  sq  ft 

Longitudinal  target  tracking  error 

The  change  in  linear  forward  velocity  resulting 
from  some  disturbance  (perturbation  forward 
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H  1 

e 

I 

[ 

velocity) 

* 

i 

1 

1 

•u 

- 

Non-dimensional  perturbation  forward  velocity 

U 

- 

Total  linear  forward  velocity 

i 

Uo 

- 

Equilibrium  forward  velocity 

V 

- 

Vertical  gust  velocity  (ft/sec) 

V 

- 

Aircraft  velocity,  in  knots 

1 

vm 

- 

Vertical  gust  peak  magnitude  (ft/sec) 

| 

vs 

- 

Velocity  of  sound,  in  knots 

XT 

- 

Aircraft  velocity  vector 

o 

Eh 

>1 

- 

Equilibrium  aircraft  velocity  vector 

w 

The  change  in  linear  vertical  velocity 
resulting  from  some  disturbance  (perturbation 
vertical  velocity) 

i 

h 

V 

- 

Total  linear  vertical  velocity, 

positive  downward 

def ined 

H 

tj 

Wo 

- 

Equilibrium  vertical  velocity,  defined 
downward 

positive 

wg 

- 

Vertical  gust  velocity,  defined 

downward 

positive 

li 

"s 

- 

Vertical  velocity  of  the  angle¬ 

sensor,  defined  positive  downward 

of-attack 

Wa/A 

- 

Vertical  velocity  of  the  aircraft  relative  to 
the  surrounding  air 

\/E 

- 

Vertical  velocity  of  the  aircraft  relative  to 
the  Earth 

WA/E 

- 

Vertical  velocity  of  the  surrounding  air 

relative  to  the  Earth 

X 

- 

Command  prefilter  eigenvalue 

i 

X,  Y,  Z 

- 

Orthogonal  disturbed  aircraft  axes 

Xo»Yq»Zo 

- 

Orthogonal  equilibrium  aircraft  axes 

N 

Xe»Ye»Ze 

- 

Orthogonal  Earth  reference  axes 

4 

z 

- 

HT  P+I  cascade  compensator  zero 
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L 

■  II 

.  - 

-  •  .i  —y-.. —  ^  .  .  i 

I 

“ 

I"* 

4 

r 

% 

1 

1  -  cos 

'a 

a 


cto 

°6 

Om 

Os 

at 

aa/A 


Ve 

aA/E 


P 

6 

&e 


6f 


Discrete  vertical  gust  formed  by  displacing  an 
inverted  cosine  wave  a  unit  distance  above  the 
time  axis 

Perturbation  attack  angle  (~w/U) 

Angle-of-attack,  defined  as  the  angle  between 
the  velocity  vector  of  the  aircraft  and  the 
chord  of  the  wing 

Angle  between  the  aircraft  X  axis  and  the  chord 
of  the  wing 

Effective  attack  angle  generated  as  a  result  of 
the  vertical  velocity  of  the  surrounding  air 
relative  to  the  Earth,  defined  positive  for  an 
upward  gust  (-aA//E) 

True  perturbation  attack  angle  ('a  +  Og) 

Attack  angle  sensor  output 

Trim  angle-of-attack  (at  equilibrium) 

Effective  attack  angle  generated  as  a  result  of 
the  vertical  velocity  of  the  aircraft  relative 
to  the  surrounding  air 

Effective  attack  angle  generated  as  a  result  of 
the  vertical  velocity  of  the  aircraft  relative 
to  the  Earth 

Effective  attack  angle  generated  as  a  result  of 
the  vertical  velocity  of  the  surrounding  air 
relative  to  the  Earth 

TEF  servo  eigenvalue 

Control  surface  angular  deflection 

Elevator  deflection  (used  in te rchangably  with 

&h  ) 

Flaperon  deflection  (used  interchangably  with 
6tef  > 

Horizontal  tail  deflection 


6htL  -  L/H  horizontal  tail  deflection 
6htp  “  R/H  horizontal  tail  deflection 

6-tef  -  Trailing-edge  flaperon  deflection 
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6tefL 

6tefR 

Y 


Yo 

e 


L/H  trailing-edge  flaperon  deflection 
R/H  trailing-edge  flaperon  deflection 
Air  density  ratio 
Flight  path  angle 
HT  servo  eigenvalue 

Initial  (equilibrium)  flight  path  angle 
Perturbation  pitch  angle 

Initial  (equilibrium)  pitch  angle  (angle 
between  the  horizontal  and  the  aircraft 
axis)  , 


Appendix  D 

Derivation  of  Stability  Derivatives 


To  make  use  of  the  equation?  of  motion  presented  in 
Chapter  II,  it  ir.  necessary  to  have  certain  aerodynamic  data 
available  on  the  aircraft  to  be  modeled.  Generally  speakinq, 
stability  derivatives  are  usually  not  available,  and  they 
must  be  determined  from  lilt,  drag,  and  moment  data.  In  the 
next  several  pages,  non-dimensional  stability  derivatives  lor 

t 

six  flight  conditions  (KC s)  are  presented  as  derived  from 
sources  of  corrected  wind  tunnel  data  (Deferences 
10,19,25,26).  The  six  flight  conditions  are: 


F.nch 

A 1 1 i t  ude 

(1.6 

Den  I.evel 

(SL.) 

0.6 

30,000 

ft 

0.8 

SI. 

0.8 

30 ,000 

ft 

0.9 

SI. 

0.9 

30 ,000 

ft 

Since  the  YF-.16  is  unstable  only  in  sub-sonic  flight  (Kef  9), 
no  supersonic  flight  conditions  are  included.  In  addition, 
the  aircraft  is  most  susceptible  to  w i nd  gust  disturbances 
over  the  range  of  flight  conditions  chosen  (Kef  8). 

The  procedure  used  for  drtermininq  the  derivatives  and 
coefficients  i  r-  illustrated  for  the  flight  condition  used  in 


the  analyses  in  Chapters  III  and  IV  (0.8  Mach  at  EL).  Data 
fcr  the  remaining  FCs  was  determined  using  the  same  methods. 
Lata  for  all  FCs  is  summarized  in  Table  6. 

The  first  step  in  the  process  is  to  determine  how  much 
lift  is  necessary  to  maintain  the  aircraft  in  equilibrium 
flight.  This  occurs  when  the  force  generated  by  lift  is 
exactly  equal  and  opposite  to  the  force  of  gravity.  Lift  is 
given  by 


L  -  CLSq 


where  S  is  the  planform  area  of  the  wing,  q  is  dynamic 
pressure,  and  is  the  coefficient  of  lift.  Setting  lift 

equal  to  the  force  of  gravity  yields 


CLSq  =  mg 


where  m  is  the  mass  of  the  aircraft,  and  g  is  gravitational 
acceleration,  defined  positive  downward.  Solving  this 
equation  for  results  in  an  expression  that  determines  the 
necessary  coefficient  of  lift  to  maintain  equilibrium  flight 
at  any  given  altitude  or  airspeed.  This  expression  is 


UL  =  s? 

In  this  analysis,  the  GV.  (grcss  weight)  of  the  aircraft 
is  assum.ed  tc  be  16,519  lbs  at  EL.  As  g  =  32.81  at  EL,  the 
mass  of  the  aircraft  is  determined  tc  be 


m  =  —  =  16,519/32.18  =  513*4  slugs 

D 
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For  altitudes  other  than  Sea  Level,  the  weight  (mg)  will  vary 
with  g.  A  relation  for  gravitational  acceleration  as  a 
function  of  altitude  is  (Ref  14), 


g  =  32.177  -  (O.OOOO9406) H  (49) 

where  H  is  the  particular  altitude,  in  feet. 

The  dynamic  pressure,  q,  is  a  function  of  both  altitude 
and  airspeed.  It  can  be  determined  from  the  expression 


q  = 


(50) 


where  Y  is  the  air  density  ratio  at  the  altitude  of  interest, 
and  V  is  the  velocity  of  the  aircraft,  in  knots.  The  air 
density  ratio,  defined  as  the  ratio  of  air  density  at  a 

specified  altitude  to  air  density  at  Sea  Level,  is  a  function 

/ 

of  altitude  only,  and  its  value  may  be  determined  from  a 
standard  altitude  table.  To  express  q  as  a  function  of  Mach 
Number,  recall  that  Mach  Number  is  defined  as 


Mach  Number  =  M  =  —  (51) 

vs 

where  Vs  is  the  velocity  of  sound,  in  knots.  Combining  Eqns 
50  and  51,  dynamic  pressure  can  now  be  determined  from  the 
final  expression 


q  = 


(52) 


For  the  flight  condition  0.8  Mach  at  SL,  q=  949.9,  since  Y=  1 
and  Vs  =  661.7  knots.  The  planform  wing  area  of  the  YF-16  is 
S  =  280  sq  ft  (Ref  10).  The  required  lift  coefficient  can  now 
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be  readily  determined  as 


c  _  mg  _  (32.18)(513.4) 

L  Sq  (280)(949.9)  ' 

Given  the  value  of  C^,  the  attack  angle  at  which  the 
aircraft  will  trim  can  be  found.  This  establishes  the 
equilibrium  flight  condition.  Typical  corrected  wind  tunnel 
data  for  this  determination  is  shown  in  Figure  42  (Ref  26) . 
The  value  of  this  trim  angle-of-at  tack ,  at=  0.81  degrees,  can 

be  read  directly.  Actually,  tfye  CL  vs  a  curve  of  Figure  42  is 

a  family  of  curves,  and  the  curve  chosen  depends  upon  the 
amount  of  horizontal  tail  deflection  required  to  maintain 
equilibrium  flight.  Due  to  the  aircraft  velocity  at  this  FC, 
however,  the  horizontal  tail  deflection  is  almost  negligible, 
and  therefore  the  =  0  curve  is  appropriate.  With  an 

established  equilibrium  condition,  the  required  stability 
derivatives  and  coefficients  may  now  be  determined. 


Given  airspeed  and  at »  this  derivative  can  be  determined 
directly  from  available  data.  A  typical  vs  Mach  curve  is 
shown  in  Figure  43  (Ref  19).  For  the  FC  of  interest, 

%  =  -4.05  (54) 


b.  C. 


Data  for  the  determination  of  this  derivative  is  also 
available  plotted  as  a  function  of  Mach  Number  for  a  given 
at*  A  typical  curve  is  given  in  Figure  44  (Ref  19).  For  the 
FC  of  interest, 

Cma  =  _1,45  (55) 

c. 

/ 

This  derivative  may  sometimes  be  approximated  by  (Ref  1) 

%  =  (cAt)cmq  (56) 

where  It  is  the  distance  from  the  wing  MAC  to  the  horizontal 
tail  MAC.  For  the  YF-16,  It  =  15.66  ft  (Ref  10),  and 
calculated  in  this  way,  = -2.83.  However,  this  expression 
is  an  approximation,  and  is  generally  valid  only  for 
conventional  aircraft.  Due  to  the  placement  of  the  horizontal 
tail  relative  to  the  wing  of  the  YF-16,  this  method  of 
determining  may  not  be  as  accurate  as  desired.  It  is 

assumed  that  a  more  accurate  estimate  can  be  made  using  the 
relation 

%  *  "Clq  (57) 

Plots  of  Cjq  vs  Mach  Number  are  provided  in  Reference  19,  and 
for  the  FC  of  interest,  Ci^=3.86;  therefore, 

-  -3.86  (  58) 
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This  derivative  may  sometimes  be  approximated  by  (Ref  1) 


The  relationship  is  the  same  as  given  in  Eqn  56.  Calculated 


in  this  manner,  C 


airframe,  it  is  assumed  that  a  more  accurate  estimate  of  this 


derivative  can  be  made  using  the  relation 


vs  Mach  Number  (Ref  19),  the  value  at  the 


FC  of  interest  is  determined  to  be  C- 


In  this  case,  the  values  obtained  using  both  approximations 
are  fairly  close. 


This  term  is  simply  a  gravity  term,  and  is  equal  to  -Cj, 
For  the  FC  of  interest. 


The  only  unknown  in  this  quantity  is  the  forward  velocity 
of  the  aircraft,  U,  in  ft/sec.  Knowing. Mach  Number,  and  given 
that  the  velocity  of  sound  at  SL  is  661.7  knots,  the  velocity 
of  the  aircraft  in  knots  can  be  determined  from  Eqn  51.  As  1 
knot  is  equal  to  1.6878  ft/ sec,  the  velocity  of  the  aircraft, 
in  ft/sec,  is  determined  to  be 

U  =  (i.6878)(0.8)(66l.7)  =  893  (63) 


and  the  value  of  the  coefficient  can  now  be  found: 


_  (513.4)(893)  _ 
Sq  (280)(949.9) 


(64) 


Given  the  Mean  Aerodynamic  Chord  length,  c 
(see  Appendix  E),  the  value  of  this  coefficient 


=  10.937  feet 
can  be  found: 


_c_  _  10.937 
2U  2(893) 


0.00612 


(65) 


h. 

Sqc 

The  unknown  in  this  quantity  is  ly  ,  the  mass  moment  of 
inertia  about  the  Y  axis.  From  Reference  19,  Iy  = 39,199; 


therefore , 


upon  the  distance  from* the 
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'•vwAciaviun  oasea 
aircraft  CG  to  the 


39199 _ 

(280)(949.9)(10.937) 


0.0135 


(66) 


i. 

This  derivative  is  simply  the  change  in  CL  as  a  result  of 
a  change  in  a-  Therefore,  its  value  can  be  determined  by 

taking  the  slope  of  the  curve  in  Figure  42  at  the  operating 

/ 

point.  This  slope  is  approximately  0.0758  per  degree.  Since 
the  stability  derivatives  are  non-dimensional,  this  figure 
must  be  multiplied  by  57.296  degrees/radian.  This 
multiplication  results  in 

Clq  =  (-0758  X  57.296)  =  4.34  (  67) 


j  • 

This  derivative  in  the  change  in  Cm  (moment  coefficient) 
as  a  result  of  a  change  in  a.  Unfortunately,  a  vs  a  curve 
is  not  generally  available.  However,  Cma  can  be  broken  down 
into  two  terms  as  follows: 


Sx 


_  >Cb  _  JCjj  ;cL  _  iCn 


(6  8) 


«ia  JCj,  ia  ICl 
There  are  Cm  vs  CL  curves  available  (Hef  26),  and  for  the  FC 
of  interest,  the  slope  of  the  curve  is  approximately  0.0331. 
Since  both  Cjq  and  C*,  are  non-dimensional,  there  are  no 
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■me  c nange 
from  some 


in  linear  forward  velocity 
disturbance  (perturbation 


resulting 

forward 


•14? 


( 


conversion  factors  to  be  concerned  with.  The  product  of  this 
derivative  and  Cj^  is 

cma  =  (0.033l)(4.34)  =  0.144  (69) 

k.  cd 

This  term  is  the  coefficient  of  drag,  and  is  equal  and 

/' 

opposite  to  the  total  thrust  in  equilibrium.  Its  value  can  be 
read  directly  from  a  drag  polar  (Ref  26),  which  is  a  family 
of  cD  vs  Cj,  curves  that  vary  with  Mach  Number.  For  the  FC  of 
interest, 


CD  =  0.024 


(70) 


This  derivative  can  be  determined  from  the  relationship 
(Ref  1)  s 

Gza  =  -Cd  -  (71) 

Since  Cq  and  Ci^  have  already  been  determined,  the  expression 
is  useful,  and  this  derivative  is 

Cza  =  -0.024  -  4.34  =  -4.364  (  72) 
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This  derivative  can  be  determined  from  the  relationship 


I  Ref  1) : 


%=CL-^  (73) 

As  Cd  vs  a  curves  are  not  standard  aerodynamic  data,  this 
expression  must  be  broken  down  into 


r  -  icDacL/_  jcd 

Cxa  ~  CL  Jci ~  CL_^CLa  (74) 

The  only  unknown  term  is  the  change  in  CD  for  a  change  in  CL. 
This  information  may  be  obtained  from  a  drag  polar  (Ref  26), 
by  finding  the  slope  of  the  curve  at  the  operating  condition. 
This  slope  is  determined  to  be  approximately  0.0025,  and  is 
dimensionless.  Substitution  into  Eqn  74  yields: 

CX(X  =  0.0621  -  (.0025)(4.34)  =  0.0513  (  75) 


n.  C„ 


This  derivative  can  be  determined  from  the  relationship 


(Ref  1): 


CZU  =  -2cl  -  VdiT  (76) 

Since  aerodynamic  data  expressing  a  change  in  for  a  change 
in  V  is  not  generally  available,  one  method  that  can  be  used 
to  determine  the  derivative  is  to  collect  sufficient  data 
points  to  make  a  plot  of  vs  V.  From  such  a  plot,  the  slope 
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6htR 

&tef 


R/H  horizontal  tail  deflection 
Trailing-edge  flaperon  deflection 
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of  the  curve  at  the  operating  point  yields  the  required 
derivative.  The  available  data  is  in  the  form  of  vs  a 
curves  as  a  function  of  Mach  Number  (Ref  25).  Equation  76, 
therefore,  can  not  be  used  directly.  However,  since  Mach 
Number  is  defined  as  vehicle  velocity,  V,  divided  by  the 
velocity  of  sound,  Vs  ,  the  following  equivalence  can  be 
shown: 


dM 


(v/v^) — —  =  (v/v  ) — ^ _ 

^/  s)  d(v/vs)  y/v^  (lAs)du 


=  V 


dC 


dV 


(77) 


Therefore,  Eqn  76  can  be  written  as 

„  dCL 

Cz,  =  -2CL 


(78) 


which  can  be  used  directly  with  the  available  data. 

Values  of  as  a  function  of  both  Mach  Number  and 

angle-of-attack  are  plotted  in  Figure  45.  Since,  for  any 
flight  condition,  the  values  of  Mach  Number,  Cj, ,  and  attack 
angle  are  all  known,  any  two  of  the  three  will  locate  the 
operating  point.  The  slope  of  the  curve  at  that  operating 
point  is 


dCT  .075 

Slope  =  — -  =  — -  =  .0625 
dM  1.2 

Substituting  this  value  into  Eqn  78  yields: 


(79) 


CZu  =  -2(.062i)  -  ,8( .0625)  = 


(80) 


o. 
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- 


-  .174 


Figure  1*5.  Plot  of  as  a  function  of  Mach  Number  and  angle -of -attack. 


This  derivative  can  be  determined  from  the  relationship 
(Ref  1) : 

°xu  =  -2CD  ~  ^  (81) 

As  in  the  previous  case  for  data  expressing  a  change  in 

drag  coefficient  for  a  change  in  forward  velocity  is  not 
generally  available.  Using  the  equivalence  given  in  Eqn  77, 
it  can  be  shewn  that  Eqn  81  is  equivalent  to 


Cx  =  -2Cn  - 
xu  D  dM 


(82) 


Values  of  CD  as  a  function  of  both  Mach  Number  a  l 
angle-of- at  tack  are  plotted  in  Figure  46.  The  slope  of  the 
curve  at  the  operating  point  provides  the  required 
derivative.  For  the  FC  being  considered, 

dCD  .027 


dM  1.2  ~  *0225 

Substituting  this  value  into  Eqn  82  yields: 


(83) 


"xu  =  -2( .024)  -  .8( .0225)  =  -.0660 


(84) 


P*  Cm&ht 

This  derivative  is  given  as  cm6ht=“*01  per  degree  in 
inference  10.  Multiplying  by  57.296  degrees/radian  to  put  it 
-d  inrns  ional  form, 

C»6ht  =  ( - . 01 ) ( 57 . 296 )  =  -.573  (85) 
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Cm6tef 

This  derivative  is  given  as  cm6tef  =  -.0015  per  degree  in 
Reference  10.  Multiplying  by  57.296  degrees/radian  to  put  it 
in  non-dimensional  form, 

Cm6tef  =  (-.0015  X  57.296)  =  -.0859  (  86) 

r*  Cx6ht 


Reference  26  provides  a  family  of  curves  that  describe 
the  changes  in  the  coefficient  of  drag  that  result  from 
various  magnitudes  of  horizontal  tail  deflection  for  various 
Mach  Numbers.  The  approximate  change  in  total  drag  resulting 
from  a  horizontal  tail  deflection  can  be  determined  using 
linear  estimation.  Assuming  that  an  average  horizontal  tail 
deflection  at  this  flight  condition  is  approximately  5 
degrees  or  less  for  a  perturbation  analysis,  two  data  points 
can  be  extracted  from  the  given  information: 


ACD  6ht 

0  0 

.00291  5 


Assuming  that  cxtht  changes  linearly  for  horizontal  tail 
deflections  between  0  and  5  degrees,  the  slope  of  the  line 
indicates  the  desired  derivative.  The  non-dimensional  form  of 
the  result  is 
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.  ^•D 

Cx6h+  =  -~(57.296)  =  -(.000581X57.296)  =  -.0333  (  8  7) 

nx  6ht 

s  *  Cx6tef 

Reference  25  provides  a  family  of  curves  that  describe 
the  changes  in  the  coefficient  of  drag  that  result  from 
various  magnitudes  of  trailing-edge  flaperon  deflection  for 
various  Mach  Numbers.  The  same  technique  as  in  the  previous 
case  for  Cx$h.j.  is  used  to  determine  this  derivative.  In 
non-dimensional  form, 

Q 

Cx6.  =  -  —(57.296)  =  -  (.000722) (57. 296)  =  -.0414  (8  8) 

tel  &tef 

Cz6ht 

From  Reference  10,  it  is  given  that  CL6ht=  0.0085/degree. 
If  it  is  assumed  that  —  -Cl^^  #  then,  in  non-dimensional 

form, 

CZ6ht-  -  (.0085X57.296)  =  -.487  (89) 

U*  Cz6tef 

»■ 
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From  Reference  10,  it  is  given  that  Cz6tef  =  0.011/degree. 
If  it  is  assumed  that  cz6tef  “  _CL6tef'  then,  in 
non-dimensional  form. 


czfttef  -  -  (.011X57.296)  =  -.630 


(90) 


These  stability  derivatives  and  coefficients,  and  those 
for  the  remaining  five  flight  conditions  of  interest,  are 
summarized  in  Table  6.  / 


( 
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TABIE  6 


Stability  Derivatives  and  Coefficients 


L 


J 


Mach 

S  0.6 

0.8 

0.9 

Altitude 

SL 

30K 

SL 

30K 

SL 

30  K 

at 

1.6 

4.5 

0.8 

2.4 

0.5 

1.6 

U  (fps) 

670 

597 

893 

796 

1005 

895 

CL 

.1104 

•  339 

.0621 

t 

.191 

.0491 

.151 

CD 

.026 

.036 

.024 

.028 

.026 

.027 

mU/  Sq 

2.300 

6.899 

1.725 

5-175 

1.533 

4.600 

Iy  /  Sqc 

.0240 

.0807 

.0135 

.0454 

.0107 

.0359 

c/2U 

.00816 

.00916 

.00612 

.00687 

.00544 

.00611 

cw 

-  .1104 

-.339 

-.0621 

-.191 

-.0491 

-.151 

Omq 

-3.60 

-3.65 

-4.05 

-4.12 

-4.80 

-4.95 

Ozq 

-3.25 

-3.27 

-3.86 

-3.92 

-4.76 

-4.91 

CLa 

4.2? 

4.41 

4.34 

4.55 

4.86 

4.92 

cma 

.183 

.135 

.144 

.0919 

.153 

.0974 

cma 

-.950 

-1.15 

-1.45 

-1.47 

-1.90 

-2.20 

cxa 

.0216 

-.113 

.0513 

-.0140 

.0406 

.0647 

Cza 

-4.30 

-4.45 

-4.36 

-4.58 

-4.89 

-4.95 

Cza 

-1.20 

-1.22 

-1.08 

-1.04 

-1.36 

-1.35 

cxu 

-.061 

-.0811 

-.0660 

-.0767 

-.112 

-.108 

Czu 

-.240 

-.722 

-.174 

-.477 

-.219 

-.525 

Cm6ht 

-.590 

-.607 

-.573 

-.602 

-.573 

-.619 

Cm6tef 

-.0544 

-.0607 

-.0859 

-.100 

-.118 

-.154 

Cx6ht 

-.0305 

-.0616 

-.0333 

-.0644 

-.0362 

-.0514 

Cx6tef 

-.0384 

-.0688 

-.0414 

-.0722 

-.0426 

-.0593 

EE 

-.536 

-.544 

-.487 

-.516 

-.498 

-.521 

HBS81 

-.688 

-.802 

-.630 

-.802 

-.544 

-.796 
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Appendix  C 


State  Space  Modeling  of  the 
Longitudinal  Equations  of  Motion 

When  analyzing  or  designing  a  system  with  state  variable 
feedback,  as  in  Chapter  IV,  it  is  most  straightforward  to 
model  the  system  in  matrix  format.  In  this  way,  each  state  is 
readily  identified  as  a  feedback  source,  and  the  closed-loop 
system  can  be  obtained  with  simple  matrix  manipulations.  This 
technique  becomes  especially  attractive  when  the  added 
complexity  of  servos  and  dynamic  compensation  are  considered. 
Modeling  the  equations  of  motion  of  the  aircraft  and  all 
peripheral  system  dynamics  individually  in  "state  space" 
format  permits  the  composite  system  to  be  obtained  by  matrix 
augmentation  (see  Appendix  G)  . 

The  thr ee-degree-of- freedom,  non-dimensional,  longitud¬ 
inal,  stability-axis  equations  of  motion  used  in  this  inves¬ 
tigation  are  presented  in  Chapter  II  as  Eqns  1,  2,  and  3. 
Assuming  0O=  0 ,  and  substituting  a  unique  upper  case  letter 
for  each  coefficient  of  each  state  or  input  variable,  the 
equations  can  be  written  in  a  more  convenient  form  as 
follows: 


1 

3 

1 

•3 

< 

c'a  - 

D  0  =  X  6 

(91) 

-E'u  +  F'a  - 

G'a  - 

H  0  =  Z  6 

(92) 

-I'd  -  J'a  + 

K0  - 

L  0  =  M  6 

(93) 

The  states  of  this  system  of  equations  are  perturbation  in 
forward  velocity,  'u  ,  perturbation  in  angle-of-a  ttack,  'a. 
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perturbation  in  pitch  angle,  6,  and  perturbation  in  pitch 
rate,  0  or  q.  The  input,  6,  is  representative  of  any  general 
input  (horizontal  tail  or  trail ing-edge  flaperons,  or  wind 


gust  disturbances . as  used  in  this  report).  Equations  91,  92, 
and  93  can  now  be  solved  for  the  first  derivative  of  each  of 


the  system  states: 

•u  =  (B/A)’u  +  (C/A)  'a  +  (d/a)  0  +  (X/a)  6  (94) 

=  (E/F)’u  +  (G/F)  ’a  +  (H/F ) q  +  (z/F)  6  (95) 

e  =  q  '  (96) 

4  =  (J/K)’a  +  (I/K)  'a  +  (L/K) q  +  (m/k)  6  (97) 


Because  Eqn  97  contains  an  "a"  term,  Eqn  95  must  be 
substituted  to  allow  the  derivative  term  on  the  left  of  each 
equation  to  be  expressed  in  terms  of  system  states  only. 
Performing  this  substitution,  the  new  set  of  equations  is 


•u  =  (b/a)  'u  +  (C/A)  ’a  +  (d/a)  0  +  (X/A)  6 
’a  =  (E/F)  *u  +  (G/F)  'a  +  (H/F)q  +  (Z/F)  6 
0  =  q 

q  =  (IE/KF)  'u  +  [(JF  +  IG)/KF]  'a 

+  [(HI  +  LF)/KF]q  +  [(IZ  +  MF)/KF]6 

These  equations  can  now  be  put  into  matrix  form 


V 

U1  U2  U3  o 

•u 

Uin 

•d 

A1  A2  0  A4 

•a 

^in 

• 

— 

+ 

0 

0  0  0  1 

0 

0 

• 

9 

- 1 

& 

o 

& 

<y 

_ i 

9 

172 


(98) 

(99) 
(100) 

(101) 


(102) 


where  the  terms  U1  ,  .  etc .  represent  the  state  and 
input  variable  coefficients  in  Eans  98  through  101.  Defined 
in  terms  of  the  original  parameters  in  Equations  1,  2,  and  3, 
the  matrix  entries  of  Eqn  102  are  given  in  Eqns  103  through 
114  as  follows: 


U1  = 


(cxu)Sq 

mU 


U2  = 


U3  = 


u. 


in  = 


(Cxa)Sq 

'mU 

(Cw)Sq 

mU 

(CxjJSq 

mU 


A1  = 


jzu 


(mU/Sq)  -  (c/2U)Cz- 


A2  = 


EflL 


Q1  = 


(mU/Sq)  -  (c/2u)Czi 

(mU/Sq)  +  (c/2u)G^ 
(mU/Sq)  -  (c/2u)Czi 

_ 2^6 _ 

(mU/Sq)  -  (c/2u)Gz^ 

{(c/2U)Cmg}CZu 

(ly/Sqc){mU/Sq  -  (c/2u)Czi] 


A4  = 


^in  ~ 


Q2 


{mU/Sq  -  (c/2U)Cz^}Cma  +  (c/2U)(Cm^)CZg 
(ly/Sqc)£mU/Sq  -  (c/2U)CZ(i} 


(103) 

(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

(110) 

(111) 

(112) 


{mU/Sq  +  (c/2u)CZq3(c/2U)cffi^  +  (mU/Sq  -  (c/2u)cZ(g(c/2u)c^ 
(ly/Sqc) {mU/Sq  -  (c/2U)Cz£[} 


(113) 
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(c/2U)(Cmfr)Cz6  +  QnU/Sg  -  (c/2u)Cz&}Cm6 
ln  (ly/Sqc){mU/Sq  -  (c/2u)CZ(i3 


(114) 


The  "Short  Period"  model  for  the  aircraft  system  is 
simply  a  reduced  form  of  Eqn  102.  It  is  obtained  as  a  result 
of  two  basic  assumptions: 

1.  Small  perturbation  in  pitch  have  a  negligible  effect 
on  aircraft  forward  velocity. 

2.  Control  surface  deflections  introduce  significantly 
small  changes  in  total  ^rag  so  as  to  cause  negligible 
changes  in  aircraft  forward  velocity. 

These  assumptions  permit  the  exclusion  of  the  Eqn  98  from  the 
system  state  matrix.  Also,  since  neither  Eqn  99  or  Eqn  101 
contain  pitch  angle  terms,  deletion  of  Eqn  100  does  not 
result  in  the  loss  of  any  information.  The  two-state  Short 
Period  state  space  model  may  be  expressed  as 


'a  A2  A4  'a  Ain 

;  =  + 

_4J  |_Q2  Q4j[q.  J  |_Qin 


(115) 


where  each  term  is  as  defined  in  Eqns  108,  109,  110,  112,  113 
and  1 14  . 


The  scheduled  gains  used  in  the  YF-16  CCV  longitudinal 
FCS  are  functions  of  impact  pressure  (qc),  the  ratio  of 
impact  pressure  to  static  pressure  at  altitude  (qc/Ps)r  or 
the  ratio  of  static  pressure  at  altitude  to  static  pressure 
at  sea  level  (Ps/P0).  To  permit  evaluation  of  these  gains. 
Figures  47,  48,  and  49  provide  values  of  qc  ,  qc/Ps •  and 
Ps/P0»  respectively,  as  functions  of  known  flight  condition 
parameters.  Graphs  and  equations  for  the  scheduled  gains  are 
given  in  Figures  50  through  55. 
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2.25 


Figure  48.  Ratio  of  impact  pressure  to  static  pressure  at  altitude  as  a  function  of 
Mach  Number  (Ref  10). 


ALTITUDE  (Thousands  of  feet) 


Figure  49.  Ratio  of  static  pressure  at  altitude  to  static  pressure  at 
Sea  Level  as  a  function  of  pressure  altitude  (Ref  10). 
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AD-A080  520  AIR  FORCE  INST  OF  TECH  VRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/6  1/3 
USING  VERTICAL  GUST  ALLEVIATION  TO  IMPROVE  THE  TARGET  TRACKING  —ETC 
AUG  79  E  R  MOLNER 

UNCLASSIFIED  AFIT/G6C/EE/79-5  M( 


Figure  50.  Scheduled  gain  F3  as  a  function  of  q.c  (Ref  10) 


-1 


F4  =  -i.l905(qc/Ps  )  +  1.1310  (  .53  <  Ic/Ps  £  1.79  ) 


Figure  51,  Scheduled  gain  F4  as  a  function  of  q.c/  P8  (Ref  10) 


Figure  53.  Scheduled  gain  6e/6f  as  a  function  of  qc/  (Ref  1°) 


165.0(^0/1*8  )  +  4°.°  (0.0^tc/l*s  <  0.2) 


(  0.2  £  qc/Ps  <  0.7  ) 


17.7l(qc/l*s  )  -  5.40  (  0.7  £  (ic/Ps  £  3*0) 


1845 /afpe 


1845/afps  =  -0.26(Ps/Po  )  +  1.91 


■2 . 00(  qc/Ps  )  +  2.35  (  0.0  <  q.c/Pg  <  0.165  ) 


1.93(  qc/Ps  )  +  2.34  (  0.165  <  qc/Ps  <  0.6  ) 


438(  qc/Ps  )  +  1 .44  (  0.6  <  qc/Ps  <  1.4  ) 


133(  qc/Ps  )  +  1.02  (  1.4  <  qc/Ps  <  6.0  ) 


Figure  55.  Scheduled  gain  1845 /afps  818  a  function  of  Ps/P0  and 
determination  of  1  /  M  as  a  function  of  qc  /  Ps  (Ref  10) 
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Appendix  E 

Determination  of  Sensor  Locations 

As  the  position  of  the  normal  accelerometers  and  attack 
angle  sensors  must  be  determined  relative  to  the  center  of 
gravity  (CG) ,  the  location  of  the  CG  must  first  be  known.  To 
determine  its  location,  certain  structural  data  must  be 
available  to  describe  the  wing  of  the  aircraft.  The  following 
information  is  known  (Ref  10): 


1) 

Most 

forward 

point  of  wing:  FS  218.80 

2) 

Wing 

sweep: 

40  degrees 

3) 

Wing 

tip  location:  SS  173.90 

4) 

Wing 

tip  chord:  42.97  inches 

5) 

Taper 

ratio: 

0.2275 

6) 

Aspect  ratio 

:  3.00 

7) 

Wing 

area: 

280  square  feet 

8) 

Most 

forward 

CG:  FS  319.70 

9) 

Most 

aft  CG: 

FS  326.20 

10) 

Wing 

vertical  centerline:  WL  91.0 

Since  the  wing  tip  is  173.90  inches  from  the  centerline, 
the  wing  span  is  twice  that  measurement,  or  347.80  inches. 
Referring  to  Figure  56,  the  distances  dj  and  d£  are 
determined  to  be 


*  tan  (50*) 


=  145.92  inches 


d2  =  d^  +  42.97  =  188.89  inches 


(116) 
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For  reassurance,  some  of  the  given  data  may  now  be  checked. 
The  wing  area  is  given  by 

.  A  =  (145.92)(173.9)  +  (2)(42.97)(173.9)  =  40320.5  in2 

=  280  ft2  (117) 

The  wing  taper  ratio  is  given  by 


T 


42.97 

188.89 


O.2275 


(118) 


Lastly,  the  wing  aspect  ratio  is  given  by 


/P  _  (347.80)2 

40320.5 


3.00 


(119) 


From  the  above  checks,  it  is  fairly  certain  that  the  given 
dimensional  data  is  consistant. 

As  a  final  step  before  locating  the  CG,  the  wing  Mean 
Aerodynamic  Chord  (MAC)  must  first  be  known.  This  can  be 
determined  from  the  known  dimensional  data  (Hef  1): 


c 


(2/3)  188.89  +  42.97 


(188.89X42.97) 

188.89+42.97 


=  131.24  inches  (120) 

As  the  CG  has  been  defined  to  be  at  35%  MAC  for  the  analysis 
in  this  thesis,  the  CG  location  is  determined  to  be 

CG  =  FS  218.8  +  (188.89  -  131.24  +  ( .35)(131.24)) 


FS  322.4  ,  WL  91.0 


(121) 


This  value  falls  within  the  limits  of  the  most  forward  and 
most  aft  CC-  locations  previously  given,  as  it  should. 

The  location  of  the  angle-of-at tack  sensors  is  known  to 
be  FS  55,  WL  88.51  (Ref  19).  In  the  longitudinal  plane,  the 
distance  from  these  sensors  to  the  CG  can  now  be  readily 

fou nd: 

da  =  J(322.4  -  55)2  +  (91.0  -  88. 51)2'  =  26?. 4l  inches 

=  22.28  feet  (122) 

The  location  of  the  normal  accelerometers  is  known  to  be  FS 
167,  WL  90.0  (Ref  19).  In  the  longitudinal  plane,  the 

distance  from  these  accelerometers  to  the  CG  is  determined  to 
be 

dAj|  =  J(322.4  -  167. 0)2  +  (91.0  -  90. 0)2'  =  155.4  inches 

=  12.95  feet  (123) 


( 
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Appendix  F 

Rejected  System  Configurations 


As  is  mentioned  in  Chapter  IV,  three  basic  system 
configurations  were  investigated  in  an  attempt  to  determine 
which  is  most  effective  in  meeting  the  defined  design 
criteria.  Structured  as  G-command,  C*  command,  and  pitch  rate 
command  systems,  these  configurations  make  use  of  both  the 

horizontal  tail  and  the  tr ai li ng-e dge  flaperons  for  control. 

/ 

From  a  comparison  of  the  transfer  functions  for  these 
systems,  the  pitch  rate  command  system  described  in  Chapter 
IV  was  determined  to  be  potentially  the  most  capable  of 
controlling  aircraft  pitch  responses  to  vertical  gusting. 
This  decision  was  based  primarily  on  the  form  of  the  lowest 
order  term  in  the  numerator  and  the  denominator  of  the 
transfer  function  relating  pitch  rate  output  to  a  disturbance 
input.  In  the  case  of  the  G-command  system,  shown  in  Figure 
57,  these  two  terms  are  almost  identical  -  regardless  of 
system  gain  levels  -  indicating  absolutely  no  freedom  to 
reduce  the  overall  magnitude  of  the  numerator  relative  to  the 
denominator.  In  the  case  of  the  C*  command  system,  shown  in 
Figure  58,  the  addition  of  a  pitch  rate  term  to  the  normal 
acceleration  feedback  used  to  provide  the  controlled  system 
variable  results  in  improved,  but  still  insufficient, 
flexibility.  Only  a  small  reduction  in  the  relative  magnitude 
of  the  numerator  is  possible  compared  to  that  available  with 
the  selected  pitch  rate  command  system  shown  in  Figure  26. 

For  reference,  the  transfer  functions  derived  for  the 
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Figure  58.  C*  command  flight  control  system  configuration. 


i 


G-command  systeir.  are  given  in  Eqn6  124,  125,  and  126;  and 
those  derived  for  the  C*  command  system  are  given  in  Eons 
127,  128,  and  129.  As  both  of  these  systems  are  configured  to 
control  some  degree  of  normal  acceleration,  it  is  probable 
that  each  is  capable  of  providing  reduced  normal  acceleration 
response  to  vertical  gusting  when  compared  to  the  performance 
of  the  selected  pitch  rate  control  system. 
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q/p  Numerator  =  [  122.64Kf  ]  s3  +  [  398 . 74Kf  +  122.64-yKf  +  846.34KhKc  ]  s3 

+  [  398.74vKf  +  846.34pKhKc  +  2201.7KhKe  -  5^66.6KhKfK1  +  2649.5KhKfKc  -  2649.5KhKfKcK6 ]  s 
+  [22O1.70KhKc  +  .33393KhKfKQ  -  .33393KhKfKcK6  +  5466.6KhKfKcK2]  (126) 


[  22O1.70KhKc  -  .33393KhKfKcK6  +  5k66 .6KhKfKcKz  +  . 33393KhKfKc  ]  (129) 


Determination  of  the  System  Equation 


The  system  configuration  chosen  for  the  design  presented 
in  Chapter  IV  is  shown  again  in  Figure  59  for  easy  reference. 
Beginning  with  the  input  to  the  aircraft  equations  ("A/C 
DYNAMICS"  block  in  the  figure),  the  transfer  function  of  the 
horizontal  tail  servo  and  the  preceeding  gain  is 


&h_  =  -20Kh 

^  s  +  y 


(130) 


Cross-multiplying,  shuffling  terms,  and  removing  the  Laplace 
notation  results  in 

6h  =  -Y&h  "  20Kheh  (131) 

which  is  the  desired  differential  equation. 

The  transfer  function  of  the  trailing-edge  flaperon  servo 
and  the  preceeding  compensation  is 


il.  =  _ ~  20Kf  S  (132) 

ef  s2  +  (R  +  B) s  +  R3 

Again,  simple  manipulations  result  in 

6f  -  -(R  +  p)6f  -  Rp6f  -  20Kfef  (133) 

To  remove  the  input  derivative,  this  equation  is  integrated, 
resulting  in 

bf  =  -(R  +  P)6f  -  RB$6f  -  20Kfef  (134) 


Defining  ' x' ,  such  that 


x  =  6f  (135) 

and  substituting  Eqn  135  into  Egn  134,  the  desired 
differential  equation  is  obtained: 

6f  =  -(R  +  3)6f  -  Rpx  -  ZOKfej-  (136) 


Combining  the  differential  equations  given  in  Equations  131, 

/ 

135  and  136  into  a  single  matrix  expression,  a  three-state 
system  is  obtained  as  an  input  to  the  aircraft: 


(137) 


Using  matrix  augmentation,  Eqn  137  can  be  combined  with 
the  aircraft  Short  Period  matrix  equation  given  in  Eqn  38  to 
form  a  system  with  the  e-vector  as  the  input.  This  is  the 
inner  portion  of  the  system  shown  in  Figure  59  with  no 
feedback.  This  new  system  state  equation  is  given  as 


Additionally,  Eqn  137  is  combined  with  Eqn  39  to  produce  the 
new  system  measurable  output  equation 


(139) 

Backing  up  another  step,  such  that  the  r-vector  is  the 
system  input,  the  following  equations  result: 

(140) 

(141) 


(14  2) 

Letting  the  system  state  equation  (Eqn  138)  be  expressed 
as 

x  =  Ax  +  Be  +  NOg  (14  3) 

the  system  measurable  output  equation  (Eqn  139)  be  expressed 
as 

y;  =  Cx  +  M  ctg  (14  4) 

and  the  inner-stage  feedback  loop  equation  (Eqn  142)  be 
expressed  as 

e  =  r  -  K£  (14  5) 

Then,  combining  Eqns  144  and  145  results  in 

e  =  r  -  K(Cx  +  Motg)  =  r  -  KCx  -  KMOg  (14  6) 


eh  =  rh  "  Kl«s  " 

ef  =  rf  -  K2Os  -  K4qs  -  K6ANg 


These  equations,  written  in  matrix  form,  become 


«h  =  ^  _  K1  k3  0 
ef  rf  K2  K4  K6 
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Combining  Eqn  143  with  Eqn  146 ,  and  collecting  terms  provides 
the  desired  final  result: 

x  =  (A  -  BK  C)x  +  Br  +  (N  -  BKM)cig  (147) 

This  expression  allows  a  straightforward  calculation  of  the 
closed-loop  system  with  the  r-vector  as  the  input.  Performing 
the  matrix  multiplications  and  manipulations,  the  resultant 
state  and  output  equations  are  given  as  Equations  148  and 
149. 

Backing  up  the  next  step,  the  transfer  function  of  the 
proportional  plus  integral  compensation  is 


fh  _  Kc(s  +  z) 
uh  s 


(150) 


Cross-multiplication  and  removal  of  Laplace  notation  results 
in  the  following  differential  equation: 


rh  =  Kcuh  +  zKcuh  (151) 

To  eliminate  the  input  derivative,  this  equation  must  be 
integrated,  resulting  in 

*h  =  Kcuh  +  zKc^uh  (152) 

If  'y'  is  defined  such  that 


y  =  ^ 

then  Eqn  152  can  be  expressed  as 


(153) 
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+  ^cy 


(154) 


XV,  = 


as  desired. 

The  proportional  input  path  to  the  trailing-edge 

flaperons  contains  a  simple  gain,  and  the  resultant  transfer 
function  is  simply 

rf  =  KpUf  (155) 

Equations  148,  153,  154,  and  155  combine  to  produce  the 

/ 

open-loop  state  equation  with  the  u- vector  as  input.  This 
equation  is  given  in  Eqn  156.  Equation  157,  the  measurable 
output  equation,  is  changed  from  Eqn  149  only  by  the 
additional  state. 

Closing  the  final  feedback  loop, 

uh  =  P  "  k5<1s  (158) 

and 

Uf  =  p  (159) 

Equation  147,  defining  the  closed-loop  system,  may  be 
utilized  again  -  with  different  matrices,  of  course  -  to 
solve  for  the  system  with  scalar  input,  F.  The  final  system 
state  equation,  exclusive  of  the  command  prefilter  and  sensor 
dynamics  only,  is  given  in  Eqn  160.  The  measurable  output 
equation  is  unchanged,  and  is  given  as  Eqn  157.  The  command 
prefilter  need  not  be  modeled  in  the  matrix  equation,  as  it 
can  be  easily  incorporated  at  a  later  time  in  the  design 
process. 
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endix  H 


A££ 

Solution  of  the  System  Equation 
tor  Transfer  Functions 

The  closed-loop  system  equation  derived  in  Appendix  G  and 
given  as  Eqn  160  is  of  the  general  form 

x  =  Ax  +  Bp  +  Nag  (161) 

where  'A'  is  the  system  dynamics  matrix,  'b'  is  the  system 
control  matrix,  and  'N '  is  the  disturbance  input  matrix.  The 
characteristic  equation  (CE)  for  the  system  can  be  derived 
from  the  system  dynamics  matrix  using  the  relationship  (Ref 
6) 

CE  =  |  £(  s )  |  =  |  [si  -  A]|  =  det[  s  I  -  A]  (162) 

The  resulting  matrix,  Q^(s ),  is  given  in  Equation  163.  The 
determinant  of  C_(s)  can  be  determined  by  expanding  by  minors. 
Choosing  the  entries  in  column  four,  rows  four  and  five  as 
"pivots",  |q(b)|  can  be  expressed  as  the  difference  of  two 
two  5x5  terms  given  in  Eqn  164.  The  process  is  continued  in 
this  manner  until  no  further  reductions  are  possible.  The 
final  result  is  given  in  Eqn  165. 

To  determine  the  numerator  equation  for  the  transfer 
function  relating  pitch  rate  output  to  a  pilot  command  input, 
the  system  B  vector  is  substituted  for  the  second  column  in 
the  Q(s)  matrix,  as  per  Cramer's  Rule.  The  result  is  given  in 
Eqn  166.  Solving  for  the  determinant  of  £/(S)  will  produce 
the  desired  numerator  equation.  The  final  result  is  given  in 
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Eqn  167. 


The  numerator  equation  for  the  transfer  function  relating 
pitch  rate  output  to  a  disturbance  input  is  found  in  a 
similar  fashion.  Substitution  of  the  system  N  vector  for  the 
second  column  in  the  Q(s)  matrix  results  in  the  matrix  0"(s), 
given  in  Eqn  168.  Again,  solving  for  the  determinant  yields 
the  desired  numerator  equation,  given  in  Eqn  169. 
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CE  [-(•♦d)  ]  •*’  ♦  [  (ad  -  be)  -  (»  ♦  d)(P  ♦H»y)  +  RS«-Py 

♦  »Y  -  20  hKjX4  ♦  .5hKf»2  -  20hpXfXfc  -  20fX,*2  '  ZOfnKfK^  ♦  20(a  ♦  d)vKfH£  -  20^*^ 

-  20gKhKcK5  -  20*KhK1  -  20gXhX3  ♦  .^hK,  ]  •*  ♦  [  (ad  -  bc)(p  ♦  I  ♦  t)  -  (a  ♦  d)(H0  ♦  By  ♦  By) 

♦  Rfrl  '  ZOgdthKjKj  -  20(cf  -  alOKfKj,  ♦  ,5(cf  -  •*i)KfK2  -  20(cf  -  ahJpXjKg  -  20hYKfX4( 

/ 

♦  .^>1*^2  -  20hpY*fK6  ♦  20(df  -  th)X|X2  -  20fYKfX2  ♦  20(df  -  bhJnX^  -  UOOhuKhXf-K^Xj 

-  20fnYKfK6  ♦  400(gv  -  ♦  400(ay  -  fu)KhXfK6El  ♦  10{hu  -  gy^Xj^X,  -  20*0  ♦  R)XhXj 

-  20(ad  -  be )vKfK£  ♦  20(4  ♦  djvy*^  -  20(g  ♦  R^X^j  ♦  40QgyXhXfKcX<XJ  -  20g(@  ♦  H)KhKj 

♦  .Self  *  B)XhX,  +  20(d4  -  tg)KhKj  ♦  20(ag  -  ca^XjItj  ♦  ,5(c.  -  agfehXj  ♦  20(ag  -  ca)KhX3]aJ 

♦  [(id  -  bc)(Rfl  ♦  PY  ♦  «V)  *  (»  ♦  <OBPy  -  ZOgRBXhKcItj  -  20gR(iXhX3  -  20aR(KhXl  ♦  .SeRBR,^ 

♦  20(ag  -  MirthKeEj  -  20gi(KhKcKj  -  20gRzXhXcXj  ♦  400(gy  -  hu)«KhXfrcK6Kj  -  20(cf  -  t*h*t*k 

♦  Met  -  ah)yXfX2  -  20(ef  -  4h)pYXfX6  -  400(gf  -  •h)KhKfK4X,  -  400(gf  -  •h)pKhKfK6K1 

♦  20(df  -  bh)YKik2  *  *  •b)KhKfKcK3X2  .♦  400(gf  -  eh)XhKf*3K2  ♦  400(gf  -  ah  inX^fK^Kj 

♦  20(df  -  W>)nYKf*6  ♦  4<>0(gr  -  ahJnXnXfXfeX-j  ♦  1°<4*>  -  gOnXhX^Xj  *  400<dr  "  bh)“*h't««6*j 

♦  400(ah  -  cf)uKhKfKcK6K3  ♦  400(ah  -  cf )uKhKfX6IC3  ♦  10(ef  -  «b)ul(hXfX6X1  20(ad  -  be^^ 

♦  20(d.  -  bg)(p  ♦  R)XhK1  -  400(d.  -  bg)vXhXfX6X,  ♦  20(ag  -  c.)(R  ♦  B)XhXcXi 

-  400(ag  -  fr0»XhXfXo,t6l!j  ♦  2°<»«  -  «)(B  ♦  B)XhX3  -  400(44  -  c.)vXhX,*6X3 

♦  .5(C4  -  ag)(p  ♦  R)XhXj  -  10(c#  -  tgJyXhXjXjXj  ]  a»  ♦  [  (ad  -  be)RPY  ♦  20(da  -  bg)R@KhX, 

♦  20(44  -  cajRflXfcXeXj  ♦  20(44  *  ea)R0XhX3  ♦  .J(c4  -  agjRpXhXj  -  20gRaSXhXcXJ 

♦  20(44  *  ca).*,^  ♦  20(44  -  ca)RriChXcXi  ♦  400(gf  -  ah)rthXfXcX3X2  ♦  400(4 *  -  ahJntfhXjXeXjXj 

♦  400(ah  -  efJuriChXfXcXfcXj  ♦  400(ca  -  ag^aX^^Xj  -  400gy*XhX1KcX6XJ  ]  a 

♦  [  20(ag  -  ca)*ReXhXcX}  ] 


(165) 
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20(ag  -  ce)RpKhKc]  s  +  [  20(ag  -  ce)zRpKhKc  ]  (167) 


20(hi  -  fm)nvKfK6  +  400(gf  -  ehJxKbKfKgKj  +  400(eh  -  gf )nKhKfK6K1  ]  ss 
[20(gi  -  mORP^i  +  20(ce  -  ag)R3KhK1  +  (cl  -  am)RPY]s  +  [o]  (l69) 


Appendix  I 


Design  of  a  First-Order  Observer 


To  be  useful  in  controlling  aircraft  motion,  the 
electrical  outputs  of  a  flight  computer  must  be  converted 
into  control  surface  deflections.  In  the  YF-16,  these 
electrical  signals  drive  electromechanical  command  servos, 
which  in  turn  drive  hydraulic  power  actuators.  As  previously 

discussed,  the  command  servo  used  with  both  the  horizontal 

/ 

tail  and  trailing-edge  flaperons  has  two  outputs  available. 
One  output,  obviously,  is  mechanical  in  nature,  and  is  the 
controlling  input  to  a  power  actuator.  The  remaining  output 
is  a  variable  resistance,  the  amount  of  resistance  being 
dependent  upon  the  servo  position.  Such  an  output  can  be  used 
to  provide  a  voltage  or  current  that  is  proportional  to  servo 
position,  permitting  servo  position  to  be  input  to  the  flight 
computer. 

If  it  is  desired  to  include  the  second-order  servo 
dynamics  in  the  derivation  of  a  full  state  feedback  FCS,  an 
additional  output  must  be  derived  to  allow  feedback  signals 
representing  each  of  the  two  states  (Ref  6).  The  electrical 
input  to  the  servo  and  the  electrical  output  described  above 
can  be  used  to  derive  this  required  output  through  the  use  of 
a  simple  first-order  observer  that  can  be  easily  implemented 
within  the  flight  computer. 

If  the  dynamics  of  the  command  servo  are  expressed  as 


(Servo)  sa  +  ?2.8s  +  2704 


(170) 


I 


then  the  system  differential  equation  is  given  by 


y  =  -72.8  y  -  2704 y  +  2704  z 


(171) 


Defining  the  two  states  as  xj  and  x2 ,  such  thatxj=y  and 
x2  =  y»  then 


X1 

x2 


0  1 
-2704  -72.8 


X1 

x2 


0 

2704 


and 


■  0  •] 


X1 

x2 


which  can  also  be  written  as 


(172) 


(173) 


x  =  Ax  +  Bz 
y  =  CTx 


(174) 

(175) 


A  first-order  observer  state,  w ,  is  now  defined  such  that 
(Refs  12,  15) 


w  =  Tx 

(176) 

w  =  Dw  + 

ECtx 

+  TBz 

(17  7) 

TA  -  DT  = 

ect 

(178) 

where  D  and  E  are  scalars 

(in 

this  first-order 

case) 

representing  the  observer  time  constant  and  a  system  output 
matrix  multiplier,  respectively,  and  the  vector  T  is  a  system 
state  matrix  multiplier.  As  the  definable  observer  time 
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constant  will  determine  how  readily  the  observer  output  is 
able  to  settle  to  the  actual  value  of  the  observed  state,  a 
fairly  fast  time  constant  is  desirable.  This  time  constant  is 
chosen  to  allow  a  response  of  e  ~100t^  and  thus 

D  =  -100  (179) 

Also,  for  simplicity,  "E"  may  assume  a  value  of  unity. 

All  values  in  Eqn  177  are  now  known,  with  the  exception 
of  the  vector  T.  if  this  vector  is  expressed  as 


i  =  t2] 

then  Eqn  178  can  be  written  as 


(180) 


[\  *2] 


0  1 
-2704  -72.8 


+  100  tg]  =  [l  0] 


(181) 


The  values  of  tj  and  t2  can  be  found  by  expressing  this 
equation  in  the  form  of  two  equations  in  two  unknowns: 


100  tj  -  2704  12  =  1 


(182) 


+  27.2 12  =  0 


(183) 


Solving , 


t1  =  -.000184366 

t2  =  .00501475  (184) 

Equation  177  can  now  be  written  as 


w 


-100 w  +  y  -  0.49853z 


(185) 


An  output  vector  can  now  be  formed  from  Eqns  175  and  176, 
such  that 


or 


(186) 


-  — 

w 

.00501475 

-.000184366 

X1 

y 

l 

0 

/ 

x2 

(187) 


It  is  clear  that  the  two-by-two  matrix  made  up  of  the  vectors 
T  and  CT  is  a  transformation  matrix  allowing  the  new  output 
vector  to  be  expressed  as  a  linear  combination  of  the  two 
system  states.  If  this  transformation  matrix  has  an  inverse, 
then  the  two  system  states  can  be  expressed  as  a  linear 
combination  of  the  outputs  w  and  y  .  This  matrix  does  indeed 
have  an  inverse,  which  is 


F"1  = 

.00501475 

-.000184366 

-1 

0 

1 

1 

0 

-5424 

27.2 

The  final  system  equations  are  given  by 


(188) 


w  =  -100 w  +  y  -  .49853  (185) 

x1  =  y  (189) 

x2  =  -5424  w  +  27.2  y  (190) 


A  block  diagram  of  the  complete  servo  and  observer  system  is 
shown  in  Figure  60. 


A  total  system  may  now  be  easily  formed  introducing  a 


forward-loop  gain,  Kx,  and  two  feedback  gains  Kj  and  K£.  such 
a  system  is  shown  in  Figure  61.  The  transfer  function  for 
this  system  is  also  given  in  the  figure. 


Figure  60 .  Servo  and  first-order  observer . 


*  * 


Appendix  J 

C*  Control  Derivation 


The  concept  of  C*  is  based  upon  the  active  control  of  a 
linear  combination  of  longitudinal  aircraft  states;  more 
specifically,  a  weighted  combination  of  pitch  rate,  pitch 
acceleration,  and  the  time  rate  of  change  of  angle-of-attack . 
Introduced  in  1966  by  Tobie,  Elliot,  and  Malcoir  (Ref  22),  one 
definition  of  the  C*  performance  criterion  is 

I 

C*  =  KNUQ(q  -  'a)  +  Kqq  +  q  (191) 

where  KN  ,  Kq ,  and  are  gain  constants  used  to  weight  each 
term,  and  UQ  is  the  forward  velocity  of  the  aircraft.  In 
Chapter  IV,  Eqn  28  defines  normal  acceleration,  positive 
downward,  as 

Az  =  M'a  -  9)  (28) 

which  reveals  that  the  first  term  in  Eqn  191  is  actually 
weighted  normal  acceleration,  defined  positive  upward.  Making 
this  substitution,  and  letting  AN  =  -Az,  the  definition  of  C* 
looks  a  bit  more  familiar: 


C*  =  KnAn  +  l^q  +  K^q 


(192) 


The  C*  control  concept  was  conceived  as  a  result  of  a 
desire  to  improve  aircraft  handling  qualities  based  upon 
pilot  sensory  perception.  It  was  hypothesized  that  at  low 
aircraft  velocities,  a  pilot  responds  more  to  pitch  cues  than 
to  normal  acceleration,  as  large  changes  in  pitch  attitude 
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i  ^  >L  *- 


result  in  only  small  changes  in  normal  acceleration.  At 
higher  velocities,  a  pilot  tends  to  sense  normal  acceleration 
cues  more  than  pitch,  as  slight  pitching  results  in  large 
acceleration  changes.  Therefore,  the  idea  of  C*  control  is  to 
blend  pitch  rate  and  pitch  acceleration  with  normal 
acceleration,  thereby  allowing  the  pilot  to  control  the 
dominant  mode. 

The  quantity  C*  is  generally  dimensioned  in  G's,  implying 
that  Kq  and  must  be  chosen  ,in  such  a  way  as  to  permit  each 
term  in  Eqn  192  to  be  dimensioned  in  G's.  As  a  means  of 
standardization,  KN  normally  assumes  a  unit  value,  and  and 
Kq  are  then  selected  on  that  basis.  Making  this  assumption  is 
in  no  way  restrictive,  as  it  is  only  the  relative  weighting 
of  feedback  quantities  that  is  important.  Equation  192  can 
now  be  expressed  as 

C*  =  AN  +  K^q  +  K^q  (193) 

As  it  is  pitching  acceleration  at  the  pilot  station  that 
is  of  interest,  the  value  of  is  selected  based  upon  the 

physical  distance  from  the  aircraft  CG  to  the  location  of  the 
cockpit.  In  the  YF-16,  these  pitching  accelerations  are 
sensed  by  the  normal  accelerometers,  as  they  are  located 
directly  under  the  pilot's  seat.  This  occurance  is  no 
accident,  however,  as  the  YF-16  has  a  somewhat  modified  C* 
flight  control  system  using  washed-out  pitch  rate  to  respond 
to  pitch  transients  only.  For  the  YF-16  then,  the  C* 
definition  is  reduced  to 


c*  =  ANg  +  Kq  q 


(194) 


where  Afls  is  the  combination  of  normal  and  pitching 
accelerations  at  the  normal  accelerometer  location.  This  is 
the  definition  used  in  the  C*  FCS  shown  in  Figure  58  in 
Appendix  F. 

The  quantity  Kq  represents  a  chosen  "crossover  velocity", 
or  the  aircraft  velocity  at  which  normal  acceleration  and 
pitch  rate  will  be  controlled  equally.  Assuming  that  velocity 
is  dimensioned  in  feet/seccnd,  and  pitch  rate  is  measured  in 
degrees/second,  the  dimensional  product  of  these  two 
quantities  will  be  f t-deg/sec-sec.  As  described  in  Chapter 
IV,  Eqn  32,  a  dimensional  conversion  factor  must  be  applied 
to  permit  a  dimensional  product  of  G's.  If  "Vco"  is  defined 
as  the  crossover  velocity  in  feet/second,  Eqn  194  becomes 


C* 


Ans  + 


vco 

1843.6 


(195) 


Typically,  Vco  is  selected  as  400  feet/second  for  jet  fighter 
aircraft  (Ref  22) .  For  the  YF-16  then,  the  final  form  of  the 
C*  definition  is  given  as 

C*  =  ANs  +  0.217  q.  (196) 

The  performance  of  a  C*  commanded  FCS  is  generally 
assumed  to  be  acceptable  if  the  system  transient  response  to 
a  step  input  falls  within  a  specified  envelope,  depending 
upon  the  particular  mode  (category)  of  flight.  There  are  four 


such  categories  defined  in  MIL-F-8785B  (Ref  13),  and  they 
have  been  directly  applied  to  the  transient  performance  of  a 
C*  FCS.  The  envelopes  are  depicted  in  Figure  62,  and  the 
categories  of  flight  are  as  defined  below. 

I  Optimal  response  for  critical  flying  conditions 
(aerial  combat,  ground  attack,  penetration). 

II  Response  for  less  critical  conditions  (cruising, 
air  refueling)  . 

/ 

III  Response  for  conditions  not  covered  by  I,  II,  or 
IV  (loiter)  . 

IV  Response  for  power  approach  (landing) . 
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YFCCV  Computer  Prcgrair 


YFCCV,  a  FORTRAN  Extended  software  package  developed  as  a 
part  of  this  thesis,  and  designed  to  run  on  a  CDC  6600  series 
computer,  was  created  for  the  purpose  of  simulating  the 
multi-input,  multi-output  longitudinal  flight  control  system 
of  the  YF-16  CCV.  Although  the  program  was  designed 
especially  for  this  aircraft,  it  is  capable  of  simulating  the 
longitudinal  dynamics  and  flight  control  system  of  almost  any 
other  aircraft  having  a  maximum  of  two  active  longitudinal 
aerodynamic  control  surfaces.  The  user  is  free  to  specify  the 
dynamics  of  a  normal  accelerometer,  a  pitch  rate  gyro,  an 
angle-of-at tack  sensor,  the  servos  and  actuators  for  the 
horizontal  tail  and  tr ailing-edge  flaperons,  two  forward-loop 
cascade  compensators,  and  up  to  six  feedback  compensators. 

The  aircraft  aerodynamic  data  is  conveniently  input  in 
the  form  of  non-dimensional  stability  derivatives.  The  pitch 
rate  gyro,  angle-of-  at  tack  sensor,  normal  accelerometer,  and 
all  compensators  are  input  as  transfer  functions.  Pilot 
inputs  include  step,  pulse,  ramp,  and  polynomial  functions  of 
up  to  5th  order.  Disturbance  inputs  include  step  and  "1-cos" 
discrete  vertical  gusts,  and  a  uniformly  distributed 
pseudo-random  vertical  gusting.  The  magnitudes  (and  duration, 
where  applicable)  of  both  pilot  and  disturbance  inputs  are 
user  definable.  In  addition,  both  pilot  and  disturbance 
inputs  may  be  used  in  the  same  simulation.  In  this  manner, 
the  user  can  observe  aircraft  response  to  a  pilot  input  in 


the  presence  of  turbulence. 

Program  output  is  available  as  a  listing  or  in  plotted 
form.  In  most  cases,  output  data  will  be  routed  to  a  line 
printer.  Using  this  output  device,  listings  of  measurements 
taken  at  21  points  in  the  flight  control  system  are  printed 
out  at  user  definable  time  increments.  In  addition,  if 
desired,  the  program  will  generate  line  printer  plots  for  the 
same  simulation  to  be  printed  with,  or  instead  of,  the 
listings.  Any  of  the  21  measurements  can  be  plotted  (scaled 
automatically),  but  a  maximum  of  15  may  be  requested  in  any 
one  run.  As  an  alternative,  the  user  may  request  that  the 
measurement  data  be  formatted  for  routing  to  a  Calcomp  (or 
equivalent)  plotter.  The  computer-generated  plots  in  this 
report  were  created  in  this  manner. 

If  desired,  .  the  horizontal  tail  and  trailing-edge  flap 
servos  may  be  rate  limited  and/or  displacement  limited.  When 
the  deflection  of  either  of  these  two  surfaces  is  plotted  as 
a  line  printer  plot,  the  maximum  rate  of  deflection  achieved 
is  calculated  and  printed  out  on  the  upper  right-hand  corner 
of  the  plot.  If  either  limiter  is  in  use,  the  flag  "LIMIT, 
ON"  will  also  be  printed,  saturation  is  specified  in  degrees 
per  second  or  in  degrees. 

The  program  may  be  used  interactively  or  in  batch  mode. 
Having  been  designed  for  interactive  use,  all  data  input  is 
prompted.  However,  a  data  file  may  be  used  if  the  same  input 
format  is  maintained.  In  either  mode,  all  specified 
parameters  are  retained  within  the  program  until  execution  is 


terminated.  This  feature  allows  the  user  to  run  several 
simulations,  perhaps  changing  a  compensator  eigenvalue 
location,  without  having  to  re-enter  the  unchanged  input  data 
each.  time. 

All  simulation  output  data  (except  for  Calcomp  use  as 
mentioned  previously)  is  formatted  for  the  line  printer.  This 
was  done  because  of  the  large  volume  of  output  generally 
produced.  However,  the  user  may  view  this  output  at  an 
intercom  terminal  if  desired./ The  flow  of  simulation  output 
is  controlled  by  a  switch  command.  Simply  typing  "COPY, ON" 
will  cause  all  simulation  output  to  be  written  to  a  local 
file  called  "COPY".  Typing  "COPY, OFF"  will  cause  future 
output  to  be  written  to  the  terminal.  All  output  other  than 
simulation  data  is  still  connected  to  the  terminal  in  either 
switch  position.  For  this  reason,  the  switch  may  be  turned  on 
at  the  beginning  of  an  interactive  session  and  still  allow 
normal  program  operation  with  the  advantage  of  the  user  not 
having  to  watch  data  being  listed.  At  the  conclusion  of  the 
session,  the  COPY  file  may  be  routed  to  the  line  printer  or 
viewed  at  the  intercom  terminal. 

In  addition  to  the  simulation  routine,  there  are  a  few 
other  options  available  to  the  user.  Originally,  these 
options  were  intended  to  be  used  for  program  de-bugging,  but 
were  found  to  be  useful  in  other  ways  and  thus  retained. 
Since  the  program  uses  a  state  space  representation  of  the 
aircraft  in  its  simulation  process,  printout  of  these 
matrices  (in  either  continuous  or  discrete  form)  can  be 


requested.  Given  a  set  of  non-dimensional  stability 
derivatives,  the  A,  B,  F,  and  G  matrices  may  be  obtained  for 
the  dynamics  of  any  aircraft.  Also,  the  F  and  G  matrices  for 
any  general  A  and  B  matrices  (up  to  5th  order)  may  be 
obtained.  In  addition,  the  F  and  G  matrices  may  be  obtained 
for  a  given  transfer  function  of  up  to  5th  order.  The 
remaining  options  include  setting  up  a  state  space 
representation  of  aircraft  dynamics  generated  from  stability 
derivatives,  the  capability  tp  "convert"  a  root  location  in 
s-plane  to  a  corresponding  location  in  z-plane,  and  executing 
a  polynomial  curve  fitting  routine  using  up  to  10th  order 
polynomials  to  fit  data  points. 

Theory  of  Operation 

During  a  simulation,  all  integration  is  performed  using 
the  state  transition  method.  The  aircraft  dynamics,  pitch 
rate  gyro,  angle-of-attack  sensor,  normal  accelerometer,  and 
all  compensator  transfer  functions  are  converted  to  discrete 
state  space  form  (Ref  2) .  With  each  iteration,  the  values  of 
the  states  propagate  throughout  the  entire  system,  as  shown 
in  Figure  63.  Since  the  state  transition  method  performs  an 
exact  integration  of  a  linear  system,  accuracy  is  maintained 
within  each  "subsystem"  (allowing  for  some  round-off  error). 
Some  loss  of  accuracy  occurs,  however,  when  interfacing  these 
subsystems  to  produce  the  entire  system.  As  the  integration 
of  states  is  being  performed  in  separate  parts,  there  is  the 
underlying  assumption  that  the  input  to  each  subsystem 
remains  constant  during  the  integration  process.  Except  in 


Figure  63 •  Basic  functional  YPCCV  Computer  Program  block  diagram  showing  propagation  of 
subsystem  states. 


steady-state,  this  never  truely  happens.  It  is  assumed, 
however,  that  if  the  iteration  rate  is  much  greater  than  the 
maximum  rate  of  change  of  any  system  state,  the  approximation 
is  valid. 

To  test  this  scheme,  the  system  in  Figure  64  was 
simulated  for  a  step  input  using  the  program  at  various 
iteration  rates.  The  same  system  was  simulated  using  an 
inverse  Laplace  technique,  or  actual  time  function  (Ref  23). 
The  results  of  these  simulations  are  plotted  on  one  axis  in 
Figure  65.  Note  that  as  the  iteration  rate  is  increased,  the 
discrete  approximation  rapidly  converges  to  the  continuous 
response.  Although  a  rate  of  200  iterations  per  second  has 
been  found  to  be  suitable  for  simulating  most  aircraft  flight 
control  systems,  a  rate  of  2000  iterations  per  second  was 
used  for  the  majority  of  the  simulations  in  this  thesis. 

There  are  two  main  reasons  for  chosing  this  integration 
scheme.  The  first  is  that  it  allows  the  simultaneous 
application  of  several  inputs  to  a  system  in  a 
straightforward  manner  (matrix  operations).  The  second  is 
that  it  allows  each  subsystem  to  be  changed  independently  of 
the  others,  and  avoids  the  alternative  of  reconstructing  and 
processing  a  system  matrix  that  may  contain  up  to  59  states 
(this  is  the  maximum  number  of  states  available  in  YFCCV).  In 
addition,  this  method  allows  the  simulation  of  certain 
digital  controllers.  Since  the  FCS  signals  are  applied  to  the 
aircraft  at  each  iteration,  and  held  constant  between 
iterations,  the  program  operates  using  the  zero-order  hold 
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concept.  By  allowing  the  servo/actuator,  aircraft,  and  gyro 
subsystems  (having  continuous  dynamics)  to  operate  at  a 
faster  iteration  rate,  and  the  feedback  to  operate  at  a 
slower  rate  (a  sampling  rate),  the  program  is  capable  of 
simulating  a  digital  controller. 

Commands 

After  the  program  has  been  attached  for  interactive  use, 
simply  entering  "YFCCV"  will  begin  execution,  and  result  in 

i 

the  prompt  " <COMMAND> " .  At  this  point,  the  user  may  enter  a 
command,  a  "?"  to  obtain  a  listing  of  valid  commands,  or 
"STOP"  to  terminate  the  program.  At  any  time  that  the  program 
requests  an  alphabetical  response,  entering  a  "?"  will 
provide  an  explanation  of  the  requested  response.  In  addition 
to  '*?"  and  "STOP",  there  are  21  other  valid  commands.  These 
commands  are  listed  on  the  next  several  pages  with 
explanations  of  their  functions. 

(1)  AC DAT  Using  this  option,  the  user  specifies 
certain  physical  properties  of  the  aircraft,  defines  the 
flight  condition,  and  inputs  the  non-dimensional  stability 
derivatives  for  that  flight  condition.  For  the  YF-16,  some 
quantities  are  defined  internally  (wing  area,  aircraft  mass, 
MAC,  etc),  and  in  addition,  if  a  flight  condition  of  0.8  Mach 
at  SL  is  specified,  the  required  stability  derivatives  are 
specified  internally.  The  location  of  the  normal 
accelerometer  ("AZL",  the  distance  in  feet  from  the  CG  to  the 
accelerometer  location)  must  also  be  specified  in  this 
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option.  The  input  of  this  data  if-  necessary  prior  to  any 
simulation  attempt. 


(2)  ACSTM  This  command  will  cause  the  aircraft  state 
transition  matrix  and  input  matrices  to  be  calculated  and 
printed  out.  Before  requesting  this  option,  the  user  must 
enter  aircraft  data  and  specify  a  flight  condition  using 
ACbAT. 

(3)  ALPHA  Using  this  command,  the  user  specifies  the 
feedback  from  the  attack  angle  sensor  to  the  HT  (horizontal 
tail)  and  TEF  ( t r a i 1 i ng-edge  flaperon)  channels.  Each 
feedback  may  be  zero  (default),  a  gain,  or  a  compensator  (up 
to  5th  order).  To  accommodate  a  wider  variety  of  system 
configurations,  cither  attack  angle  or  its  derivative  nay  be 
specified  for  feedback  by  entering  the  responses  "NCFMAL"  or 
"DLK1VAT1VE"  ,  respectively,  M  the  requested  time.  Input  of 
this  feedback  element  is  for  use  in  GIN SIM  simulations  only. 

(4)  AN  Using  this  command,  the  user  specifies  the 
feedback  from  the  normal  accelerometer  to  the  HT  and  TEF 
channels.  Each  feedback  may  be  zero  (default),  a  gain,  or  a 
compensator  (up  to  5th  order).  It  should  be  noted  that  the 
position  of  the  normal  accelerometer  is  specified  by  data 
input  using  the  ACDAT  command.  Input  of  this  feedback  element 
is  for  use  in  GENS1M  simulations  only. 

(5)  COPY ,0N  or  COPY , OFF  The  COPY  switch  command 
controls-  the  file  to  which  simulation  output  data  is  written. 
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Turning  this  switch  on  will  cause  simulation  data  to  be 
written  to  the  local  file  "COPY".  Turning  it  off  will  cause 
the  data  to  be  written  to  the  "OUTPUT"  file  connected  at  the 
intercom  terminal  (default) .  The  status  of  this  switch  does 
not  affect  any  other  program  operation  or  user  prompting. 

(6)  CSA  This  command  is  an  acronym  that  stands  for 

Compens a tor/Ser vo/Actuator.  The  combined  dynamics  of  these 

devices  are  entered  as  a  single  transfer  function  (one  each 

/ 

for  HT  and  TEF  channels)  using  this  option.  Input  of  this 
feedforward  element  is  for  use  in  GENSIM  simulations  only. 

(7)  CURVE  This  command  executes  a  polynomial  curve 

fitting  routine.  It  is  limited  to  10th  order  polynomials. 
Calcomp  plots  and  numerical  listings  are  available  to  display 
both  input  data  points  and  the  result  of  the  curve  fitting 
process. 

(8)  GENSIM  This  command  is  used  to  initiate  a 
simulation  process  that  is  referred  to  as  "GENeral 
SIMulation".  This  is  the  standard  simulation  of  a 
continuous-time  flight  control  system  of  the  general  form 
depicted  in  Figure  66.  The  user  specifies  the  desired 
iteration  rate  (0.0005  is  typical),  the  form  of  the  output 
desired,  and  the  location  at  which  normal  acceleration, 
velocity,  and/or  position  will  be  measured  (if  measurement  is 
requested)  during  the  simulation.  Although  the  normal 
acceleration  used  for  feedback  within  the  FCS  is  always  from 
the  specified  accelerometer  location  (using  ACDAT) ,  this 
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state  and  its  integrals  (normal  velocity  and  position)  may  be 
measured  at  either  the  CG  or  the  accelerometer  location.  This 
permits  observation  of  these  quantities  at  one  location, 
while  the  system  feedback  may  represent  another  location.  The 
actual  simulation  begins  at  the  close  of  the  command  sequence 
immediately  after  the  requested  command  "PUN"  is  entered. 

(9)  GSTM  With  this  option,  the  user  may  input  any 

continuous  state  space  system  (matrix  form)  of  up  to  5th 

/ 

order,  and  obtain  a  state  transition  matrix  and  input  matrix 
for  a  specified  sampling  rate.  In  addition,  a  time  simulation 
of  the  system  can  be  obtained  for  any  number  of  iterations, 
with  any  state  as  output.  This  output  is  available  in  listed 
form  only  (formatted  for  intercom  display). 

(10)  GYRO  This  command  allows  the  dynamics  of  the 
pitch  rate  gyro,  the  angle-of-attack  sensor,  and  the  normal 
accelerometer  to  be  specified.  A  minimum  of  a  first-order  lag 
must  be  specified  for  each  of  these  devices  that  will  be  used 
in  the  simulation,  or  zero  output  (default)  will  result  from 
that  device.  Any  transfer  function  of  up  to  5th  order  is 
permissible. 

(11)  This  command  allows  the  user  to  specify 
initial  conditions  on  (the  perturbation  quantities)  forward 
velocity,  pitch  angle,  pitch  rate,  and  angle-of-attack.  This 
option  is  often  useful  when  designing  and  simulating  a 
regulator  without  pilot  or  disturbance  inputs.  Initial 
conditions  can  be  specified  only  for  use  in  GENSIM 
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simulations 


(12)  INPUT  This  command  permits  definition  of  pilot 

and  disturbance  inputs.  If  either  of  these  two  inputs  is  not 
defined,  default  is  zero.  Cither  one  or  both  inputs  may  be 
specified . 

(13)  OPTSIM  This  command  is  used  to  initiate  a 

simulation  process  that  is  referred  to  as  "OPTimal 
SIMulation".  This  option  ,  provides  a  more  specialized 
simulation  capability  (relative  to  GENSIM),  and  is  intended 
for  use  in  simulating  an  optimal  servo  class  of  flight 
control  system  configurations  using  gain  feedback  and 

forward-loop  integration  (Pef  11).  The  user  may  specify 
either  a  continuous-time  or  discrete-time  controller.  In  the 
latter  case,  controller  sample  rate  must  be  provided.  This 
option  further  differs  from  GENSIM  in  that  system  parameters 
are  specified  in  matrix  format,  rather  than  as  individual 
feedback  and  feedforward  compensators.  As  with  all  YFCCV 
commands,  input  requests  are  fully  prompted. 

(14)  PLOTTER  This  command  permits  the  user  to  plot 

simulation  results  on  a  Calccmp  (or  equivalent)  plotter.  A 
simulation  must  first  be  run  (using  GENSIM  or  OPTSIM)  to 
calculate  all  the  required  data  points.  By  entering  the 
response  "TAPE”  when  prompted  by  the  question,  "DO  YOU  WANT 
ANY  LINE  PRINTER  PLOTS?",  this  information  is  written  into 
the  local  file  "COPY"  by  the  simulation  routine,  in  a  format 
that  enables  the  PLOTTER  routine  to  extract  the  data  as 


required.  The  response,  "TAPE",  is  entered  during  the 
simulation  command  sequence,  and  simply  indicates  to  the 
routine  that  simulation  output  should  be  written  to  the  file 
COPY  for  use  by  PLOTTER.  Each  plot  requested  causes  the 
program  to  read  the  COPY  file,  reformat  and  insert  the 
necessary  labeling,  and  write  into  an  additional  local  file 
called  "PLOT".  This  PLOT  file  must  then  be  routed  to  the 
particular  plotter  that  is  to  be  used. 

/ 

(15)  SERVOFB  Using  this  command,  the  user  specifies 
four  feedback  gains  from  the  outputs  of  the  two  CSA 's . 
Essentially,  this  option  permits  a  feedback  loop  around  each 
CSA,  and  feedback  from  the  output  of  each  CSA  to  the  input  of 
the  alternate  CSA.  The  feedback  is  summed  negatively  at  the 
CSA  inputs,  and  therefore  these  gains  will  be  positive  for 
negative  feedback.  No  compensation  is  available  in  these 
feedback  paths.  As  the  concept  of  a  CSA  is  valid  only  for  use 
in  GENS1M  simulations,  this  feedback  is  otherwise  undefined. 

(16)  STSPACE  This  command  will  cause  the  continuous 
state  space  representation  of  the  aircraft  with  two 
first-order  servos  to  be  printed.  Also,  for  a  more 
specialized  problem,  this  option  will  calculate  a  state  cost 
matrix  from  an  output  cost  matrix  for  optimal  control 
problems. 

(17)  TARGET  This  option  permits  the  user  to  specify 
the  distance  (in  feet)  from  the  aircraft  to  a  target 
(defaults  to  6000  feet).  This  information  is  used  in  the 
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calculation  and  printout  or  plot  (if  requested)  of 
longitudinal  target  tracking  error  as  defined  in  Chapter  II. 

(18)  THDOT  Using  this  command,  the  user  specifies  the 
feedback  from  the  pitch  rate  gyro  to  the  HT  and  TEF  channels. 
Each  feedback  may  be  zero  (default),  a  gain,  or  a  compensator 
(up  to  5th  order).  To  accommodate  a  wider  variety  of  system 
configurations,  either  pitch  rate  or  its  derivative  may  be 

specified  for  feedback  by  entering  the  responses  "NORMAL"  or 

/ 

"DERIVATIVE",  respectively,  at  the  requested  time.  Input  of 
this  feedback  element  is  for  use  in  GENSIM  simulations  only. 

(19)  TIME  Use  of  this  command  results  in  the  current 

time  being  printed  at  the  user's  terminal. 

(20)  TRF  This  option  allows  the  user  to  input  any 

realizable  transfer  function  of  up  to  5th  order,  and  obtain  a 
state  transition  matrix  and  input  matrix  for  a  given  sampling 
rate  (F  and  G  matrices).  In  addition,  a  time  simulation  of 
the  transfer  function  may  be  obtained. 

(21)  ZROOT  This  option  performs  an  s-plane  to  z-plane 

root  location  "conversion".  Any  single  root  or  complex  pair 
of  roots  represented  in  s-plane  can  be  input,  and  the 
corresponding  root  location (s)  in  z-plane  will  be  obtained 
for  a  specified  sampling  rate. 

Any  of  the  above  commands  may  be  entered  at  the  prompt 
”<COMMAND>”.  Whenever  all  requested  data  has  been  input,  and 


any  requested  output  has  been  generated  in  a  given  command 
sequence,  the  user  is  returned  to  <COMMAND>  mode.  At  any 
point  in  the  program  where  an  alphabetical  response  is 
requested,  the  entry  of  an  invalid  command  will  result  in  the 
response,  "WHAT??",  or  "PARDON  ME??",  at  which  time  the 
proper  entry  can  be  re-typed. 

The  next  several  pages  show  a  sample  interactive  session 

with  the  program.  During  the  session,  several  of  the  commands 

/ 

required  for  GENSIM  simulations  are  utilized.  The  end  result 
would  be  an  actual  simulation  of  an  aircraft  flight  control 
system.  In  this  particular  example,  no  stability  derivatives 
are  entered,  as  all  required  aircraft  data  for  the  YF-16  at  a 
FC  of  0.8  Mach  at  SL  is  permanently  stored  within  the 
program.  Angle-of-attack ,  pitch  rate,  and  normal  acceleration 
feedback  elements  are  defined,  as  are  the  forward-locp 
dynamics  (CSA's).  The  transfer  functions  of  the 
angle-of-attack  sensor,  pitch  rate  gyro,  and  normal 
accelerometer  are  then  input,  followed  by  the  GENSIM 
simulation  command  sequence.  The  requested  output  is  four 
Calcomp  plots. 

Unfortunately,  due  to  the  length  of  program  YFCCV,  it  is 
not  practical  to  include  a  program  listing  in  this  thesis. 
Upon  request,  however,  a  listing  is  available  from  the 
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COMMAND-  YFCCV 


TYPE  '?'  FOR  <COMMAND>  OPTIONS,  'STOP'  TO  STOP... 

<COMMAND>  ACDAT 

AIRCRAFT  DESIGNATION:  YF-16 

ENTER  INITIAL  ALTITUDE , MACH  NO.:  0  .8 
SERVO  DISPLACEMENT  LIMITERS  (ON  OR  OFF):  OFF 
SERVO  RATE  LIMITERS  (ON  OR  OFF)  :  ON 
ENTER  HT  SERVO  LIMITING  RATE  ((DEG/SFC):  6>£ 

ENTER  TEF  SERVO  LIMITING  RATE  (DEG/SEC):  _56_ 

ENTER  TYPE  OF  ANALYSIS  (LINEAR  OR  NON-LINEAR):  LINEAR 
DO  YOU  WISH  TO  INPUT  STABILITY  DERIVATIVES?  NO 

<COMMAND>  INPUT 
>>>WIND 

WIND  INPUT  =  1-COS 

ENTER  VERTICAL  GUST  VELOCITY  (FT/SEC),  DURATION:  50  1 

>  >>EXIT 

< COMMA ND>  ALPHA 

SELECT  'STATE'  OR  'SENSOR'  ALPHA  FEEDBACK:  SENSOR 
>>>HT 

FEEDBACK  STATE  (NORMAL  OR  DERIVATIVE):  DEPIV 
FEEDBACK  SPECIFICATION  = COMPENSATOR 

ENTER  ORDERS  OF  NUMERATOR, DENOMINATOR  (MAX  5TH ) :  1  1 

ENTER  NUMERATOR  CONSTANT:  .035 

ENTER  2  NUMERATOR  COEFFICIENTS .. .HI  TO  LOW: 

1  8 

ENTER  2  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  5 

>>>TEF 

FEEDBACK  STATE  (NORMAL  OR  DERIVATIVE):  DEPIV 
FEEDBACK  SPECIFICATION  =COMP 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH ) :  1  1 

ENTER  NUMERATOR  CONSTANT:  -1.875 

ENTER  2  NUMERATOR  COEFFICI ENTS . . .HI  TO  LOW: 

1  B 

ENTER  2  DENOM  COEF FICIENTS . . .HI  TO  LOW: 

1  0 

>  »EXIT 
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<CCMMAND>  THDOT 

>  >>HT 

FEEDBACK  STATE  (NORMAL  OR  DERIVATIVE):  DEPIV 
FEEDBACK  SPECIFICATION  =COMP 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH ) :  2  2 

ENTER  NUMERATOR  CONSTANT:  _1_ 

ENTER  3  NUMERATOR  COEFFICIENTS _ HI  TO  LOW: 

.075  1.225  5 

ENTER  3  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

15  0 

>  »TEF 

FEEDBACK  STATE  (NORMAL  OR  DERIVATIVE):  DERIV 
FEEDBACK  SPECIFICATION  *  COMP 

ENTER  ORDERS  OF  NUMERATOR, DENOMINATOR  (MAX  5TH ) :  1  1 

ENTER  NUMERATOR  CONSTANT:  -.875 

ENTER  2  NUMERATOR  COEFFICIENTS. . .HI  TO  LOW: 

1 _ 8 

ENTER  2  DENOM  COEFFICIENTS ..  .HI  TO  LOW: 

1 _ 0 

>  >>EXIT 

<COMMAND>  AN 

>  >>TEF 

FEEDBACK  SPECIFICATION  =  COMPENSATOR 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH) :  1  1 

ENTER  NUMERATOR  CONSTANT:  187.5 

ENTER  2  NUMERATOR  COEFFICIENTS .. .HI  TO  LOW: 

1  8 

ENTER  2  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  100 
>>>EXIT 

<COMMAND>  CSA 

>  >>HT 

ENTER  ORDERS  OF  NUMERATOR  DENOMINATOR  (MAX  5TH)  :  1  3 

ENTER  NUMERATOR  COEFFICIENT:  -160 

ENTER  2  NUMERATOR  COEFFICIENTS. . .HI  TO  LOW: 

1  5 

ENTER  4  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  28  160  0 
>>>.TEF 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH) :  1  3 

ENTER  NUMERATOR  COEFFICIENT:  -320 

ENTER  2  NUMERATOR  COEFFICIENTS .. .HI  TO  LOW: 

1  0 

ENTER  4  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  43  580  2400 
>»EXIT 
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<COMMAND>  GYRO 
>»ALPHA 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH) :  0  1 

ENTER  NUMERATOR  CONSTANT:  66.7 

ENTER  1  NUMERATOR  COEFFICI ENTS . . .HI  TO  LOW: 

1. 

ENTER  2  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  66.7 

>>>THDOT  , 

ENTER  ORDERS  OF  NUMERATOR , DENOMINATOR  (MAX  5TH ) :  0  2 
ENTER  NUMERATOR  CONSTANT:  JL 

ENTER  1  NUMERATOR  COEFFICIENTS .. .HI  TO  LOW: 

98596 

ENTER  3  DENOM  COEFFICIENTS .. .HI  TO  LOW: 

1  376.8  98596 
>>>AN 

ENTER  ORDERS  OF  NUMERATOR, DENOMINATOR  (MAX  5TH ) :  0  2 

ENTER  NUMERATOR  CONSTANT:  98596 

ENTER  1  NUMERATOR  COEFFICIENTS .. .HI  TO  LOW: 

1_ 

ENTER  3  DENOM  COEFFICI ENTS .. .HI  TO  LOW: 

1  376.8  98596 
>>>EXIT 

<COMMAND>  GENSIM 

ENTER  SYSTEM  ITERATION  TIME.  IT=.0005 

ENTER  PRINTOUT  TIME  INCREMENT,  SIMULATION  DURATION:  .005  4 

DO  YOU  WANT  ANY  LINE  PRINTER  PLOTS?  TAPE 

ENTER  CROSSOVER  VELOCITY  (FT/SEC).  VCO=  400 

NORMAL  POSITION,  VELOCITY,  AND  ACCELERATION  CAN  BE 

MEASURED  AT  THE  CENTER  OF  GRAVITY  OR  AT  THE  PILOT 

STATION.  ENTER  'CG'  OR  'PS'  (DEFAULT,  PS):  PS 

TYPE  'RUN'  TO  SIMULATE,  'EXIT'  TO  RETURN  TO  <COMMAND> :  RUN 


<COMMAND>  PLOTTER 

ENTER  NUMBER  OF  PLOTS  DESIRED:  4_ 

ENTER  4  PLOT  IDENTIFIERS:  2456 


***####*###*###**###  DRAWZ  CALLED  **#***##****#####*## 


THE  NUMBER  OF  POINTS  TO  BE  PLOTTED  =  8  00 


#*#*#***###******#**  DRAWZ  CALLED  ##*###**###*#*###### 


THE  NUMBER  OF  POINTS  TO  BE  PLOTTED  =  8  00 


**####**#*##*#*##*##  DRAWZ  CALLED  ####«*<#####<<####*# 


THE  NUMBER  OF  POINTS  TO  BE  PLOTTED  =  800 


I#######**##*#*##*#*  DRAWZ  CALLED  **###*###**###*#*«#* 


THE  NUMBER  OF  POINTS  TO  BE  PLOTTED  =  800 


<COMMAND>  STOP 


GOOD-BYE.  THE  TIME  IS  13.34.36 
STOP  NORMAL  TERMINATION 
5.201  CP  SECONDS  EXECUTION  TIME 


COMMAND-  ROUTE .PLOT, DC=PT ,TIDCBB, FID=YA037 .ST-CSB 
COMMAND-  LOGOUT 


CPA  34.053  SEC.  14  .782  ADJ. 

IO  34.861  SEC.  17.440  ADJ. 

CRUS  43.695 

CONNECT  TIME  0  HRS.  16  MIN. 
01/04/79  LOGGED  OUT  AT  13.36.54 
< 
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Appendix  L 

YFSOLVE  Computer  Proqrair. 

YFSOLVE  is  a  special  purpose  FOETFAN  Extended  computer 
program  designed  to  run  on  a  CDC  6600  series  computer. 
Intended  for  interactive  use,  it  was  developed  for  the 
purpose  of  providing  certain  closed-loop  transfer  functions 
for  the  Pitch  Performance  Enhancement  (PPE)  FCS  from 
specified  input  parameters.  In  addition,  as  the  simulation 
program,  YFCCV  (see  Appendix' K),  is  not  directly  suited  to 
simulating  system  configurations  such  as  the  PPE  FCS,  program 
YFSOLVE  provides  the  required  feedback  and  feedforward  (CSA) 
dynamics  for  an  equivalent  system  model  that  fits  the  YFCCV 
system  configuration.  The  use  of  an  equivalent  model  is  due 
to  the  PPE  FCS  requirement  for  summation  of  feedback  twice  in 
the  horizontal  tail  channel,  a  forward-loop  gain  occurring 
prior  to  feedback  summation  in  the  flaperon  channel,  and  the 
inclusion  of  a  command  prefilter. 

Although  the  program  was  written  to  aid  in  the  design  of 
the  PPE  flight  control  system  for  the  selected  flight 
condition,  minor  alterations  could  enable  it  to  be  used  to 
provide  similar  transfer  functions  for  a  variety  of  systems. 
The  basic  architecture  of  the  program  would  remain  unchanged, 
with  only  the  statements  defining  polynomial  coefficients 
requiring  modification. 

Commands 


Initial  execution  of  the  program  places  the  user  in 


command  mode,  at  which  time  one  of  19  valid  commands  can  be 
entered.  Entry  of  an  invalid  command  results  in  a  response 
of,  "PARDON  ME??".  Of  the  valid  commands,  15  are  used  for 
data  entry,  three  produce  program  output,  and  the  remaining 
command,  "STOP",  is  used  to  terminate  execution.  All  system 
parameters  must  be  specified  (using  the  data  input  commands) 
prior  to  executing  a  data  output  command.  If  all  input 
parameters  are  not  so  specified,  attempted  execution  of  a 
data  output  command  will  result  in  a  response  of,  "HAVE  YOU 
INPUT  ALL  VARIABLES?",  after  which  the  user  is  returned  to 
command  mode.  Valid  commands,  and  their  functions,  are  as 
follows  (given  in  alphabetical  order): 

Data  Entry  Commands 

(1)  FILTER  This  command  permits  definition  of  the 

pilot  command  prefilter  first-order  lag.  As  this  compensator 
is  unity  gain,  specification  of  the  eigenvalue  location 
automatically  sets  the  gain  constant. 

(2)  HPOLE  This  command  permits  specification  of  the 

HT  (horizontal  tail)  first-order  lag  servo  eigenvalue. 

(3)  INPUT  This  command  permits  the  user  to  input  all 

required  system  parameters  (variables)  at  the  same  time. 
Generally,  this  option  would  be  used  at  the  start  of  a 
session,  or  when  a  large  number  of  variables  is  to  be 
modified . 

(4)  Kl_  This  command  permits  specification  of  the 

attack  angle  feedback  gain  to  the  HT  channel. 

(5)  K2  This  command  permits  specification  of  the 


attack  angle  feedback  gain  to  the  TEF  channel. 


(6)  JO  This  command  permits  specification  of  the 

pitch  rate  feedback  gain  to  the  HT  channel. 

(7)  j(4  This  command  permits  specification  of  the 

pitch  rate  feedback  gain  to  the  TEF  channel. 

(8)  1(6  This  command  permits  specification  of  the 

normal  acceleration  feedback  gain  to  the  TEF  channel. 

(9)  KC  This  command  permits  specification  of  the 

command  forward-loop  gain  in  tfie  HT  channel. 

(10)  J(F  This  command  permits  specification  of  the 

TEF  forward-loop  gain. 

(11)  j(JJ  This  command  permits  specification  of  the 

HT  forward-loop  gain. 

(12)  KP  This  command  permits  specification  of  the 

command  forward-loop  gain  in  the  TEF  channel. 

(13)  POLE  This  command  permits  specification  of  the 

washout  circuit  eigenvalue  in  the  TEF  channel. 

(14)  TPOLE  This  command  permits  specification  of  the 

TEF  first-order  lag  servo  eigenvalue. 

(15)  ZERO  This  command  permits  specification  of  the 

zero  location  in  the  proportional  plus  integral  (P+I) 
forward-loop  compensation  in  the  HT  channel. 

Data  Output  Commands 

(16)  L 1ST  This  command  lists  the  current  values  of 

the  system  input  variables. 

(17)  FUN  This  command  generates  all  available 

4 

output  from  the  routine.  The  current  values  of  the  system 


4 


input  variables  are  listed  first,  followed  by  a  listing  of 
the  two  pitch  rate  numerator  equations  and  characteristic 
equation  (polynomial  and  factored  form),  followed  by  a 
listing  of  the  feedback  and  feedforward  gains  and  transfer 
functions  to  be  used  by  program  YFCCV  for  simulation. 

(18)  SIM  This  command  calculates  and  lists  the 

feedback  and  feedforward  gains  and  transfer  functions  to  be 
used  by  program  YFCCV  for  simulation. 

As  the  intended  purpose  of  YF SOLVE  is  to  support  an 
iterative  design  process,  it  was  assumed  that  considerable 
program  output  could  result  while  zeroing  in  on  values  for 
the  system  variables.  To  enable  a  permanent  record  of  this 
output,  the  program  writes  output  to  both  the  user's  intercom 
terminal  and  a  local  file  named  "STUFF".  At  the  end  of  a 
session,  this  file  is  available  to  be  routed  to  a  line 
printer  for  a  "hard"  copy. 

A  listing  of  the  source  code  for  the  program  is  provided 
on  the  next  several  pages. 
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L 


t.vmWu»„vi 


o  o  o  o 


•»  -  »  f.  * 

SUPRrUTIK'E  FACTO* 

PFAL  ROOTR(IO)  ,R33TI ( lr  )  ,CPOLY ( ID) 
COHHCN/EXPCOM/ROOT  R,ROnTI,CPOLY, M: 
K'4  =  C 
I  =  NF+1 

IF  (  fPOL  Y  ( I )  )  3,2,  3 
!  N4=N<  *1 

FOOT  F  <R4)  =1.  0 
root:  ( K'4 >  =o,o 
1  =  1-2 

IF  <  f - -MF )  1,20,2 
I  CONTINUE  ' 

AXF  =  r  ,  c 
A  X 1  =  I  .1 
1=1 
N3  =  l 

ALPlt  =AXR 
A  LC1 1  =  AXI 
H=1 

CO  Tf  26 
PE7 1 r  =T£MR 
°E711=TEMI 
AXF  =  l  .£5 
Al_P2F=AXR 
ALP2]=AXI 


GO  TO  26 
PET2F  =  T  EMR 
pET21=TEMl 
A  XP= (  ,  F 
ALP3F  = AXR  ' 

ALF33=AXI 

R=3 

GO  Tf  25 

°ET  3  F  =T  EHR 

RET  3 1 =TEHI 

TE1=ALP1X-ALP3P 

TE2=ALP1I-ALP3I 

TES-  =  ALP3R-ALP2R 

TEb=f LP3I-ALP2I 

TEH=T  El*T  E5*TEtM  T  E  5 

TE3=  (TE1-  TE5+TE21  T  E5)  /  TEH 

TE'  =  (T  E2*  T  Ebi-TEl*  TE5)  /T  EH 

TE7=T E3*1.0 

TEF  =  I E3* I ES-TEs*  T  E  ! 

TE1G  =  2«0‘  TE3*‘TE4 

0E14  =  TE7 ■  BET  3R-T  E+  '  TFT ’I 


9 


1  0 

1 1 
1  ? 
1  * 


1 


1? 


1C 

17 

1  9 


PEl6  =  TF7‘'nET3I»T?  4  MrT  3R 
T E 1 1  =  T  E3*  B FT  2R -T F  4  M FT  ?1  ♦  BET  1 R-1 E  l  T 
TF1?-TE3  PET2I*Tr4MET?R*B  FT  11-1:15 
T  F7  =T  E9  -  1  .0 

TE1  =  T E9*BET2R-TE1  OMFT  PI 

TE  ?=T  E9*BET2I*TE10MFT  2R 

TE13  =  Tfl-f»ETl?-TF7MFT7R4-TE10*B!TTI 

7El4=TE2-BETlI-TEA  *  TET7i-TFlO«  9 F T  ?  R 

TE16  =  PE1:.>#TE3-1E15*  T  r  A 

TElb  =  PEl‘-*TE‘,+PEl:>,TF7 

lEl=TE13-TE13-T"l!<*TEl'-4.0MTEll*r:i5-TEl?,TE16) 
TE2=2.0*TE13«TE1‘  -A.0MTE12*TEl~*rEll*TEl:>* 

TFM=  5  PET  (  TE1*  TE1+TE2MP2) 

TF  ( 1 E  1 )  9,9,10 

TF i  =5  ORT ( 0.5M  TE^-TEl) > 

TF3=t ,r*7 E2/TE* 

GO  TG  l7 

TE3=SCRT  <  0.5*  (TEMTEl )  ) 

i c  nzz)  11,12,1'' 

TE3=-TE3 
TE*-.  =  t  .‘  "TE2/TE3 
T  E?  =1  E 13  +  T  E3 
TEt=,,E14  +  TFl* 

T  E*r  =  T  El  3-TE3 
TE 10  =  T FI-'  -TE‘ 

TEl  =  r.U*T  Els 
TE?=2 »  0*T Fib 

TF  (TEV*TE/  +  TE8’,lE9-TE°*TE9-TE10,rElC)  14,14,1 
TE7  =  TFC 
T  Eo  =  IE  10 

TEK  =  1 F7*T  E74TE6-1 E 3 
7Eo=  (T  El'  TE7+TF?'  TE3)  /T  EM 
TE-,  =  (TF?'  TF7  -TE1’  T  E9  )  /  TEM 
flXK=^LF3R*TE3*TZ:  -TE4»T  Eb 
fiXl  =  t LT3] +TE3*TE:  ♦TE4*~E5 
flLF4K=eXR 


A  L  p4  1  =  A  X I 

M  =  *^ 

r-0  TC  26 

IF  (APStHELL) ♦ABF(  BELL) -1.0E-20)  19,19,1^ 

TE7  =  ABS(  ALp3F-AXp)  +ABS  (ALP3I-AXH 
TF  UEV/(ABS  (AXR)4A3S  < A  X  I )  >  -  1 . 0 r- ')  19,19,16 


M  =  N?M 
AlFH=ALP2R 
ALF1TsAL°2I 
ALP2F  =  ALD?R 
ALP2T=ALP3I 
ALP7' =ALD4R 
ALF3]=ALp4I 
BET  1P  =  PFT  ?F 
nE7 1 1  =  RET  21 
BET  24  =  R FT  3R 
nc  T21  =  PET3I 


26 1 


PEI  3F  sTEMP. 

PET  3  3  =TEMI 
IF  (f.3-iu0)  t>,19,19 
11  K'4=N<  +1 

ROOT  F  (  NU  )  =  A|_P4R 
F OOTJ <N4) =  ALP4I 
N3  =  0 

IF  U^-NF)  21,20,20 

20  RETURN 

21  Ip  (APS(R00TI(N4> ) -1 . 0r-P>  4,4,2? 

22  IF (L.E0.2)  GO  TO  * 

2  3  AXk=ALPlK 

AX1=-ALD1I 
ALP1 J=-ALP1I 
«*'• 

GO  TO  26 

24  PET 1K=TEMR 
PET 1 1=T  Ell 
AXF=ALP2R 
AX]=-ALP2I 
ALF23  =  -ALP2I 

M  =  t 

GO  TO  2o 

25  PET2F  =  TE  1R 
PET2I=TERI 
AXP= ALP3R 
AXi= -ALP3I 
ALF31=-ALP3I 
L=  2 

H=2 

25  T E PK  =  C- POL  Y  ( 1 ) 

TEMI  =  C  •  0 
DO  2'  1=1, NF 
TElrTEPR^  AXR-T  E9 1 ♦  AX  I 
TEMI  =  TEMI*  AXR+TEMRMXI 

27  TEkR=TEl*CPOLy (1+1) 

«ELL=TEMR 

RFLL=TEMI 

Ip  (»4>  28,30,23 

28  DO  25  1=1, N4 
TEt'l=AXR-RCOTR(I> 

TEM2=AXI-RC0TI (I) 

T  E 1  =  T EMl  'TEMl  +  TE“2»r EM  2 
T  E2=  (TERR*TEMi*Tr3TMEM2)/TEl 
TF^  I-(TEMI'TE'11-TEM?*TpM2)^TEl 

29  TE»‘R=TF2 

30  GO  TO  (5  ,G,7 ,1E  ,  2,  ,25)  ,  M 
END 
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